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ON THE APPLICATION OF LYAPUNOV'S DIRECT METHOD 
TO DIFFERENTIAL EQUATIONS WITH 
AMBIGUOUS RIGHT SIDES* 


Yu. I. Alimov 


(Sverdlovsk) 


Translated from Avtomatika i Telemekhanika, Vol. 22, No. 7, 
pp. 817-830, July, 1961 
Original article submitted November 26, 1960 


Physical considerations are presented by virtue of which the dynamical equations of relay control 
systems may be considered as contingent equations [1,2] with ambiguous and continuous, in a certain 


sense, right sides, The theory of contingent differential equations is used to apply the direct method 
of Lyapunov to “relay” differential equations. 


1.The dynamical equations of a relay control system, 


x= f(t, x), (1) 


where x and f are n-dimensional column matrices [3] are as a rule considered as differential equations with discon- 
tinuous right sides, This point of view does not always correspond to the physical sense of the problem and leads to a 
very sharp line of demarcation between classical and “relay” differential equations as objects of mathematical inves- 
tigation, Another approach, completely justified physically, is possible in which one bases the study of relay control 
systems on the strict theory of contingent equations; this theory has been developed sufficiently well by Zaremba [1) 
and Marchaud [2], The following considerations, repeatedly put forth by many authors in a manner suitable to the 
particular case of relay systems, may serve as our guiding principle. 


In composing the equations of motion of a relay system, each relay link is often considered as some union of 
“continuous” links, described ordinarily by differential equations (possibly with retarded arguments or partial deriva- 
tives, although such cases will not be considered) and, in general, by certain inertialess relay elements with a single 
input and a single output. Then the equations of motion, being put into normal form, may be written in the form of 
a system: 


x =f (t, x) =f [t, x, 6 = o(t, x), (2) 

where f(t, x, g) and o(t,xX) are continuous functions, and ¢o)is the characteristic of the relay element. For con- 
venience in exposition, we shall meanwhile limit ourselves to the case of a single relay element, Here and in what 
follows, the continuity of functions is understood as continuity with respect to all their explicitly described arguments, 


In any working regime, both the input o = o(t) and the output ¢ = ¢fo(t)] coordinate of each real link approxi- 
mated by an inertialess relay element are continuous functions of time. Without going into a detailed description of 
the form of the characteristic g(o) near the so-called switching thresholds 


a(t, x) = Gr (k = ky, ke, (3) 


the function go) may be considered discontinuous in (3) only in regimes accompanied by sufficiently rapid traversal 
of the coordinate o(t) through the switching surfaces (3), when the magnitude of ¢(t) = ¢[.o(t)] changes exceedingly 
rapidly near (3) in comparison with other regions of the t,x space (the meaning of inexact terminology such as “suf- 


* The contents of this paper were briefly presented on June 30, 1960 at the First Congress of the IFAK, 
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ficiently rapid” is more precisely given for the class of problems considered in [4]). In the t, x phase space of model 
(2), these regimes correspond either to those portions of the surfaces (3) which are threaded through by trajectories of 


the corresponding ordinary system of differential equations 


X= Flt. x, @(s)], G= a(t, Gn (k= (4) 


in one direction, or those portions from which the trajectories of the local system corresponding to (4) diverge toward 
different sides, each into “its own", from the subspaces 


(t,x) hy, hep, ...). (5) 


But to neglect the continuous change of the output coordinate g(t) is inadmissible in describing those motions 
of a real system (sliding regimes, motions with delays) which are connected with the coordinate o(t) lingering near 
the switching thresholds, The discontinuous idealization of the function ¢(o)leads to the appearance of the surfaces 
(3) of portions threaded through by trajectories of system (4) passing from one subspace (5) into another, changing di- 
rections in so doing. The construction of the correct model within the framework of the chosen structural scheme now 
requires a more accurate expression of the dependence of ¢(o) near the surfaces (3), 


Before proceeding to this more accurate expression, it is of interest to test, directly from the discontinuous ide- 
alization, the following simple approximation to the relay characteristic at switching points: 


(Sx) = Ex (k= ha, (6) 
where £) may be an arbitrary number in the interval 


lim @ (co), lim @ 
Here and below, lim and lim are the least and greatest limit points respectively [5]. Such an approximation, already 
used successfully by Coulomb in studying dry friction, retains all the advantages of the discontinuous idealization, yet 
being in a certain sense continuous: the equations ¢ = yo) determines on the o, ¢ plane a line, any portion of 
which in general may be sketched without lifting the pencil from the paper (see, for example, Fig. 2a). 


The goal of this paper is to show how approximations of the type (6), (7) may be used in constructing a “con- 
tinuous” variant of the theory of relay control systems (2), | 


It is easy to see that approximations (6), (7) are strongly connected with a representation of the relay character 
istic go), as well as of discontinuous functions, as limits of nonuniformly convergent sequences { y(o)} of continu- 
ous functions, monotone and changing sharply near (3). The validity of such representations is discussed, and the rep- 
resentations themselves used, in many papers having to do with the theory of relay differential equations, Actually, 
the characteristic go) under approximations (6), (7) may be considered as a point set in the o, ¢ plane, the limit 
of the same sequence {¢j{0)} when each term of this sequence is also considered a set of points(o, ¢); and as a 
measure of the nearness of sets, their mutual deviation according to Hausdorff [6] is taken in the closed region of the 
0.9 plane with which we are dealing. We note that the comparison of characteristics ¢){ 0) with the aid of their 
mutual deviation in the o , ¢ plane certainly corresponds better to one's intuitive idea of the nearness of curves than 
does the comparison of the functions (0) for the same values of their argument ©, as is done in the theory of series 
of functions. 


2. Employing the usual geometrical interpretation [7], we shall say some more about the direction field of system 
(2); the directions of the field at a given point (t,x) shall be thought of as corresponding to the vector f(t, x), and con- 
sidered as a point of some n-dimensional space E(f)—the space of directions of the system (2). 


Taking the approximation (6), (7) for the relay characteristic g(o), it is natural to consider that on the switch- 
ing surfaces (3) the right sides of systern (2) give not one, but a whole cone of directions (Fig. 1); namely: the set 
{fft. x, of vectors f(t.x, corresponding to all numbers, from the interval (7). We note that at 
points (t*, x®) of the surface (3) the inequalities 


lim f(t, x) << f(t, x) S, jim 


(t, x*) 


m /f,(t,x) (i=1, 2,..., a) (8) 
t®,x°) 


x)>( 


for components of the vector function f(t,x) might not be fulfilled, An integral curve of the relay system (2) may 
now be defined as any continuous curve x(t) about which it may be said that motion along it proceeds at each mo- 
ment of time t in some direction of the cone {f(t, x (t)]}. Such a definition would seem in the present case to be a 
simple and accurate analog of the definition of integral curve of ordinary differential equations, 


The system obtained does not possess complete information for the determination of the velocity x(t) as func- 
tion of the coordinates t,x. The possible violation of uniqueness of the solution will be payment for the comparative 
simplicity of the model, reflecting the insufficiency of the information contained in the latter, The criterion of how 
successful such a simplified model will be is in the final analysis found by comparison with experiments, 


The assumed conception of a system with incomplete information naturally also includes those cases when the 
direction field determined by the functions f(t,x) is "washed away” not only on the surfaces defined, but also on 


at certain sets in phase space with nonzero volume, This latter circumstance is very essential, since for no single link 
idl in a real system may the output be characterized exactly for any lengthy period of time in the form of a single- 
di- valued functional or differential dependence of the input and time, For instance, the results of repeatedly photo- 

. oail gtaphing static characteristics of an electromagnetic relay often have an appearance like Fig, 2b rather than Fig, 2a, 
Relay switching proceeds in different experiments under differing values of the input o (the line density in the rect- 
angle 1234 is proportional to the relative frequency of relay switching at the corresponding value of o in the interval 

s ide- 12), In the first approximation, one considers a strictly determined characteristic (Fig. 2a), which is some kind of 
average or another of the real characteristic( Fig, 2b), A different approximation involves the consideration of char- 
acteristics in a form represented in Fig 2c, Here the coordinates o andg may define any point in the blacked out 
rectangle 1234, covering the hysteresis zone of a real dependence g(c), Further improvements in the model may 
consist in accounting for the statistical inhomogeneity of the region 1234, 

The above sketched principle for constructing mathematical models of relay systems reveal a definite similarity 
to the theory of equations in contingents and paratingents [1, 2, 8]. We shall give a short exposition of the fundamen- 
tals of this theory,® The exposition is given in terms as near as possible to the formalism occuring in the theory of 
ordinary differential equations. 

already 3. All spaces considered below are Euclidean, We recall that the distance p(p, q) between points p =(p,, pp..... 

on, yet Pm) and q =(q, @,...+4,,) of m-dimensional Euclidean space is given by the formula 

m 0.5 

(p, q) = [> (pi — (9) 

con- i=1 
and the distance p(A, B) between point set A and B of such a space is taken as the exact lower bound [5] of the set 

nee of quantities p(p, q where p and q are any points in the sets A and B respectively, 

ntinu- An abstract function which assigns to a point p = (p,, pz....» Pm) a nonempty bounded closed connected set 

1e rep- {f(p)} in some n-dimensional space E(f) is called a many-valued n-dimensional vector function f(p). In the scalar 

ally, case(n=1) the set {f@)} isa segment, for which we introduce the spe- 

limit cial notation [mf(p), Mf(p)}. It is further considered that the analytical 

as a expressions on the right side of systems (2), (6) and (7) provide some many- 

of the We of: fo valued vector functions f(t, x) by the superposition of the ordinary func- 
eir tions f (t,x, ¢) and the scalar many-valued functions ¢ = ¢[ o(t, x)], which 
es than S f= ong, were given in sec, 2, In accordance with this natural definition, and with 

f series account taken of one additional convention formulated in sec.4, a com- 

aga Ly ee pletely determined sense is given also to such expressions as the sum ¢,(0) 
system + of scalar many-valued functions ¢;( 0) and (superposition of 
ad call Fig, 1 the single -valued function u = v + w and the many-~valued functions v = 
and w=¢,(9)], the scalar product of two many-valued vector-func- 
tions of the same dimensionality and so on, The relations 
switch- 
set 


mV (p)>0 and MV (p) <0 (10) 
for scalar many-valued functions V(p) are customarily written in the form of the corresponding inequalities 


*A more detailed theory of equations in contingents is given in [9]. 
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V (p) >0 and V (p=) < 0 (11) 
(strict inequality in (10) will correspond to strict inequality in (11) and conversely), 


If the right side of equation (1) determines in some fashion an n-dimensional many-valued function of the 
variables t,x, then this equation will be considered as an equation in contingents, The space E(f) for the function 
f(t, x) we have already essentially agreed to call the direction space of system (1), and the sets {f(t,x)}, the 


direction cones, 
as any continuous vector function x(t)=x(tg, Xo, t) for which 


? 
a 
x(t») = X_ and whose contingents, considered as a point set in 


the direction space E(€) are included in the direction cone { f[t, x(t)]} everywhere in the interval (12): 


As is known, the contingent kont x(t) of the curve 
x = x(t) at the point (t, x(t,)) of the space t,x is defined as 
the set of all limiting directions of the chord spanning the 


? 
points (ty, x(%)Jand (t,,x(t,)) of this curve when > t,. 
g 1 2i¢ Definition 1, Following [1, 2, 8], we define a solution 
J 4 to the equation in contingents (1) passing through the point 
c 


(to, X9) and defined on some interval 


Fig. 2 t—act<ct+a (a>dd,) (12) 


kont x (¢) {f x (d)])}, (13) 


It is easily seen that definition 1 in its application to the relay equation (2), (6) and (7) corresponds very well 
with the considerations presented in sec, 1 and sec, 2, 


The ¢€-neighborhood S«, A) of the point set A is taken to mean the set of all points\p of the given space for 
which p(p,A)< 


Definition 2, The many-valued function f(p) is called 8 -continuous at the point py if for any € > 0 there 
existsa 5 = > 0 such that the set {f(p)} lies in the -neighborhood of the set {f(pp)} , provided the point 
p lies in the 6 -neighborhood of the point pp. -continuity shall require, in addition to the above, that the set 
{f (p)] in turn lie in the e-neighborhood of the set {f(p)}. 


This terminology is connected with the concepts of semideviation [6] 8 and deviation & of sets and we borrow 
them from the work of Yu. K, Solntsev [10], The term “semicontinuity from above with respect to deviation" [8] 
appears to be less convenient as it is more cumbersome, 


If at each point of the region G = G(t, x) of t, x space the many-valued function f(t, x) is 8 -continuous in the 
variables t, x and the set {f(t,x)} is convex, then we shall say that the contingent equation x = f(t, x) satisfies 
condition A in the region G, 


Following papers [1, 2 and 8] one may formulate this existence theorem. 


Theorem 1, If the contingent equation (1) satisfies condition A in a closed bounded region G = G(t, x), then 
‘at least one solution of this equation, defined on some interval of the form (12), will pass through each internal point 
of the region G, 


This theorem is very easily proved by the method of Peano, if, carrying through the discussion according to the 
scheme in [7], in the concluding portion of the proof one uses, in place of the mean value theorem of integral calcu- 
lus, the fact that the direction of the chord subtending the Euler broken arc belongs to the least convex cone which 
contains the directions of all links of this arc, It is apparently not possible to weaken significantly the requirement 
on the convexity of the set {f(t, x)}. 


Further construction of the theory of contingent equations is conveniently carried out under the following 
assumptions(integrals are to be understood in the sense of Lebesgue), 
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Theorem 2. Suppose the contingent equation (1) is given and satisfies the condition A in the closed and bound- 


ed region G = Git; x). The vector function x(t) is a solution of this equation on the interval (t,, t,) if and only ifthe 
following conditions are fulfilled for every t{, t® in the interval (t, tg): 


(t’, x (t’)) €G, =x (t") + \ @ (ED aE, (14) 
t 


where g(t) is some summable function for which 


p(t) for ti<t<te. (15) 

4. Starting from Theorems 1 and 2, it is not difficult to obtain for contingent equations a theorem on the continu- 
ability of solutions, analogous to the corresponding proposition in the theory of ordinary differential equations, and,in 
the case of continuability of solutions, to prove that the axiom of E. A. Barbashin on the generalized dynamic system 


[6] is fulfilled. The usual theorems on substitution of variables (analogous, for example, to Theorems 1 and 2 of [11)), 
are easily extended to contingent equations, 


Thus the consideration of system (2) as a contingent equation, based on approximations (6) and (7) allows one 
to obtain a correct model for a relay system without having recourse each time to a detailed construction of its phase 
space. This latter circumstance is fairly significant, since the construction indicated is cumbersome with two or 
more intersecting switching surfaces and is not too useful for “strongly washed away”® characteristics 9(¢), similar 
to that depicted in Fig. 2c. In fact, however, such an approach is admissible only under the condition that for the sys- 
tem examined the requirements of the existence theorem are satisfied. We briefly consider how natural these require- 


ments are in the case, encountered most often in practice, when the relay system may be described by equations of 
the form (2). 


Under an approximation of the type (6), (7), the characteristic g(¢) of real links (of the same number and 
similar to Fig. 2c), and also the jump-forming external disturbances are 6-continuous, so that the right side of equa- 
tion (2), being a superposition of continuous and 8-continuous functions, is always B-continuous. 


The requirement of convexity of the set { f (t, x)} is far more limiting. Actually, from the convexity of the 
set {g(o )} and the continuity (single-valued of the function f(t, x, g) the convexity of the set { f(t),x, 9(¢ )]} does 
not in general follow. In this connection, we return to the question of what one should understand by the result (t,x) 
=f[t.x, g(t.x)] of superimposing a many-valued and a single function f(t, x, ¢) and g(@) = y(t, x). Here the fol- 
lowing two limiting points of view are possible. The first, actually occurring in [12],reduces to the definition of the 
cone { f [t, x, y(t, x)]}at points (t’, x”) of the surface (3) as a set of vectors defined by the ordinary vector function 


under the condition that the numbers OG = 1, 2,...,m) take on all values from the segment (7) independently of 
one another.* At any point(t, x)the components f(t, x, ¢) (i = 1, 2,..., m) of the right side of system (2) prove to 


be not connected among themselves by scalar many-valued functions, and the set {f(t, x, ¢)} , a rectangular paral- 
lelopiped 


mf; (i, x, Qh < Mh X, (i=1, 2,,.., #) 


in the space E( f), is convex for any system (2), Definition 1 becomes near in substance to the definition of E. E. 
Viktorovskii [12]. 


The opposite point of view, taken in one form or another in almost all papers on the theory of relay equations 
or exampje, [11]) and characterized above in sec. 2, reduces to this: that only those vectors (16) for which ¢{*) 
=tk =...=€, = €, are included in the cone { f(t, x, o(t, x)]}. The cone { f(t, x, ¢)} degenerates into the 
arc of some curved line in the space E(f) and therefore is convex only when this curve is a straight line, Thus this 
approach comprises only those systems (2) which are representable in the form 


*In (16) and below the sign * denotes the transposed matrix. 


the 
ction 
fined as 
ng the 
, solution 
e point 
or which 
3) 
there 
the point 
ie set | 
borrow 
on® [8] 
ous in the 
fies 
then 
point 
bral calcu- 
which 
irement 
| 717 


x=f, (t, x) +8, (t,x)@ (0) x), 


where f,(t, x), £,{t,x) and o(t, x) are ordinary continuous functions according to Definition1, and Theorem 1 becomes 
more effective than in the preceding case, when each system (2) was often considered as a system in n relay elements, 
In particular, the method of constructing Lyapunov functions, worked out for control systems with one nonlinear link, 
is not applicable to relay system (2) with the Viktorovskii definition of solution. 


The latter point of view corresponds, strictly speaking, only to the situation when the real link actually has one 
output (one pair of contacts on the relay, one point for removing the voltage from the trigger, etc.) whereas the first 
of the limiting approaches considered assumes the existence of n outputs from the link, independent at moments of 
switching. Such an assumption will apparently be natural in most cases; however, there may occur cases (when the 
relay has several groups of contacts, the output voltages of the trigger are taken at various points in its scheme, etc.), 
when it may prove useful to have some intermediate point of view, employing not one and not n, but two or three 
characteristics, identical everywhere with the exclusion of the sliding regimes, for the description of the real link. 
Then a compromise is attained between the necessity of describing the link by several characteristics, and the endeav- 
or to simplify the equations of motion. Allthe approaches given to determining the cones { f(t, x)} involve alike the 
theory of contingent equations, Results near in substance to the results [11, 12], are obtained from this theory, as was 
already noted above, in the form of limiting particular cases, We note that the theory developed in [11, 12], was 
presented in a form difficult for its application to the problem of automatic control, and the Solntsev definition [10] 
does not correspond to the physical considerations presented above. 


Uniqueness may be attained in determining the direction field for relay systems of the form 


x=f [t, (5), P2 (62), Om (Sm)], 
(t,x) =4,..., m), (17) 
under the condition that the second, more effective of the limiting approaches examined takes place for those and 
only those functions ¢(o) which are designated in (17) by undistinguished symbols. Thus, for example, the definition 
of solution according to Viktorovskii means that each system (2) is considered as a system 


x = (f, [t, x, 9: fn [t, x, (0)))’, 
and in the examples introduced below (see secs, 5, 6) the cones { f(t, x)} psenape cena 
E(f) (compare with [11], page 105). 
5. In the extension to contingent equations of Lyapunov's direct method, it is useful to use the following assertion, 


Theorem 3. Suppose the contingent equation (1) satisfies condition A in a closed bounded region G(t,x), and 
the continuous and continuously differentiable function V = V(t, x) is such that for the scalar many-valued function 


V 


the following inequality is satisfied everywhere in G(t, x): 
W (t, x) < 0. (19) 


Suppose, further, that some arc of some integral curve x(t) of system (1), corresponding to the segment 


(20) 


is completely situated in Gct, x). Then the function V(t) = V{t, x(t)] does not increase in (20), and we have the fol- 
lowing estimate: 


V (ts) —V (ts) < (te — 4) sup (MW Ir, x (21) 


The paoof of Theossm.$ is givenin Appendix 1, The scalar product { (¢, x) of the row matrix = 
3; the set ( f(t, x)} is a segment formed by numbers which are scalar products of n-dimensional vectors 


av / x and E(t, x) under the condition that €(t, x) is an arbitrary vector of the cone { f(t, x)}. The computation 


of the function (18) for system ie is accomplished completely analogously to the method in the case of ordinary 
differential equations, 


If stability is understood in the sense of the usual definition [13], then on the basis of Theorems 1-3, considering 
the many valued function (18) in place of the ordinary derivative of the Lyapunov V-function, it is easy to establish 
the validity of Lyapunov's theorems on the stability and instability of motion in the supplemented formulation intro- 
duced in [13] for contingent equations. The proof is a verbatim repetition of the proof of the classical theorem. 


In the way of an example, we investigate, by the method of Lur'e, the stability in the large of the equilibrium 
of the system of contingent equations 


x =Ax-+aq(s), c= k’x, (22) 


describing normal modes of motion of a control system with one nonlinear link (A is a square, and x, a, k column n- 
dimensional matrices), We assume that the linear part x = Ax is asymptotically stable, and that the function o(¢) 
is 8-continuous and satisfies the condition g(¢)0 = 0, i.e. has, for example, the form represented in Fig. 2c (the 
blacked out regions cover the hysteresis zone or the possible rapid fluctuation and drift characteristics of the nonlinear 
link). Taking the Y-function as a quadratic form V(x) = x'Hx(H =H"), we find that the many-valued function (18) is 
determined in the given case by the equations 


W (x) = — x’Gx + 2a’Hxg (oc), G=—HA—A’H, H=L(@). 
Suppose the relation 


2L (G) a = — k, 


holds, where G is some symmetric matrix for which 


x’Gx>0 for x0. (23) 


W (x) = — x’Gx — og (0), (24) 


Wi(x)<0 for x+0. (25) 


Since the form x'L(G)x is positive definite (Lyapunov's theorem), the stability of the equilibrium x= 0 in the 
large will follow, in accordance with the above, from Lyapunov's theorem [13]. Inequality (23) is satisfied if 
x’Gx = x’Gix + (u’x)®> (G = Gi + w’), (26) 


where 2G,, is some positive definite form, and y a real column matrix. Representation (26) is possible under the 
condition that the equation for wu 


2L (uu’)a + k + 22 (Gi) a=0 (27) 


admits a real solution, This is known to take place for any matrix G, with sufficiently small elements, as well as 
for those which correspond to the positive definite form x"G,x, if the shortened resolving equation 


2L (uw’ ja +k = 0 (28) 
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possesses simple roots. The method of investigating equation (28) is indicated in [14]. We have generalized the 
results obtained in [15] for systems of the form (22) with g(o) = g@,t): in the example analyzed, the nonlinear 
characteristic may depend in a different manner on time along different trajectories of the system (this occurs, for 


example, in the reversal of magnetization in magnetic relays). 


In order to extend the contingent equations, the theorems of Barbashin and Krasovskii ({13], Theorem 14.1), 
Krasovskii ({13], Theorem 15.1), Wechsler [16], and other applications of the direct method of Lyapunov, which 
permit one to weaken the requirements on the V-function and their derivatives, it is necessary to use also the fact 
that contingent equations under certain natural restrictions (see sec, 4) determine a generalized dynamic system. An 
illustration of this is the proof of the following theorem given in Appendix II. 


Theorem 4, Suppose the contingent equation 
x=f(x), (0) (29) 


satisfies condition A over the entire t, x space, and all its solutions x(t) are continuous for -« < t< +. If there 
exists a continuous positive definite function V = V(x), infinitely large for x - @, not increasing along any solution 
x(t) and such that no one of the surfaces V(x) = c for c> 0 contains any solution x(t) completely, then the equilibrium 
‘x= 0 of the system (29) is asymptotically stable in the large. 


6. As is shown by the investigation of system (22) carried through above, significant information on the solutions 
of equation (2) may be obtained comparatively simply if the function (18) is considered in place of the V in the 
theorems of the direct method of Lyapunov. However, the limitations imposed here on the V-function are not known 
to be excessively rigid except under the condition that motion resulting from the system studied actually may pro- 
ceed at each point (t, x) in any direction belonging to the cone { f(t, x)}. If the right side is ambiguous only on 
certain surfaces of the t, x space, then this latter condition is violated. In accordance with Definition 1, the motion 
often proceeds in one definite direction, and to make the stability condition more precise, the measurements of the 
V-function should be restricted to this direction only. 


A possible procedure for making the investigation of stability in the given situation more precise will be con- 
sidered in connection with the contingent equations 


x = F(x) + a(x), F(O)=0, 0€ (0), (30) 


given over the entire phase space x Here P(x) is a continuous single-valued vector function, a is a constant column 
‘Matrix, g(x) is a 6 -continuous scalar function, ambiguous only on the sufficiently smooth surfaces Sy, determined by 


the equations 


(x)= 0 (k=...,—1,0,4,...) 


and separating the x space into nonintersecting regions. We denote the regions separated by the surface S, as Gy, and 
Gx, supposing that ok(z)>0 in G,, and o}(x) <0 in Gy. 


The equations of a controlled system with one link with ambiguous characteristic may be brought into the form 


(80). 


By definition 1, motion in the regions Gt proceeds in accordance with the equations 
x = F (x) + agp (x) = ff (x), 
(x) B GY, 
9, (x) BG, (31) 
(the superscript signs are chosen to be the same in Gy, and (x) ), and within the set 


=| 


(x)= 0, Of (x) <0, of (x) (x), 


0 
(x) = UF (x) + (x)] 


(32) 


in accordance with the system 


| 
720 


X= fh (x) = F(x) + af; (x), 


lor 
03, 
)s 3} (x) = aT (x) =0. (34) 
= From (34) it follows that 
06 
(x) = — (Fea) (85) 
me (x) < (x) < Mg (x). 
ti 
anil It is easy to see that from (32) the relation 
tions 
follows. 
known 
pro- If, on the surface Sy 
on 
jotion 
f the mo (x) = min (Gf (x), (x) = max (x), (38) 
which is always fulfilled for the characteristic of a real link, then (36) is valid only in the set (32). 
con- On those portions of the switching surfaces where 
ot (x) (x) <0, of 
(39) 
Definition 1 leads to the system 
+ 
G,, and x = F(x) + af, (x) = f, (x), 
he form 


Oy 0, ay = (x). (41) 
Since on the set (39) necessarily 0 and 
0s 
& axo, 
then on this set 0, (x) = 6, (x) = 0 and mp (x) < & (x) < Me (x), 


(42) 


= — 2 F(x) (43) 


(33) 
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When inequalities (42) are satisfied over all space, the motion on those portions of the surfaces S, which enter 
into none of the sets (32) or (39) proceeds under condition (38) along the semitrajectories of the local system (31). 


The structure of the phase space of system (30) on the boundaries of the sets (32), (39) is not explained in de- 
tail here, since the existence and continuability of the solution of this system is always guaranteed in practice by the 
general theorems of contingent equations. Definition 1 leads in the given case basically to those representations of 
the trajectories which are usually used in writings on the theory of relay systems. The difference consists in the fol- 
lowing: it is usually considered that, moving along the system's trajectories, it is impossible to fall into the zone of 


dynamic indeterminacy 


(x) >0> (x), (44) 


from which the trajectories of system (31) diverge to different sides. According to Definition 1, there is a trajectory 

) of system (33) lying on the surface S;, passing through each point of zone (44) (see Fig. 3a). Such an approach is 

better founded since the phase picture of a relay system near the zone (44) should be considered a limiting case of 

the phase picture of an ordinary dynamical system with sharp dependence of the solution on initial conditions (Fig. 
3b). This difference presents no obstacle to carrying over direct- 
ly to contingent equations many results (in particular, those 


RO DP et, obtained in [3]) of the theory of differential equations with dis- 
Sy continuous right sides, 
a 


ieee os Suppose the equation of the switching surface S, for system 
(30) is represented in the form 


Fig. 3 


and there is constructed a continuous positive definite function V(x), growing to infinity asx -+@, with continuous 
partial derivatives V/ dx in the regions Che and the partial derivatives dx of the function Vi(x;, Xg,..., 
Xp -;),obtained from V(x) by the substitution (45),are also continuous. We introduce the notation: 


where ff, fj, 1” are the right sides of systems (31), (33) and (40) respectively. If the inequality V(x) <0 is valid 
over the entire space x for the function 


(*) in the region GE. 
v (x) = \¢ (x) in the region Si. (47) 
lo (x) in the region Sy 


where Sk and Sf? are given by the relations 


(x) ~ Sta =0, 2 (Ft @) a+ 0, (48) 


and V(x) = 0 only on some set M not containing whole trajectories of system (30), different from the equilibrium 
x = 0, then one may consider that Theorem 4 is applicable. We use this fact for transferring Less algorithm [14] 
to the particular case of system (30): 


x=Ax+ag(), = b’x— ppd), 


p>0, T?=p+b’/Aa>d0, (50) 
assuming that the eigenvalues of matrix A have negative real parts, and the function 9(o)is such that 9(¢)o > 0 for 
o #0 and 51 

= 00 (51) 
0 


| 6, (x) = 0, ra +0, 


valid 


3) 
rium 


um [14] 


0) 
yo > 0 for 


1) 


We shall prove that if the equation in the column matrix u 


(p>0), (52) 
where the matrix L(G) is determined for a symmetric matrix G by the relation 
G = — L(G) A— A'L (@), 


has a simple real solution, then the equilibrium x = 0, o = 0 of the sytem (50) considered is asymptotically stable in 
the large. 
If equation (52) is soluble in the indicated sense, then the equation 


2L (uu’)a +2 Vp u+b+2L(G,)a=0 (53) 


is also soluble, where G, = Gj is any square matrix with sufficiently small elements. We choose G, such that the 
form x"Gx is positive definite. Then for the positive definite Lyapunov function 


V(x,0) = x'L(@)x+ (e@d, +e 
0 


the function (46) assumes the form 


V, (x, 0) = V, = Vy (x, 0) = — x’Gix— (u’x + (54) 
(x, 0) = — x'G@x — 2a’L (@) (x, 9), (55) 
(x, 0) = — x’Gx + 2a’L (G) xf, (x, 0), 
(56) 
and the regions Sf, (48) and Sj” (49) are determined respectively by the relations 
C= Ox, p> 0, (57) 
o=60,=p=0, b’a +0. (58) 


For p > 0 the set (58)is empty and in set (57) o), = b’x — pf, (x, 0) = 0, so that the function (55) is 
representable in the form 


8°) = — x'@x + 20’L (G) + — = — x’Gix — + VpE,)* 


and proves to be negative definite. According to (51) the function (55) vanishes only at the pointz =0,¢0 = 0. The 
conditions of applicability of Theorem 4 are satisfied (the set M coincides with the origin x = 0, o = 0). 


According to (53) in the case p = 0 in the set (58), 


(x, 0) — x’Gx, 
and the set (57) is empty. By the substitution x(t) = 0(t,< t < t,) we convince ourselves that the straight line x = 0, 
on which only the equations Y= 0, V = const are possible, does not contain a single arc x = x(t), @ = o(t)(t, < t< 
t,) of any trajectory of system (50) different from the equilibrium x = 0, @ = 0. Again the conditions of applicabili- 
ty of Theorem 4 hold. 


The author thanks E, A. Barbashin for guidance in completing this project. 
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APPENDIX I. 


We prove Theorem 3. Consider an arc of the solution x(t), spoken of in the conditions of the theorem, Since 
the function x(t) is absolutely continuous on the interval (20) (Theorem 2), and the function V(t, x) by virtue of con- 
tinuous differenti ability satisfies in G(t, x) a Lipschitz condition in t and-x, then the compound function V(t) = V{t, 


x(t)] is absolutely continuous on (20), since for any t'é [ty, %] 


t’ 
V(t’) = V (ts) + ( a. (59) 


On the segment (20) the vector function (=). - and the partial derivative (5), ” are continu- 


ous; therefore according to the chain rule of differentiation, the derivative W(t) = 3- V(t, x(t)] coincides with the 
function below (60) at points where the derivative x(t) exists, i.e. almost everywhere on the segment (20): 


Moreover, at these points x(t) € { f[t, x (t)]} (Definition 1), so that according to the definition of algebraic 
operations on many-valued functions (sec. 3), almost everywhere in (20) the relation 


V “av 


is satisfied, as well as the inequality 


mW (t, x (0) <V (0) < MW (t, x (0). (62) 


The assertion of Theorem 3, in particular the estimate (21), follows from (62) and (59) by the mean value 
theorem of integral calculus. 


APPENDIX II. 


We prove Theorem 4, From the conditions of the theorem, it follows that for any solution x(t) = x(x», t) [x(%, 
0) = x9) of equation (29), the limit 


[x (xo, = (xe) 
exists and the semitrajectory x(%, t), O<t< + © is situated in the region 


V (x) < V (xo), (64) 
the region (64) being bounded and shrinking to the origin as x 0. Thus stability of the equilibrium x =0 in the 


sense of Lyapunov is present, and the set Q(X») of w-limit points [6] of the trajectory x(x, t) is nonem>ty. Accord- 
ing to (63) and the continuity of the function V(x) 


V @) = if (xo). (65) 


By following papers [1, 2], it is not difficult to show that under the conditions indicated in Theorem 4, the 
mapping accomplished in the space x by trajectories of equation (29) satisfies the axiom of a generalized dynamic 

system [6], so that the theory developed in [6] is completely applicable to this equation. In particular, according to 
Theorem 3 of [6], for each point q €Q (x9) there exists a trajectory 9(q, t) satisfying the condition 


@q,9EQ (xo) + cc). (66) 


| 
| (61) 
| 
\ 
124 


ntinu- 


‘ic 


from which in correspondence with the conditions of this theorem, it follows that V(x») = 0. The validity of the 


equation jim |X (Xo, t) = 0 now becomes clear, if one bears in mind the positive definiteness and continuity of the func- 
tion V(x). 


From (65) and (66) for this trajectory, there follows the equation 


V = (x) = const 00 Kt< + o), 
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CONCERNING THE PROBLEM OF DETERMINING 
THE STABILITY-SECTOR BOUNDARIES 
IN NONLINEAR CONTROL SYSTEMS® 


A. K. Bedel'baev 


(Alma-Ata) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 7, 
j pp. 831-833, July, 1961 

Original article submitted November 9, 1960 


A method is presented for determining the boundaries of the sector in which a system of automatic 
control is stable “in the large”. 


Let a system of indirect control be described by the equations [1] 


>) +B 1,2,..., 0), 


k=1 


where bsk, Ng, j, and r are constants,and {(@) is a nonlinear function the accuracy of which is restricted by the 
relationships 


(1) 


7 (0) = 0, mo® < of (c) < Mo’. (2) 


In the inequalities (2), m and M(0< m <M< @are certain constants which are to be determined. In accord- 
ance with (2), we shall consider system (1) as being absolutely stable over the sector (m, M), provided it is asymptot- 
ically stable at any initial deviations and at any nonlinear function f(¢), as long as it satisfies the conditions (2). 

In the following we describe a method of determining the boundaries of sector (m, M) in which the system (1) is 
absolutely stable. 


Let the roots Ay, Ag,..., Ap, of equation 


(—1)" D (A) = Jou — uJ = 0 


have the negative real parts 


Rea, < 0 


(e = 1, 2,..., a) 


and be such that the condition [1] 


=r+.>) is Bo 


is satisfied. 


* We wish to state that this is not a new problem, and that it has been investigated earlier by a number of authors, 
including M. A. Aizerman, B. S, Razumikhin, V. M. Popov, V. A. Pliss, E. N. Rozenvasser, and others, 


| 
= 
(3) 
| (4) 
| 
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Then the solution of the equilibrium-position stability (ny = ng=...0, = & = 0) of system (1) can be obtained 
from an analysis of the equations 


«.= > + 1, 2,..., n), — Tp (0), 


=+ f(). 
Let us examine the positively defined Lyapunov function 


DD 


s=1 k=1 


where Gs, are determined from the formulas (2.54) of work [3]. The time derivative of this function, computed 
from (6), will be 


—v =>) + 20 >) pare + rp(o) 


s=1 k==1 k=1 
— 2p, = jx + Pep (0), = (ae + Ons) Ms, 


and ag, are arbitrary constants selected in conformity with Silvester's inequalities; 


(k = 1, 2,..., 


Consequently, the expression for —¥ will be a positively definite form of the variables x;,%, . . 
its discriminant is larger than zero: 


From here, it is readily found that [3] 


kK=1 


where A is a determinant composed of ag,, and Ag, are the determinant's minors with signs, Let 


where 
where 
e 
a3, a a 
11 41g - - - Gin Pi 
1 -- Ong Pe 
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c= Aisixs 


k=1 k=] 


Be = Dy (Ax + Abs) 


k=1 


On the basis of (7), the inequality (9) will then have the form 


2%z+c< 0 (2 =H(3)). (12) 


It is evident from (8) and (10) that a > 0 andc> 0. Hence the inequality (12) is valid only at the condition b< 0, 
i.¢., 


D 4rA. 


s=1 


This being established, we write inequality (12) in the form 


where m and M are defined by the formulas 


—b— V—ac b? — ac 5 
with the provision that the condition b* = ac, which follows from the fact that values of function #(¢) are real 
values, is satisfied. From (14) it is now easy to obtain the inequalities (2), at which, as we see, the system (1) is 
absolutely stable. 


In conclusion, we note that the expressions for m and M contain a,, values which satisfy only the condition (8). 
By varying these values, we can both extend and restrict the boundaries of the sector (m, M) in accordance with the 
conditions of practical problems. 
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THE PROBLEM OF SPEED OF RESPONSE WITHOUT SWITCHING 
FOR AN ARBITRARY NUMBER OF CONTROL FUNCTIONS 


A. N. Filatov 


(Tashkent) 
Translated from Avtomatika i Telemekhanika, Vol, 22, No. 7, 
pp. 834-837, July, 1961 

Original article received March 3, 1961 


The author demonstrates the existence of a certain class of linear systems for which the optimal 
processes with respect to speed of response are performed without switching of the control functions 
for transition of the image point from any specified point on the phase plane onto a specified 


manifold. 
The problem of speed of response in the theory of optimal processes resides, as we know, in the following [1]. 


Assume the equations of motion for the image point (x;, Xg,..., Xp) in the n-dimensional phase space X, have 
the following form: 


hi (x, Ze, ++ In, U (t)) (i == 1, . @). (1) 


Here u(t) is a piecewise-continuous control function, and | u(t)] <M. The functions fj(x,, Xg,..., Xn, u)are 
continuous over the set of variables x, Xg,..., Xp» U and are continuously differentiable with respect to x;, Xg,...+%p- 


It is necessary to choose the function u(t) in such a way that the image point in leaving a specified position 
(xf, reaches another a priori specified position (xy, x},..., x4) in the least possible time. 


If the right sides of equations (1) are linear functions with respect to the variables x;, X%g,..., X, With constant 
coefficients (the case of linear systems), then the control function u(t) must as a rule change its sign (switch) during 
the process of movement; this is, of course, provided that the initial conditions are not chosen in any special manner. 
For example, even for the very simple system 


it is impossible to proceed from the arbitrary point (x}. xg) on the phase plane to (for example) the origin in the 
shortest possible time without switching (with the exception of those cases where the initial point is located on the 
switching curve). However, if we pose the analogous problem (for linear systems) with the requirement that the 
image point of the phase space X, must make the transition from the specified position (x}, xg...» X,) onto a 
smooth manifold of dimensionality s<n specified in this space, then it is possible to indicate the conditions for which the 
control function will not switch during the entire period of movement for any initial data provided only that there 


exists an allowed control function which transfers the image point from this initial position onto the specified mani - 
fold. 


This paper is devoted to finding the above-mentioned conditions, The analysis is performed for two variables 
(we study a phase plane). A straight line is taken as the smooth manifold. 


1. Assume the system of equations according to which the image point on the phase plane x,0xg moves is 
written in the form 


Ly = + + + +... + b,,u,, 


= + + + + . . . + 


is 
(2) 
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where aj; and bj; are constant; uk(t) (k = 1,..., 1) are control functions; here | =< My forty <t=<t. We 
shall assume that the functions u,(t) are piecewise-continuous. The problem resides in the following: it is required 
to find the conditions for which the image point moving in accordance with equations (2) from any a priori specified 
position (xj, xg) will arrive on the specified straight line x; = 0 (i.e., on the Ox, axis) during the briefest possible 
time without switching of the control functions u, (t) (of course, for the condition that an allowed control function 
which transfers the image point from the specified state onto the specified manifold exists), We shall use t) to de- 
note the instant at which the image point leaves the specified state, and t, to denote the instant at which it arrives 
on the specified manifold. 


We shall solve the problem formulated above using the principle of the maximum and the transversality condi- 
tions [1]. We shall formulate the auxiliary system of equations 


= — ahi — = — — (3) 


Solving the system of equations (3), we find 


y, = — pat + 4a) ett, ap, (0) = 
12 


where A,y<A, are the roots of the characteristic equation 
A? + A + — = 
for the system (3). In accordance with the transversality conditions we have the equation 


+ = 0, 
at the instant t\, from which it follows that 


+ 


In accordance with the principle of the maximum, 


ux, (t) = My sign (diets + dex) 


We formulate the expression 


Introducing the substitutions 


bo by, 


we write (10) in the form 


~ 


| 
(7) 
Therefore 
af (8) 
| (= '.....7). (9) 
730 


eci fied 
ble 
tion 

to de- 
rrives 


condi- 


Thus, the change in the sign ofthe function u;,(t) depends only on the change in the sign of the expression 


Sox, (13) 
where 


— — > 0. 
Expression (13) does not change sign in the following cases: 


a) On. > 0, Sox < 0, b) Oi > 0, Ox < 0; 
c) On <0, Ox < 0; qd) On < 0, Ox > 0, 
e) Tx > Sx > 0, f) Ox < On <0. (14) 


We shall clarify the conditions for which inequalities (14) are valid. As an example, we shall study the 
inequality “a” in (14): 


box bu —— > 0, Dex <0. 


From this it follows that 


for 2,7 % 


6 
ok*12 1k%29 12 1 
<A 
4 


Therefore condition “a” in (14) is valid only when the inequalities 


for 


bie bon tia — 


are satisfied. 


In analogous fashion we find the inequalities on whose fulfillment the remaining conditions in (14) depend: 


b) as > 0, Ay < bx < dea; (16) 


c) 


d) 


e) 


—b 
—*>0, <n < (20) 


Note that the equality sign in conditions “e" and *f* in (14) is reached either for by, = or for Ay = Ag. The 
latter is impossible by assumption (the roots of the characteristic equation (5) are assumed to be real and different). 
For the case bk = 0 (k= 1,..., 1) it is evident that 64, = Gy, and expression (13) does not change sign. 
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Furthermore, it is evident that in finding the solution (4) we assumed a, #0. However, if ay, = 0, then the 
sysem (3) becomes 


Pr = — — = — 


from which it follows that 


We (t) = 


In view of the transversality condition, p,(t,) = 0, i. e. p,(t) = 0. 
Therefore, = ce and 


Uy (t) = My sign = My sign 


From (21) it is evident that the functions u_(t) will not change sign, 
All of the above can be formulated in the form of the following theorem. 


Theorem. If the image point on the phase plane moves in accordance with equations (2), and if the transition 
from the a priori specified position (xy, x$) on to the straight line x, = 0 is possible by means of an allowed control 
function, then for arrival on the specified straight line from the specified point by means of an optimal control func- 


tion without switching during the shortest possible time, it is necessary and sufficient that one of the following condi- 
tions be satisfied; 


1) 2) a< 


#13 1k 


< Ao, 
3) 4) be = 0, 


5) ai, = 0. 


Here k = 1,..., 1 throughout. 


Remark 1. The conditions of the theorem are also retained for the case where any straight line parallel to the 
Ox, axis is used as the specified straight line. 


Remark 2, If the straight line onto which the image point arrives is assumed to be a straight line of the form 
4% + CigXg + Oy = 0, then by performing a transformation of variables in (2) according to the formulas 


x, = X,cosp— X,sing, 7, = X,sing + X,cosq, 
where 


we reduce the problem to that just treated above. 
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ON A CERTAIN PROBLEM IN THE THEORY OF PROCESSES 
WITH AN OPTIMAL SPEED OF RESPONSE IN LINEAR SYSTEMS 
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1, We shall study the equation 


(1) 


which is defined in the n-dimensional space R", where b is a fixed vector in R" and A is a linear transformation of 
R® which does not depend on time, The class of allowed functions from which we choose the control function u(t) 
is the class of piecewise-continuous functions with | u(t)| <1. We assume that the nondegeneracy condition (cf. 


[2}) is satisfied, i.e., b does not belong to any invariant sub-space with the dimensionality = n-1 of the transforma- 
tion A. 


Assume FCR®" js a closed convex set. We shall state that the control law u(t) (0 = t = T) transfers the point 
% €R"-F to F if the corresponding solution x(t; Xp, u(t)) of Eq. (1) satisfies the conditions x (0; X, u(t)) = X», x(T; 
Xp, u(t)) EF. 


The control law u(t)(0 =t ST) which transfers the point x» into this set F shall be called optimal if for any 
allowed control law u® (t)(0=t =T®) which transfers x» into F the condition T* - T = 0 is satisfied. 


Such a statement of the problem of optimization with respect to speed of response (when the final point on 
the trajectory belongs to a specific closed set) has been studied with a varying degree of generality by the authors of 
(1, 4, 5 and 6] who dealt with the theorems for the existence and uniqueness of an optimal control law. However, in 
certain engineering respects we can only assume that the problem has been finally solved when it has been estab- 
lished that a moving point which has “passed due to inertia” through the final point on the optimal trajectory and is 
continuing its motion along a new optimal trajectory does not emerge from the vicinity of the set F, i.e., it must be 
established that the optimization is in a certain sense “stable”, In the case where the set F can be reduced to a 
point, such stability is assured by the conditions which are imposed for the possibility of synthesis (2,3), This paper 
is devoted to an investigation of this problem; the most complete results are obtained for the case where n = 2 and 
F is a straight line. 


2. If w(t) is an optimal control law which transfers xp into F, and x (T; X9 u(t)) = x, @F, then the corresponding 
optimal trajectory intersects F only at the point x, and u(t) is an optimal control law in the sense of [2, 3] that 
transfers x» to x,. From this it follows that the optimal control law under study satisfies the maximum principle 
derived by L. S, Pontryagin [3}. 


The following propositions [1, 4, 5 and 6] apply. 


A. If u(t)(0 st <T) is the optimal control law that transfers x, into F, then there exists a solution » (t) of 
the equation y= -A*p (where A® is a transformation which is conjugate with respect to A) such that u(t) = signp 
(t)-b, p (T)-(y-x,) = 0 for all y€F, where x, = x(T; Xp, u(t)); here m+n denotes the scalar product of the vectors 
m and n. 


B. The optimal control law which transfer xp into F is unique if it exists, 


C. If there exists an allowed control law u(t) which transfers x» into F, then there also exists such a control 
law which is optimal. 


5 


The propositions formulated above make it possible to synthesize an optimal control law for Eq. (1) if the set 
F is an (n-1)-dimensional hyperplane (or a convex smooth hypersurface). As we know [3], this synthesis problem con- 
sists of finding the function v(x) = u,(0) where u,(t) is the control law which transfers the point x into F; v(x) enters 
into the equation 


= Ar bo (x) 


for the optimal trajectories of Eq. (1). 
In proposition “A* the latter relationship becomes 


p(T) =a, (3) 


where a is a vector normal to the hyperplane F and directed in the opposite direction from x; therefore the vector 
&(T)isdefined in a single-valued fashion. This leads to the following method for plotting the optimal trajectories 
of Eq. (1). 


Assume that the hyperplane F passes through the point x = 0 and that the real parts of all the eigenvalues of 
the matrix A are negative. Then for any point in space there exists an “allowed” (and therefore an optimal) control 
law which transfers it into F, Assume also that the scalar products a-b # 0; moreover, it is easy to establish the fact 
that the analysis given below can also be adapted to the case a-b = 0. Assume x is some point on the hyperplane F 
that is the end of a certain optimal trajectory. If the corresponding control law is such that u(T) = sign a-b, it fol- 
lows from proposition °A™ and relationship (3) that a*(Ax+b sign a-b) = aA-x + | a-b]2=0. Analogously, if u(T) = 
-sign a+b, then a+(Ax - bsign a-b) = aA-x -|a-b |= 0. Therefore the set S of points on the hyperplane F that are 
simultaneously the ends of two optimal trajectories (located in the respective half-spaces into which F subdivides 
the space R"), is determined by the inequalities 


aA-z-+|a-b|>0, aA-x—|a-b| <0, 


with the possible exception of the points at which these inequalities become equalities. The set S is never empty 
since it contains the origin. It coincides with the entire hyperplane F if the vectors a and aA are linearly dependent; 
it reduces to a single point if a-b = 0°. 


Assume F, is one of the half-spaces into which F subdivides R", and assume a is the vector representing the 
external normal to F with respect to F,. We shall assume that u(t) = sign (t)-b(t <0), where p(t) is determined 
by expression (3). In order to be specific we shall assume that a-b = +1. 


We shall always move in the direction of negative t. Assume % < 0 is the initial instant of the step discon- 
tinuity of the control u(t). We shall use P, to denote the set of points p€F from which the semi-trajectories x(t; p, 
u (t) ) (t =0) penetrate into the half-space F,. We shall continue these semi-trajectories to the instant %, or until 
they first return to F if this occurs earlier. It is clear that each of them is optimal for any point on it. 


Assume ]j is the set of points x(%; p, u(t)), p€P, belonging to R. We find that this set is located on an (n-1)- 
dimensional hyperplane. If I, is empty, then in view of the propositions cited above and our proposition, the semi- 
trajectories which have been plotted fill the entire half-space R. If I is not empty, then our semi-trajectories form 
a set R] which is bound by P;, I, and thesemi-trajectories originating from the boundary (on F) points P,. Repeating 
the reasoning cited above, we construct the set Tj by originating from Rg ; this set is formed by the semi -trajectories 
x(t-T; q, u(t) ), tT, q@T, which are continued to the next point 73,the step of discontinuity for the control u(t), oF 
until the instant at which they leave F, -Rj, etc, Such a geometric construction leads to subdividing the half-space 
F, ay mutually non-intersecting sets on each of which the value of the function v(x) is constant (v(x) = £1 for 


Repeating this procedure exactly for the vector -a and the corresponding set P,,we obtain the analogous sub- 
division of the second half-space k. 


*From the nondegeneracy conditions cited in section 1, it follows that these two conditions ae incompatible. 


2) 
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Examples, We shall examine the following system® in the plane R*: 


—ytu, 


where we shall take the axis 0x as F, and the unit length of the Oy axis as the vector a. The optimal control law for 
the points in the lower half-plane is written as 


_u(t)=sign(—sint—cost) (t<0). 


The graphical construction described above leads to the pattern of optimal trajectories shown in Fig. 1. These 
trajectories consist of circular arcs with centers at (1,1) and (-1, -1); the step discontinuities occur on the straight 
line x+y = £#(¥2 -2). The solid lines depict the trajectories corresponding to u = 1, and the dotted lines represent 
trajectories corresponding to u = -1. Solid lines are used to draw the trajectories which pass through the points (1, 0) 


and (-1, 0) and are part of the boundaries for the regions Rj and R,. The set S coincides with the segments [-1, 1] 
of the 0x axis. 


If we choose the axis Oy of the straight line x-y = 0 as the straight line F, then the pattern for the optimal 
trajectories will be as shown in Figs, 2 and 3, respectively. 


3. We shall now proceed to a study of the “stability” of the synthesis achieved in this manner, As we have 
already said, for stability of this kind, the trajectories which originate from the vicinity of the end of an arbitrary 
optimal trajectory must remain in the vicinity of the set F. The examples cited above demonstrate that this stability 


is not always present: it is violated at external points of the set S. Thus it is natural to introduce the following 
definition, 


* Although here the matrix A does not satisfy the condition of negative real parts for the eigenvalues, this can easily 
be achieved by an obvious modification of the right sides of the system without altering the qualitative pattern of 
the optimal trajectories; this pattern is especially easy to depict for the system under study. 
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If v(x) (v(x) = 2 isa function which is defined on R"-C, where C is a closed convex set inF, and satisfies 
the following conditions: 


a) each trajectory of the equation 


dz 
= At + bv (x) 


reaches C; the first point of its intersection with C is called the end of that trajectory; 


b) for each xg€R"-C the function uyo(t) = v(x (t; X9)), where k(t; x9) is a solution of Eq. (2) with the initial 
conditions x(0; x9) = X9, is studied only up to its end p(x,) and is an allowed control for Eq. (1). 


~ 


| Such a function v(x) shall be called a synthesizing function relative to F for Eq. (1). 3 
| This synthesizing function can always be deteremined by means of the optimal control laws u, (t) which 


transfer the point x€R"-C into C, if we assume V(x) = u,(0). For example, the optimal control laws which transfer 
the points in the space into the set F determine tne synthesizing function v}(x). 


Definition. The synthesizing function v(x) is called stable if for any x€R"-C and any €> 0 there exists an? 
> 0 such that all trajectories of Eq. (2) which originate from the n -vicinity of the point p(x) remain in the €-vicin- 
ity of the set F, 

As our examples show, the synthesizing function vp(x) is, generally speaking, not stable. But if F is a hyper- 


plane passing through the origin, and if the real parts of all the eigenvalues of the matrix A are negative, then it is 
easy to see that the synthesizing function vox) defined by optimal control laws that transfer the points in the space 


‘ 
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function exists, 


y 


to the origin is stable [2, 3]. Thus the problem arises of finding the “optimal” stable synthesizing function if such a 


We shall limit our analysis to the case of n = 2, when F is a straight line passing through the origin. Certain 
of the arguments cited below can obviously be generalized to include the case of arbitrary g when F is an (n-1)-di- 
mensional hyperplane. In addition to the propositions cited above, we shall also postulate that the vectors a and aA 
are linearly independent; it is easy to deduce the fact that in the contrary case the synthesizing function vp(x) isstable. 


line F; cf. Fig. 2) for both increasing and decreasing values of t. 


3. aA-(Ax,+ b sign a*b)< 0, i.e., the trajectory passing through the point x, penetrates into the second half- 


plane F, (cf. Fig. 1)* for both directions of variation for t. 


* The case aA - (Ax; + bsign a - b) = 0 cannot arise, since the vectors a and aA are linearly independent; this is also 


obvious from purely geometric considerations. 


- 


- 


~ 


If a+ x = 0 is the equation for the straight line F and x is a point on it at which aA-x+]a-b] and aA-x-|a-b| 
have identical signs (i.e., a point at which the straight line F intersects the trajectories of the equations x = Ax+b 
and x = Ax -b in the same direction) then this straight line intersects all the allowed trajectories passing through the 
point x in the same direction. From this it follows that the point x cannot be the end of a trajectory of Eq. (2) with 
the stable synthesizing function v(x). In fact, it is easy to see than when € is sufficiently small, any vicinity of 
such a point contains points with a property such that allowed trajectories originating from them emerge beyond the 
€ -vicinity of the straight line F before they reach F, Therefore, the ends of the trajectories in Eq. (2) which are 
determined by a stable synthesizing function definitely lie on the segment S of the straight line F with the ends x, 
and x» which are respectively determined by the equations aA-x,+| a-b| = 0 and aA+xg-Ja-b] =0. Since for 
a+b = 0 the set S evidently consists solely of the origin, it follows that in this case the function v»(x) will be an 
*optimal® stable synthesizing function, However, if a-b # 0, then the following three cases are possible. 


1. Ax,+b sign a-b = 0, i.e., the point x, is singular for the equation x = Ax+b sign ab (cf. Fig. 3). 


2. aA+(Ax;,+ b sign a-b)> 0, ie., the trajectory of the equation x = Ax+b sign ab in passing through the 
point x, passes into the half-plane F, (rel ative to which a is a vector internal normal with respect to the straight 
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Since x, = -x;, it follows that each of these cases involves the satisfaction of analogous relationships for the 
trajectories of the equations x = Ax -b sign ab at the point xp. 


Simple reasoning demonstrates that in cases 1 and 2, the synthesizing function V¢(x), which is determined by 
means of the optimal control laws that transfer the points of the plane onto the segment S, is stable. Under these 
conditions it is *optimal® in the sense that no trajectory of Eq. (2) with a different synthesizing function that is 
stable relative to F is traversed during a shorter time. In fact, since the end of such a trajectory must lie on §, it 
cannot be traversed more rapidly than the optimal trajectory which transfers the initial point under study to S, 


We shall now demonstrate that incase 3 the “optimal® stable synthesizing function relative to F is the function 
Vo(x). For this purpose it is sufficient to show that if v(x) is a stable synthesizing function, then all trajectories of Eq, 
(2) end at the origin. We shall assume the contrary. Assume C is a set of terminations of the trajectories of the cor- 
responding Eq. (2) on F. The set CCS is closed. Assume y, = 0 is the extreme point in this set and lies on the semi- 
axis 0x,. At the point x, the situation is the following: it is easy to see that all the vectors Ax, + bu, -1< u = lare 
located within a certain angle which is separated from the set S by the reference straight line of this set, which passes 
through x,. Therefore there exists an €> 0 such that in each vicinity of the point x, there can be found points for 
which all allowed trajectories originating from them emerge beyond the € -vicinity of the straight line F before they 
reach S, Thus, the point x, does not belong to C. However, the same thing is also valid for y,, since the location 
of the angle containing the vectors Ay, + bu relative to C is analogous to that described above for x,. In fact, in 
the contrary case, in view of the continuity of the variation of this angle, a point z, would exist on the segment Jp 
at which the vectors Az,+ b and Az, -b would be antiparallel; this is impossible in case 3. Therefore the point y, 
cannot be non-zero; the same reasoning is also applicable to the semi-axis Ox,. 


Thus we have arrived at the following theorem. 


Theorem. There always exists an “optimal” stable synthesizing function relative to F. It is equal tov,(x) if 
the vectors a and aA are linearly dependent; it is equal to vo(x) if aA -(Ax,+ bsign a-b)< 0, and is equal to v<(x) 
in the remaining cases, 


The plotting of the functions vg(x) is achieved in a manner which is fully analogous to the procedure described 
above for plotting the function v(x). 


The author expresses his appreciation to V. M. Popov and A. Khalanai for their interesting discussion and valu- 
able remarks. 
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Statistical criteria are defined and che performance of systems with automatic optimization by 
statistical criteria is investigated. Optimality conditions for the transfer function of a compen- 
sator are found, as well as conditions for statistical invariance and, from these, the conditions 
for non-interaction of the system. 

This investigation permits one to select the most efficient method of searching for the 
optimal regime of a system, and guides one in designing the necessary control apparatus. 


I, Formulation of the Problem 

Automatic optimization systems form an important class of self-adjusting systems [1]. The use of such systems 
allows one to solve the problem of complex automation of production processes connected with the maintenance of 
the most suitable (optimal) regime of operation. In the case of many manufactured objects, random disturbances 
play an important role, which complicates the conditions for an automatic optimization system. The question of in- 
vestigating a system of automatic optimization in the presence of random hindrances is the subject of a fair number 
of papers [1-5]. The investigation of closed systems of automatic optimization is carried through in the paper of A. 
A. Fel'dbaum (see, for example, [1]). In these papers, the question of automatic search for extrema and the stabili- 
zation of the system in the presence of random disturbances acting on the automatic optimizer is considered, and the 
error determined, as well as the optimal magnitude, for various types of disturbance probability distribution functions, 


The influence of obstacles acting on the automatic optimizer A (Fig. 1) may be significantly lowered when 
the system is optimized under statistical criteria of the form 


@=0{> M lan) Fy 


j=1i=1 


+3 DM Bz = min, 


(1) 

j=1i=1 
where M is the symbol for mathematic expectancy, # is the symbol for the functional transformation; F, the values 
of the controlled disturbances acting on the system; aj and 8j, certain weight functions; xj, the controlled variables 
of the system; and x the value of the regulated variable. By optimization of the system (Fig. 1) one understands 
minimization of the magnitude of ¢, formed by the block for the formation of the minimized quantity BFMQ ac- 
cording to formula (1). The minimization of ¢ is accomplished by the automatic optimizer A, which carries out an 
automatic search for the necessary values Z; of controlled actions on the object O, and values 1j. The computing 
device BFMQ and the optimizer A constitute the computing device for optimization of the regime of the system 
CDORS, 


The main purpose of this article is the investigation of questions conceming the automatic optimization of a 
system (Fig. 1) by different statistical criteria gy. The investigation is limited to the case of the search for the opti- 
mal regime in several variables, for a class of teady random disturbances, Methods are shown for designing an ap- 
paratus for the automatic optimization by statistical criteria. 


ll. Equations and Methods for Investigating the System 
An illustration of the method of investigation and a derivation of the equations of the system in Fig. 1 are out- 
lined for the example involving the correction of the parameters ky, T, and T, of the compensating device K,(p) of 
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a linear system with control by a disturbance F, (Fig. 2). Here Fy is the uncontrolled equivalent disturbance; K,(p), 
K,(p), Kg(p) are the corresponding transfer functions for the system; and x is the controlled quantity. The results ob- 
tained and the method of investigation are applicable to a wide class of systems, including automatic control sys- 

tems with automatic control by deviation, and combined systems. 


= 
‘ 
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We consider the problem of forming various statistical criteria g and demonstrating the system's extremal 
characteristics obtained as functions of the control variables (Z,, 2,..., 2m; Ti» Ter--+» Tp) in the presence of 
functions Hj (24, 2,- +++ Tis limiting the variables mentioned. 


* The knowledge of similar characteristics allows one then to choose the most efficient (for example, in the 
sense of rapidity in searching out the extremals) method of searching for the optimal regime, and also to select prin- 
ciples for constructing the control apparatus CDORS. 


The control systems, in which the conditions of invariance are realized, reduce, as shown in many articles (see, 
for example,[5}) to the structural scheme of Fig. 2, with two channels for the propagation of the controlled disturbing 
action F,. In the case when the transfer functions K,(p) and K(p)Ko(p) in the two channels are equal in magnitude 
and opposite in sign, complete invariance (independence) of the regulated output of variable x in relation to the 
disturbance F, occurs, If the mentioned transfer functions are not subjected to this condition, then a shift in the con- 
trolled quantity x appears at the output, this shift being due to the controlling signal Fy. Moreover, the quantity x 
shifts under the action of uncontrolled disturbances, reduced to the equivalent disturbance F,. Adjustment of the 
parameters of the compensating device with transfer function K,(p) allows one to reduce the magnitude of the statis- 
tical criterion ¢ to a minimum, In the case of a linear operator Ko(p) the variable parameters adjustable by the 
automatic optimizer A are the values of the amplification coefficient k,, the delay time rg» and the constant times 
Ty (Fig. 2). The magnitude of ¢ is formed by the computing ‘levice BFMQ in accordance with expression (1), the 
argument 1; of which may in general be corrected by the automatic optimizer A. 


We consider the definition of the criterion g as the mean square of the regulated quantity x for the class of 
stationary random input disturbances F;, assuming that uncontrolled disturbances are absent, i.e. F, = 0. In this case 


the expression 


T 


T 
= M = tim [2 (OP dt = Sp (@) | K (jo) do 


Sp (w)|Ky (jo) — Ky (ja) K; (jo) do, (2) 


is valid, where Sp(w) is the spectral density of the input disturbance, |K(jw) |is the modulus of the amplitude -phase 
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characteristic (a.p.c.) of the system; K,(w), K(jw),the a.p.c. of the system; and Kg(jw),the a.p.c. of the compensat- 
ing device. 


We consider the adjustment of the parameters of a compensating device Ky(p) with the following a.p.c.: 


K, (jo) = re (3) 


k 
K, (jo) = e— sets, (4) 
The quantity ¢ is determined from the expression 


1 (k,)? A + (k,)® k2B (@) — 2ky kak Q (w) 
1\"2 1"2 
r (@) C (@) 


dw, (5) 


A () = (1 — @*7;7,)* + + TY), 
B (@) = (1 — + (7, + 
(w) = + + (747 + (777 + (12 + + TH + 4. 


Here asterisks indicate parameters of Kg(jw), and Q is determined by the formula 


Q = (cos wt; cos wt + sin wt; sin wt) [w? (77 + T.) (7; + T3) 
+ (t — o°7T,) (1 — @*7,72)] + © (cos wt; sin wt — sin @T, COS wT) 
x + (4 — — + T,) (1 — 


where T = Ty - %+ Tp. 
We consider further particular cases obtained from the equation (5). 


1. For the case of an inertialess object with coefficients ky and ka, the compensating device is an amplifier 
with variable amplification coefficiem ky. 


Considering that T, = T,* = T, = T,* = 7, = T, = 0, and setting ky = k,* /ky®, we arrive at the expression 


= 11 Sp (w) (ky — do = Der (1) 


where DF is the dispersion of the input signal. 


An investigation of the process by which such a system moves to its optimal regime with an automatic optimi- 
zer of step type, and a formula for choosing the necessary time of averaging T and the mathematical expectation of 
the required steps of the search are given in [8]. The expressions indicated may be used for investigating the process 
by which a system moves to its optimal regime for other particular cases considered below. 


2. For the case of an inertialess object T, = Ty* = T, = Tj® = 0 in the presence of pure retardation T) = TP - T° 
with ky = kp,we find 


(kg) 
= 2k [Der — Rep (At)), 
where Ref@T) is the correlation function of the input signal AT = 7 - T. 


(8) 
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Setting Rpp(AT)= “HAT and ky = 1/ ¥2, from formula (8) we obtain 


@ = 1 — 
Approximating a portion of the correlation function by a parabola, we arrive at the expression 


= Der [i — (t1 — — %)*), (10) 
where y is the scale factor coefficient, and 1) is bounded by the values 7, - T+ 1//y. 


We find the equation of motion for the system in this case, using the expression considered in [8]. Introducing 
the notations @ = 1, - T, - T, we may write 


Tos k 
50; = “fa. Aq, dt = fi. -2A070,Derr = BD 


B = 2 A0, 


80, = = BD — = BD pr (1 — BD pr) 
— BD pr, (1 — = BD (1 — BD rp)’, 


0 [n] => 0 (i _ BD gp)". 
Substituting the value of 6, we find 


[n] = — (4 — BD + (0) (1 BD pr)". 


The condition for convergence of the search process may be written out in the form 


<1 


— BDpp| <1. (16) 


Considering the condition of boundedness, we find the equation of the search process of the optimal regime 
from its bounds in the form 


[n} = — [4 — (1 — (1 — t+ (1 — (7) 


The corresponding mathematical expectation of the number of steps in the search is determined from the 
expression 


At, 


3. For the case of an inertialess object (ky * kg) in the presence of pure cppenaen and considering that T, = 
Ty® = Ty = Ty® = 0 and ky = k,* we obtain the expression 


where 
a2 
(14) 
| or 
| 
| 
| 
: 
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cing 


= Sp (0) + — (C08 wt; cos or 
0 


+ sin wt, sin wt)]do = Der (ky — + 2kikghk (Depp — Rep (ti —1)), 


where Dpp is the dispersion, and Rep (7 - 7) the correlation function of the input signal. Thus formula (19) deter- 
mines a surface y = (ky, T%) with a minimum at definite values 7, and ky. In the case when Rep(AT) = Dppe MTR 
expression (19) may be written in the form 


= Deg — + (1 — (20) 


Figure 3 illustrates the shape of a similar surface, whose sections at ky = cons are curves described by correla- 
tion function equations of the form (9) and whose section at T, = const are parabolas analogous to the curves described 
by equations (7). The bound functions Hj are planes bounding the admissible values of ky and 7) in the sytem. 


4. For the case of an inertial object, setting 7, = T = T, = T;® = 0 in formula (5), we obtain the expression 
co 


+ 1) + 4) (21) 


Expression (21) determines the form of the surface ¢ = (Kk , m), where m = T)/T,. For Sp(w) = So, after ex- 
panding the fractional integrand, we obtain the expression 


ons 


deo 


S, 
+ Ti (kym — | — hee)? + = (kum — 


Fig. 4 illustrates the shape of the surface ¢ = ¢(ky, m), corresponding to the expression (22). The bounding 
function is represented by planes bounding the admissible values of m and ky. 


The definition of the statistical criterion g by the magnitude of the mean square of the regulated quantity is 
inefficient in those cases when the system is under the action of additional noncontrolled disturbances F, causing 
changes in the regulated quantity x. The correction of the parameters of the compensating device Ky(p) according 
to the criterion of mean square of the quantity x leads to a “yawing” of the automatic optimization system near the 
minimum point of g due to components conditioned by the influence of the uncontrolled disturbance Fy. 


The formation in the BF MQ of a statistical criterion g as the magnitude of the mutual correlation function of 
the controlled disturbances F, and the output x is more efficient in this case. 


We return again to the diagram in Fig. 2, Denoting by W,(s) and Wo9(s) and Ws,(s) the weight functions of the 
system corresponding to the channel for the controlled disturbance F,, to the channel for the controlling action 
formed by the compensating device and the second object channel, and to the channel of uncontrolled additional 
disturbance F, respectively, we write the value of the regulated quantity x as 


(0) = 8) (6) ds— Fy 8) Wea (8) ds + | — (9) ds. 


The mutual correlation function Rp, (7) of the controlled disturbance F, and the system's output x is deter- 
mined by the expression 


T 


co 
\ Wn (s) (¢ — s) Fy (t+ 


T 
+1i W, (s)ds Fi(t— (t+ dae 
—T 


W, (8) (t + 8) ds+ (s) Rep, (t + s)ds 


—oo 


(8) + 8) ds. 
(24) 


The mutual spectral density Sp, x(w) is determined by the expression 
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Srdo) = | + 8) de | W, (s) eds 
+ \ Rp.p, 8) de \ W, (s) ds 


— Rew, + 8) | Woy (25) 


—oo 


Bearing in mind that the second integral of each of the terms represents a quantity conjugate to the correspond- 
ing transfer function of the object and compensating device, we arrive at the expression 


(@) = Sr,r,(@) (jo) + () (Kj (jo) — Ke (jo)). (26) 


Since jw) = Ko(jw) Ke(jw), expression (26) may be written in the form 


Spx (@) = Spyr,(o) K; (jo) + Sp,r,(w) (Kj (jo) — Ke (jo) K; (jo)I. (27) 


The mutual correlation function Rp, (7) is determined by the expression 


(t) = \ (jo) + (Kj (jo) — Ky (jo) Kj (jo))) eda. (28) 


Equations (27) and (28) allow a definite conclusion to be made as regards fulfillment of the condition of com- 
plete invariance in a system having additional uncontrolled disturbance F,, The condition of complete invariance 
of the quantity Rp, (7) with respect to perturbations acting on the system will be fulfilled if the equation 


K; K; (j@) (j@) (29) 
Ki (jo) K, (jo) Sr,r, 


= 
is satisfied, 
Equation (29) determines the required value of the a.p.c. of the compensating device, proceeding from the 


condition of complete invariance of Rg. ,(7) with respect to F, and F,. Thus K#(jw) is the optimal a.p.c. for a linear 
control system with disturbance F, (Fig. 2). 


In the case when there is no correlation between the perturbations F, and F,, equation (29) may be written in 
the usual form 


K — i lie) 30 
0 (jo) Ks (jo) ( ) 

One may consider that the perturbation F,, correlated with F,, forms an additional channel for the basic con- 
trolled disturbance F,, determined by the a.p.c. of Kg(jw) and the mutual spectral density Sp,p(jw). Therefore in 
the diagram of Fig. 2 the parameters a.p.c. Ky(jw) will now already be determined by the parallel union of the di- 
rect and supplementary channels of the basic controlled disturbance Fy. 


Aside from the condition of complete statistical invariance (29) formulated above, one may consider the prob: 
lem of choosing definite classes of functions Kg(jw), Sp.p,(jw) and values of the argument T which will guarantee 
fulfillment of the equation 


( W, (s) Re,r,(t + 8) ds = 0. (31) 
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The fulfillment of equation (31) for definite values of the argument 7 means preservation of the condition of 
complete invariance, determined by equation (30), and the independence of Re, x(7) with respect to additional un- 
controlled disturbance F,. In other words, in this case each of the terms of expression (27) attains the value zero. 


The equation defining the statistical criterion g as the magnitude of the mutual correlation function of the 
controlled perturbation F, and the regulated quantity x is written for a system of automatic control by F, (Fig. 2) in 
the form 


9 = (0) Ks (jo) + (01K (jo) —Ko (jo) K; (jo)}) do} min, 
(32) 
where represents some function of Rp, x(7 )- 


The condition of physical realizability of the optimal a.p.c. (29) is determined by the general methods of the 
theory of optimal filters of Wiener and Kolmogorov [2]. 


We consider a concrete particular case when the fulfillment of these conditions is comparatively simple. 
1. Case of an inertialess object. Setting in equation (32) 


Ky (jo) = ky, Ky (jo) = ky, Ko(jo) =k, K; (jo) = ke 


we arrive at the following expression, determining the criterion ¢: 


Suppose » is a quadratic transformation;.then we obtain the following formula determining the surface ¢ = 
= T): 


34 
In the case when an uncontrolled disturbance F, exerts an influence on the system, one may guarantee minimi- 
zation of the chosen criterion g by accounting for the change of two parameters of the system: the coefficient of the 
compensating device k, and the argument of the correlation function T [see formula (31)]. 
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Fig. 5 illustrates the character of the surface (kp, T) obtained for correlation functions of the type e~#T* 
having a minimum on the surface (kp, T). The bounding functions Hj are represented by planes bounding the admis- 
sible values of ky and T, 


2. Case of an inertialess object with pure retardation. Setting 


K; (jo) = ks, 

K, (jo) 
K, (jo) = e—iets, 
Ky (jo) 


we arrive at the expression 


= (t) + Rep, (ti + t) — Ree, + + (35) 


(t) + Rr,p, + t) — Re,r, + + (36) 


By way of example, the shape of the surface obtained for correlation functions of exponential type is given in 
Fig. 6. This surface is constructed according to expression (36). The boundings planes are shown there also. In the 
case when the correlation functions have the form eHT*cosat, the surfaces in Figs. 3, 4, 5 and 6 take on a more 
complicated form and several minima are attained for gy. This complicates the problem of searching for the least of 
the values ¢, and leads to the required use of some automatic optimizer of a global type [1]. 


By analogous methods, one may construct formulas and surfaces for many different cases determined by expres- 
sions (1), (31) and (32) and answer the question of the physical realizability of the optimal operators obtained, Thus 
the investigation of the problem of defining a statistical criterion by the magnitude of the mutual correlation func- 
tion between the input, representing controlled disturbances, and the regulated quantity of the system allows one to 
clarify the conditions of complete statistical invariance of the criterion g with respect to the controlled and uncon- 
trolled disturbances acting on the system. The determination of the conditions of complete statistical invariance of 
the system allows one to solve the problem of non-interaction and compensation for the influence of disturbances con- 
ditioned by the control in a connected system with many inputs and outputs, since [1, 5] the conditions for non-inter- 
action follow from the conditions of invariance. 


We consider, by way of example, a system consisting of the two diagrams of Fig. 2, the output of which is the 
sum of outputs of the individual diagrams, Then, writing the disturbance, weight functions and output signals of each 
diagram as in Fig. 2, and denoting the output quantity by x, we obtain the following expression for the mutual corre- 
lation functions Rp,,.x(7) and Rp, x(2): 


Regs (1) = War (9) — (8 + %) ds 


+ Wis (8) Rr, (8 + t1)ds + [Wor (8) — (8)] Rear, (8 + ts) ds 
+ Wes (8) Rear, (8+ as, en 
Regs = (Was (8) — Woes (8 + %) ds 
+ Wes (8) (8 + %) ds + [Wir (8) — Wis (8)) Re, r,, (8 ds 
+ (8) Re, (8 + %) ds. (38) 
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For further simplification of expressions (37) and (38), one may neglect the influence of Fj,(t), since the action 
of one circuit on the other is equivalent to a similar disturbance, In this case, the necessary conditions for independ- 
ence in the construction of parameters of the compensating device of each circuit may be written in the form 


(s) — W 202 (s)] Ry, (8 + ds = 0, (39) 


(8) — Re, (8 + ds = 0. (40) 


i 


The fulfillment of conditions (39) and (40) is determined, as for (31), for definite classes of integrands and 
values of the arguments 7; and 7». 
The arguments given may be generalized to the case of n circuits, which leads to the expression 


— Woe (9) (8 + ds = 0, (41) 


i=1—oo 


where j = 1, 2,..., 0, and 64j is defined by the conditions 


= 0 for i= j, = 1 for i j. (42) 


The reasoning may also be generalized to other structures and object diagrams, including systems with negative 
feedback. 


Equations (41)—(42) are analogous in structure to the system of equations associated with ordinary autonomous 
control systems (see [1]). However, equations (41) determine conditions for non-interaction of not the regulated 
variables themselves, but rather the statistical criteria, representing the mutual correlation functions of the con- 
trolled disturbances and the quantity x. The question of the physical realizability of similar conditions, depending 
on the choice of necessary classes of integrand and values of the arguments, is, as for ordinary coupled automatic 


control systems, an independent problem. 


Ill, Certain Possibilities for Application of the Results Obtained 
| As an example of the correction of the parameters of the compensating device of a linear control system, a 
bi method is shown of defining statistical criteria. The characteristics of an automatic optimization sytem have been 
| found and illustrated by a series of examples, and also the condition of complete statistical invariance and a method 
) for determining the optimal transfer function of a compensating device have been found, as well as the condition of 
) non-interaction following from them, The investigation which has been carried out allows one to choose the most 
i efficient method of searching for the optimal regime of a system, and also of designing the necessary apparatus. 


1, Since the characteristics of the system examined are to a definite extent known, in order to hurry the proc- 
ess it is expedient to organize the search in the form of a combination of the method of extrapolation and the meth- 
od of gradients, The application of the method of gradients is most useful, as is known [1], when the system is ini 
tially not too far from the extremum point. 


2. In IAT AN SSSR an apparatus for the automatic optimization of a system of statistical criteria, named th 
computing device for optimization of the regime of a system (CDORS) [6] is described. Such a device, as is evident 
from Fig. 1, consists of a block for the formation of the minimized quantity ¢j on the basis of the statistical criterion 
taken, and an automatic many-channeled optimizer performing a search for such values of the variables of the system 
Zi and quantities Tj for which minimization of the criterion gj chosen is guaranteed. 


The application of CDORS for a continuous search is illustrated in Fig. 7, where the structural scheme of the 
control system with object (OP) having an extremal static characteristic of the form Y = -kx* is depicted. In Fig. 7 
the notation Y is used for the regulated quantity, and X for the regulating quantity, The variation of the quantity X 
may be considered slow (secular) in comparison with the variations of the quantity AX, which represents perturbation 


obeying a normal distribution law with zero mathematical expectancy. 
wae eeeeenny We suppose the control process starts with the value X not in the optimal 
regime. In this case the following expressions are valid; 


= M[AXY] = — 2kXM [AX%), 


X=m M (m Ye k). 


Since the rate of change of the controlling signal is proportional to 
the output signal of BFMQ, the equation dX/dt = —2kaM([AX*}X, 
is valid, whose solution determines the exponential decay curve for the 
Fig. 7 variation in X, i.e, the searching process for the optimal regime of the 
system converges, 


SUMMARY 


1. In this article, the questions of automatic optimization of a system by statistical criteria are examined, 
the advantages of forming the criterion by the mutual correlation function of controlled input perturbations and 
regulated output quantities are shown, and the characteristics of an automatic optimization system and conditions 
for statistical invariance of the system are found. The optimal transfer function for the example of a system of con- 
trol by perturbation with correction of the parameters of the compensating device is found, 


2. The investigation carried out here allows one to determine the most efficient method of searching for the 
optimal regime and of designing the control apparatus. 


The author thanks N. S, Taubman for reading the manuscript and N. S, Varaks for carrying out the numerical 
calculations and construction of the ¢-surfaces, 


LITERATURE CITED 


Fel‘dbaum, A. A. Computing Devices in Automatic Systems [in Russian], Fizmatgiz, 1959. 

Pugachev, V. S, Theory of Random Functions and Its Application to the Problem of Automatic Control 

{in Russian}, Fizmatgiz, 1960. 

Chelustkin, A. B. The Design and Application of Correlation Control. Automatic Control. May, 1958, 
Krasovskii, A. A. Dynamics of continuous systems of extremal control with random search signals. 

Izv. AN SSSR, OTN, Energetika i avtomatika, No. 3, 1960. 

Ulanov, G. M. Control by Perturbation (Compensation of Perturbations and Invariance) [in Russian] 
Gosénergoizdat, 1960. 

Doganovskii, S. A. Automatic correction of the parameters of compensating devices for objects with changing 
characteristics, Proceedings of the First Congress of IFAK, Vol. 1, 1960. 

Doganovskii, S, A. Automatic adjustment of the parameters of compensators in systems involving control by 
perturbation. Avtomatika i Telemekhanika, Vol. XX, No, 8, 1959. 

Doganovskii, S. A. Automatic optimization of the regime of one class of systems by means of a statistical 
quality criterion, Avtomatika i Telemekhanika, Vol, XXI, No, 8, 1960. 


ction 
pend- CDORS | 

' 

- 

= AX 
ymous 
| 
n- 
ding 1, 
tic 2. 
3. 
4, 
1,a 
» been 5 | 
nethod 
ion of 6. 
most 
1s. 
e proc- 
e meth- 8. 
s ini« 
od the 
-vident 
sriterion 
e system 
the 
Fig. 7 
atity X 
‘urbations 
7149 
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For the extremal system considered in [1], a method is described for designing a coordinate converter. 
The effects of lag elements, dry friction, backlash, and dead-zones on the stability of the search 
process are evaluated. Some considerations on the influence of random noise are presented. 


1, Coordinate Converter 

In [1] we have considered automation of the process of the space distribution of some physical quantity in a 
three-dimensional space, using an automatic search system. As the basic example, we consider the process of the 
automatic distribution of the cutting allowance for blanks of complicated parts, In the method chosen, a number 
of parameters of the mechanical parts of the setup and random effects have an appreciable influence on the search 
process, This gave rise to the need to find the conditions permitting stable operation of the system and small 
search error. 


The present article considers solutions to these questions. 


It was proved in [1] that the process of searching for the extremal of a function characterizing the position of 
a blank piece in space depends essentially on the directions of the test motions. 


In searching for the extremal of a function Q of the variables x,,..., Xp, it is frequently possible to approxi- 
mate this function in the neighborhood of the extremal by a quadratic form 


The search for the extremal is substantially facilitated if this quadratic form has a canonical form 


Dida (2) 


In this case, a single optimization cycle over each variable yj individually is sufficient. The search is thus 
substantially accelerated. The problem then arises of converting the quadratic form (1) to the canonical form (2). 
However, since the coefficients aj; are originally unknown, this transformation must be realized automatically dur- 
ing the search for the minimum itself. 


For n = 3 the transformation to the coordinates (y,, yz, ys) of a three-dimensional y-space from the coordinates 
(Xq. Xge Xs) Of a three-dimensional x-space is realized by the formula 


3 
Ye= Dj (k=4, 2, 3), 
t=1 


where a), are the angles between corresponding axes in the old and new coordinate systems. 
Since we have the relations 


2 
a) 

\ 
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COS? = 1 l=4, 2,3 
2 (4) 


>) cos Cos = 0, 
k=1 


> COS & 4% COS Hy, = 0, 
k=1 


COS COS & = 0, 
k=1 


(4°) 


the nine factors cos ay, are connected by the six relations (4) and (4°). Consequently there are three independent 
factors, It is therefore suffient to carry out the process with respect to three parameters, denoted below by q, qe, qs. 
A supplementary automatic optimizer is used for the selection of these parameters, 


A block diagram of a setup including an automatic coordinate con- 
verter is given in Fig. 1. The system consists of the object O, at whose out- 
put the function Q of the input quantities y,, yg and ys is obtained, the co- 

y ty ordinate converter CC, the basic optimizer AO,, carrying out automatic de- 
a 0 cc a0, + termination of the variables x,, xg and xs, minimizing Q, with a search 
time meter STM and the auxiliary optimizer AO,. In this circuit, aside 
41% from the basic optimization process, i.e,, minimization of the function Q, 
a second, slower optimization process takes place, consisting in the deter- 
AO, mination of those parameters q, dg, and qs for which the search time is 
minimized. After a series of parts of a single type has been worked, the 
o second optimizer automatically finds that coordinate transformation for 
which the search time is reduced, 


Fig. 1 The automatic coordinate converter consists of a circuit transforming 
the three input parameters q, qa, qs into the nine coefficients cos &y, ac- 
cording to formulae (4) and (4"), These coefficients correspond to the gain 


factors of circuits transmitting the x, (1 = 1, 2, 3). Summing the outputs of the three circuits according to formula 
(3), the quantities y, are obtained. 


2. The Effect of Inertial Elements in the System on the Search Process 
In the present system, the drive elements, the carriage holding the part, and the relay for measuring the 


direction of motion of the drive organ are, as a rule, inertial elements, Position indicators may be considered prac- 
tically inertialess. 


A block diagram of the optimization system is given in Fig. 2, Here 
O is the object of search, IN is a linear inertial part, N is the nonlinear 
inertialess part, AO is the automatic optimizer, 


1 Consider a drive organ consisting of a separately excited dc motor. 
rt] IN N , To a sufficient degree of exactness, such a motor may be described by a 
| 

J 


third-order equation [4]. The relay in the automatic optimizer introduces 
delay into the system. The stability of the search for partial extremes with 
40 respect to each of the variables with this search method depends primarily 
on the parameters of the inertial part of the system. For the search to be 
stable, the time intervals T between measurements of increments of Q must 
be correctly chosen. This question is considered in [5]. Unstable search can 
Fig. 2. also result from incorrect choice of the time intervals necessary to pass 
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from the search over the ith to the (i+ 1)st variable. Let us denote the durations of these intervals by Tj. The search 
will be stable if during the time Tj the increment to Q due to the change in x; does not exceed the increment to Q 
caused by the change in xj,,. The inertial part of the system (drive organ and carriage) for the ith channel can be 
represented, as a rule, in the form of a chain of inertial circuits with output coordinate mj and an integrating circuit 
with output coordinate xj. 


The change in pj after disconnection of the drive organ (curve 1) and 
a | ak the change in p1j4; after connection of the drive organ in the channel of 
variable xj,, (curve 2) are shown in Fig. 3. 


The value of Tj providing stable search can be found from the condi- 


Y tion 
yp, 
0 
\ 
> Here S,, S, are the areas under the curves shown in Fig. 3 (with €; = 
Where are the coefficients of the integrating networks in the 
ith and (i + 1) st channels, respectively). 


Expanding (5), we have 


T T 
— ne dt |—| dt |< 0. 


(6) 


In expression (6) nj(t) and nj+,(t) are the transfer functions of the inertial circuits in the ith and (i+ 1)st channels, 


Let the inertial part of the ith channel of the system be constituted by the series connection of an inertial cir- 
cuit with time constant 1j and an integrating circuit with coefficient €;. Then expression (6) may be rewritten in 


the form 
‘at| — 1) dt| <0. (7) 
Substituting 6, we have 
(r+ 0 *— 1) <0. (8) 


Eitis then 


a 


1—e-"|—|y—2e —2 
(9) 


From this, for the boundary of the region of stable search, y = 2.25. For smaller values of y , the search will 
be unstable, 


| According to (7) the required relationship of dQ/dxi to dQ/dxi+, can be found from the static characteristics 
of the search object. The calculation should be oriented on the best case, when the ratio is maximum, Since the 
parameters of the inertial elements in the various channels are not identical, the interval Tj should be found from 
formulas (6), (7) and the formulas presented in [5] individually for each of the channels. 
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$8, Estimating the Influence of Dead Zones 


ate Let it be necessary to find the extremal value of the function Q = Q(x) using an automatic optimizer with dead 
iscull zone 6q. Motion of the system in an incorrect direction will be detected only in a region for which the increment of 
Q(x) during the time T exceeds 6 y, i.e., 
(z) 
=A ba, 
1) and (10) 
of 
condi - At = 
From (10) and (11) we obtain 
Umax!” > (12) 
hej= From (12) we can find the values of the amplitudes xX» and the period of hunting Ty caused by the presence of 
in the &q- Let, for example, 
= az*?+c. 
Q (x) + as) 
Then, according to (12), the system will deviate from the extremal point for 
< 
Hence, the hunting amplitude 
els, 
ial cir- 
en in 
Fig. 5 
The presence in (15) of the term ptmaxT is explained by the adoption of a discrete method of measuring the 
increments, in which the time of measuring the increment does not coincide with the time of reversal of the in- 
equality sign in (12). According to (15) the hunting period 
ch will 
T, = 4 _4 ag, as) 
Pmax 2apfnaxT 
ristics 
‘ol These formulae permit finding the dependence of amplitude and period of hunting on the magnitude of the 


dead zone. In those cases where the amplitude and period of hunting (and consequently, the error and time of 
research for the optimal position) become unacceptably great due to the presence of the dead zone, the latter can 
usually be decreased at the price of a somewhat more complicated automatic optimizer circuit. 
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4. The Effect of Backlash 
As a rule, backlash arises in the mechanical transmission between the drive organ and the carriage. 


A circuit representing the search dynamics in the presence of backlash is given in Fig. 4. In the diagram we 
have: O) the search object; AO) the automatic optimizer; IN) the inertial circuit; I) the integrating circuit; B) the 
backlash circuit; N) a nonlinear element whose characteristic presents an extreme. 


Following the method adopted in [5], it is easily shown that this system can operate if the inequality 


T 
| [2n(t) — 1] dt —8>0 
0 


is respected. 


Here max is the maximum value at the output of the inertial circuit; € is the gain factor of the integrating 
circuit, 0 (t) is the transfer function of the inertial circuit; § is the magnitude of backlash. 


t 
H(t) = 1—e 


Substituting (18) in (17), we obtain 
T 


T 
a te —i— > 0. (19) 


If the inertial part consists of the series connection of a delay network with delay time tg and an integrating 
network, condition (17) may be written in the form 


T — 2ty— > 0- (20) 


From expressions (17), (19), and (20) it follows that the backlash increases the time T and thus the total time 
for the search for an optimal position. 


In [5] it was shown that in the presence of delay and the absence of backlash an extremal system is stable with 
> 2tg (21) 


where T, isthe time interval between measurements of increments to the minimized function. 
Substituting (21) in (20) we obtain a formula for determining T in the presence of backlash: 


5. The Influence of Drive Friction on the Search Process 
In the present extremal system,drive friction acts in the mechanical transmission between the drive organ and 
the carriage. : 


Let the drive organ employ a separately excited dc motor. 
For such a system the following equations are valid: 
Mg — Mg =J2, 
Mg =A sign 


(18) 
(23) 
| 
> 


am we 
) the 


(17) 


(18) 


(19) 


an and 


(23) 


M « = Cela (25) 


(26) 


Here Mg, is the drive moment; Mg is the moment of friction; w is the angular velocity of the drive; J is the 
moment of intertia reduced to the drive shaft; A is the absolute value of the moment of friction; I, is the armature 
winding current; U is the voltage across the armature winding; cy, o are constants. 


Employing (23) - (26), the equation of motion is written in the form 


I,= 


dw u m 
fr o<0, an 
JR, 


where Tem = the electromechanical time constant, A*= 


Integrating (27) for two segments of the trajectory, we obtain 


t 
— Temp 4A" for w<0, 
o= (28) 


‘em+- for @>0. 


From (5) it follows that the present extremal system can function (has hunting of the first kind) when the in- 
equality 


is observed. 
Here t, is the time corresponding to w = 0. 
We obtain the equation for the boundary of the region of hunting of the first kind by integrating (29), 


where v= T'/% Putting A’c/u= Mg, / (My max =’, where (Mg, is the starting moment of 
the drive, (30) takes the form 


1—A)je~” = 0. 


( 


Let us consider the influence of friction in a system whose inertial part contains, aside from inertial and in- 
tegrating networks, a delay element with a delay time tg. This circuit takes into account the delay occurring in 
the drive organ commutator. The motion of such a system is described by the expression 


t, 

tote 


ty T tes 
0 ‘ 
1 time 
(21) 
(24) 
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Integrating in (32), we obtain 


T 


Putting tg/Tem = 9, we write (33) in the form 


(34) 
+ = 9 


The curve v = f(A), obtained as the solution of equation (31) is given in Fig. 5. From the graph it follows 
that drive friction broadens the region of hunting of the first kind in which the extremal system can function. Other 
conditions being equal, drive friction therefore permits decreasing T, and thus the total search time. 


However it does not follow from this that in such systems it is alway useful to increase the friction, since the 
latter is connected with the need for a more powerful drive, having a larger value of T.7. 


6. Considerations of the Influence of Random Disturbances 

In [1] it was noted that the basic source of noise affecting the search process is the unevenness of the surface 
of the blank to be positioned. The effect on the search process of protuberances and indentations in the surface 
depends substantially on the choice of criteria for the search. 


In search according to criterion Q; (minimization of the maximum deviation),indentations do not as a rule 
lead to false motions of the system. Protuberances exceeding the increment of the optimized function AQ can cause 
a false direction of motion in the system; however, in ac- 
cordance with the adopted criterion, in the next search cycle 
the system leaves the protuberance. In search according to 
criterion Q, (maximization of the minimum deviation) 
protuberances have no effect on the search process; the 
presence of indentation, as a rule, does not lead to an in- 
correct result, although connected with false movement of 
the system, which increases the search time. 


The criterion Q, (minimization of the difference 
Q;—Q,) is poorer from the point of view of noise resistance, 
since in this case both protuberances and indentation affect 
the search process. Comparison of criteria Q, and Q, with 
respect to noise resistance can be done only if data on the 
statistics of protuberances and indentations on the surface of the blank are available. 


The proposed system for automatic search of an optimum position of the blank is distinguished by relatively 
high noise resistance. 


If the dimensions of the protuberances (indentations) do not exceed Ax;, the optimum position will be missed 
only in cases of very low probability, where the two indicators fall simultaneously during the search process on 
protuberances (indentations) and find themselves on opposite sides. In this case the system will hunt about the 
protuberances or indentations exactly as about an extremal point. However, even in this case it is possible to leave 
the false extremum in any other direction, as a result of the regular shift in the direction of motion. 


The permissible search time for the best position is usually given by technological conditions. 


The basic problem of investigating the effects of noise in this case consists in determining the magnitudes of 
change in each of n variables (Ax j) between measurements of increments (magnitudes of the steps in the stepping 
canes While the coanch time docs nat exceed on acceptable valus. 


This problem can be solved directly on the automatic search system, using an automatic synthesizer. 


Fig. 6 


| | 


An automatic search system including the search object O and an automatic optimizer AO, as well as an 
automatic synthesizer consisting of an automatic optimizer AO,, a search time meter STM and a hunting amplitude 
meter HAM, is shownin Fig. 6. The automatic optimizer of the synthesizer AQ, minimizes the hunting amplitude 
(and therefore the search error) by acting on the magnitude of the increments Ax; given by the automatic optimizer 
AOQ,. Search is carried out under the condition that the mean duration of the search process (for differing initial con- 
ditions) not exceed a prescribed quantity (T,)pres- A voltage H, proportional to T, is applied to AO, and a decision 
is made to carry out the search only in regions for which Ts < (Ts)pres. 


SUMMARY 


The formulas presented permit estimating the effects on the search process of individual parameters of an 
extremal system. Such estimates are necessary for the design and adjustment of systems. In more complicated 
cases, for example, estimating the influence of random noise, the influence on the search process of the simultaneous 
action of several system parameters, and others, analytical methods become practically inapplicable. 


This indicates that for the solution of these problems it is necessary to develop a special computer, which we 
have termed the automatic synthesizer. It will be described in the succeeding article of the present series. 
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The problem of optimum speed control in second-order pulsed relay (number) systems is 
investigated. The optimum control in such systems is discussed for random disturbances 
that are given in the form of a transfer matrix. An optimizing control algorithm for cer - 
tain types of matrices is determined. 


1. Optimum Control in Continuous Second-Order Pulsed Systems 

Optimum speed control in application to continuous second-order systems (see, for example, [1]) evolves down 
to the construction of an optimum control function F(s, v). Optimum control response is determined by the expression 
w= -N sign F, changing its sign when the function F converges through zero. For a system which controls a member 
described by the equations 


(1) 


at a restriction | w| <N, the optimum control function has the form 


F (sv) signy (2) 


(the finite point coincides with the origin of the coordinates). Figure 1 shows the block diagram of such a system. 
Optimum control can be effected by plotting on the phase plane an optimum switching curve (see, for example, [1] 


s=— + v* sign v. (3) 


It is readily conceived that for the restriction |w]| = N, the optimum control in such a system does not differ 
in any way from the case | w{ =N. 


As has been shown by N. N. Krasovskii [2], in the optimum control of pulsed systems it is not possible to 
limit the control response only to its extremal values, because the amplitude of control response depends upon the 
initial state of the system, assuming all possible values within the range | w]| = N. Optimum control in pulsed 
relay systems of the type | w| = N differs therefore considerably from that in pulsed systems in which the modulus 
of acceleration has a limitation of the type | w| =N. Moreover, from the engineering point of view, it is better 
to have a control response that depends upon the coordinates of a system at a given moment of time than a control 
response tha: is a function of the initial state of the system and time. According to the results of [2] such a 
dependence is not easily obtained even for second-order systems, while in the particular case of pulsed systems 
using switching curves of the type described in (3), such a dependence is not feasible at all. 


For a pulsed control system in which the controlled member is described by equation (1), it is possible to ob- 
tain a system of difference equations valid for any repetition period, by making use of the concept of a grid function 
[3]. At T = 1, this system has the form 


sin + 4] — s(n] = v [np +2). 
2 (4) 
v(n + 1) — = w [n] 


dv ds 
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(4) 


where the limitation of the value w[n] with respect to the modulus is preserved. 
Let us introduce the following notations at arbitrary values of T: 


6) 


Then the equation (4) takes the form 


u 
Tei Yeti— Yr (6) 


For small repetition periods it is possible to obtain another system of equations by employing the well-known 
relation 


dt)t—nT T . 


In this case the difference equations would have the form 


(7) 


The system (7) has been investigated in the works [4] and [5], and the obtained results can be readily applied 
to system (6). In the following, the system (6) is investigated for a restriction |u| <1. 


2. Optimum Control in a Pulsed Relay System in the Absence of Random Disturbances 


In all practical cases, the accuracy required in performing any operation or process is limited by a certain 
value €. Hence it is natural to formulate the problem of optimum control in regard to speed as a problem of 
minimum -time displacement of the image point in the system phase-space from an arbitrary initial position to the prox- 
imity of the origin of coordinates, limited by the value ze. Here it is necessary that the final velocity of the 
point be strictly equal to zero, so as to allow the image point 
to maintain its motionless state after reaching the above proximity. 


o--tg t : k When discrete equidistant values of the control response are 
used, the number of the (2n + 1) levels is defined as a function of 
* ~F(s,0) the value of ¢ by the equation 
sign cD i 


i 
Qe ’ if 3, is an integer, 
R= 
Fig. 1. Optimum second-order relay els-|+1 in all other cases *) 
(Here denotes the integral part of o.) 


If 26 = 1, then two or three discrete controkresponse levels may be used in the pulsed relay system. The 
optimum control of system (7) at two values of the control response is investigated in work [4]. In work [5], an expres- 
sion for the optimum control in such a system has been obtained at 2¢ = 1 and three levels of control response, and 
the advantages of such a method of control have been demonstrated. 


The presence of only discrete values of control response at ambiguous optimum -control characteristics permits 
the plotting of switching curves analogous to (3), which divide the system phase plane into two or three zones having 
values of u equal to +1 and -1 or +1, 0 and -1, respectively. The equation of the switching curve obtained in [4] 
for the system (7) may be written for system (6) in the following form: 


Fy) =— sign y. (9) 


In this case, the optimum control has the form 
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= — sign — F (y,)}. 


At ternary control (for three discrete control-response levels), the phase plane is divided into three parts of the 
same type as shown in [5]. The boundaries between them are the following two switching curves: 


if y<0, ae 
1), if y>0 


x= F, (y) 


1), if 


r=F,(y)= 
y<o. 


xz > F, (y) 
Fi (y) > 


= 


~ 


Fig. 2. Optimum switching curves for the systems (7) and (4) shown by the solid 
and dotted lines, respectively. 


A switching curve (and a phase trajectory) is understood to be the discrete points(of the corresponding equations) 
that are obtained at integral values of the coordinates, connected by straight lines. Here the x coordinate is assumed 
to be continuous while the y coordinate may take on only integral values. For a control response having two or three 
values, such an assumption is natural. Similarly, in the uniform division of the integral | u] <1 into 2n parts, it is 
assumed that y takes on valves that are multiples of 1/2n. In this case, it is these points that should be connected 


by straight lines. 

Optimum control when using 2n + 1 discrete levels may be constructed in the same manner as shown in [4] and 
[5] for the case of two and three levels. In performing this construction, the points having the coordinates (Xo, 0) and 
(Xo, 1/2n) are considered first, and then the obtained results are generalized to include arbitrary points of the phase 
plane. Speaking generally, optimum control is ambiguous, and among all the possible algorithms of control it is 
necessary to select the one that leads to 2n optimum switching curves. These divide the phase plane into 2n + 1 zones, 
in each of which the control-response value is constant. 


and 
Here 
a° = 0 (11) 
F,(y) >. 
b 
a 4 ~ 
2 u ‘=--] === 
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10) 


(11) 


In the case of Fig. 2d the optimum trajectory begins in an arbitrary point (x, 0) on the x axis of the phase 
plane, and intersects twice the zone of continuous dependence of u upon the system coordinates. As a result, at X9< 0, 
the optimum control has the form shown in Fig. 3. Here the dotted lines denote the moments of the first and 
second intersection of the trajectories with the | u| < 1 zone. (In this case, obviously, ¢ = 0.) An exception are the 
trajectories which begin in the points (xp = n’, Yo= 0). For such points, the optimum control has the “canonical” 


form given in [4]: 


With increasing n of a switching curve, the dividing neighboring zones approach each other so that in the 
limit defined by n - ©, the phase plane is found to be divided into three regions. In two of these the control 


response assumes limit values, while in the third region, located between them, there takes place a continuous 
dependence of the control-response value upon the system coordinates. Figure 2, a, b, c, d show the optimum 
switching curves for a number of levels equal to 2, 3, 5, 8, and ~, respectively. 


uy = — sign 


of the system. 


At a limited n (and corresponding € > 0), the optimum trajectory also intersects twice this region (x # n”), and 


u assumes discrete values in conformity with the switching curves. 


Thus the optimum control over nearly the entire duration of the transition process (with the exception of two 
steps, at the most) consists in maintaining the limit values of the control response. It should be noted thatas h- 0, 
the zone of continuous dependence of u on the system coordinates grows narrower, and the optimum control is reduced 
to the retainment of limit values for u over the duration of the entire transition process. 


The duration of the optimum process (number of K° steps) as a function of the initial deviation of the image 
point from the origin of the coordinates at zero value of the initial velocity (yg = 0) is defined in Fig. 4, where 
the dependence of K° upon X¢ is shown for the cases: ¢ = 1/2, n= 2 (thin line), n= 3, 5, 7,..., @ (thick line), 
and € = 0, n= @ (dotted line). It can be seen from the figure that the pulsed relay system having three control - 
response levels has an exceptional position among the totality of such systems. A system having three control- 
response levels provides for a transition process shorter than that in any system with an even number of levels (in 
the figure, n = 2), while a further increase in the number of levels relative to n = 3, for a given accuracy, does not 
result in decreased duration of the transition process. With decreasing ¢, the number of levels will increase in 
conformity with (8), so that the duration of the transition process will slightly increase. In the limit case, this 


dependence will take the form shown by the dotted line in Fig. 4. 


control having two and three control-response levels is discussed in [5]. A pulsed relay system with three levels of 
control response is the easiest to realize (the system with n = 2 being an exception). The following discussion deals 


mainly with a system of this type. 


Fig. 3. Optimum control u,”. Fig. 4. Dependence of K° upon the initial conditions 
in pulsed relay systems. 


Here, K° denotes the duration of the optimum transition process in a system with arbitrary values of n. In 
this case, none of the (x, yy) points are located within the region of continuous dependence of u upon the coordinates 


The realization of pulsed relay systems of 
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3. Optimum Control in the Presence of Random Disturbances 

In a control analysis of pulsed relay systems in the presence of random disturbances, the latter are readily 
described by the transfer matrix || Pij ||, which defines the probabilities of the error-induced transitions from the 
ith to the jth value of the control response. Fortuitious distortions at each step are assumed to be independent. 


The control analysis of a pulsed relay system in the presence of random disturbances of the type described by 
the matrix shown in Fig. 5a _leads to results typical for such systems. A comparison analysis of the matrices 
shown in Fig. 5b, c leads to similar results. 


PN | 0 | 0 -1 o | 1 

a b c 


Fig. 5. Pj; tansfer matrices. 


A matrix of the type shown in Fig. 5a is obtained from the following considerations: Let the control response 
at the input of system (6) be combined with a random disturbance. We denote through q the probability of un- 
distorted passage of a symbol that corresponds to any control-response level. Suppose, further, that the disturbance 
effect is reduced to an independent distortion of the control signals, so that, with the same degree of probability 
(equal to p/2),the resulting signal can resolve into a limit value of the control response. The matrices of the type 
shown in Fig. 5b, c are obtained upon similar assumptions. 


If in the absence of random disturbances, the optimum trajectory (at the selected algorithm) was uniquely 
determined in the system phase space, then, in the presence of disturbances, instead of a single trajectory being 
produced from the initial point, the latter produces a whole “bunch” of such trajectories, each of which can have 
its own probability defined by the transfer matrix. This bunch can be characterized over each mth step by a cer- 
tain point in the phase space of the system, the point coordinates being determined by the expressions 


Im = [Ex] m Um;P), 


(12) 
Ym = [Ey]m Ug, Um, 


where [Ex], and [Ey],, denote the ensemble average of the coordinate values after m steps. 


A sequence of (x;, yj) points upon the phase plane forms an “ensemble-average trajectory." The character of 
this trajectory is dependent on the algorithm of uj-production, the form of the transfer matrix, and the value of p. 


We will determine the optimum control in a pulsed relay system as follows: 


At an initial morfent of time let the image point of the system be situated in an arbitrary point of the phase 
plane with the coordinates (x9, yg). Together with the control response, at the system output arrives a random dis- 
turbance which causes control distortion in correspondence to a certain transfer matrix. It is required to satisfy the 
equality 


[Ex}x,, = [Eylx,, = 9 (13) 


for a minimum number of K,,, steps elapsed. K,, we shall denote as the duration of the optimum process, and the 
method of u;-production that provides for such duration we shall term the optimizing algorithm of control. 


For the system under discussion having disturbances described by any of the matrices given in Fig. 5 the 
optimizing algorithm of control evolves down to the application of two switching curves shifted with regard to their 
position in the absence of disturbances. 


Such a formulation of the problem differs fundamentally from the well-known formulation in which the criterion 
by for high-speed optimization is the determination of the minimum of the ensemble average of steps necessary to reach 
a given point. Let us, therefore, examine this criterion in more detail. 


Previously, we have said that in the presence of disturbances the initial point produces a whole bunch of trajec- 
tories instead of one. If under the effect of disturbances there can occur a jump in the produced command, then 
after the first step, the image point, under the effect of the disturbance and 
control, can make the transition to one of the three phase-space points. The 
totality of these points we shall denote through D,. After the second step, the 
region of possible states can already contain 9 points. Generalizing, the re- 
gion of possible states after the kth step (let us denote it through D;,) can con- 
tain a number of points less than or equal to nX, where n is the number of 
the discrete levels of the control response, i.e., the number of commands. 
Figure 6 shows a schematic of these regions. By optimum control we 


Ty understand one that provides an “overlapping” of a given point in minimum 
Fig. 6. Schematic of the op- — 
timum process on the phase plane. Speaking generally, the selection of such an algorithm does not neces- 


sarily provide for the passage of the trajectories realized after K,,, steps 
through the given point in phase space. 
The probability, P;,(i, j),of the passage of the trajectory, at the kth step, through an arbitrary point (i, j) of 
sponse the Dy region, is defined by the chosen control algorithm and, subsequently, by the set uy, ty,..., Uy. It can be 
readily conceived that as p ~ 0, the distribution Km (i, j) takes the form 


ce 
(is j) = 8 (x) 

pe 

and the problem resolves to the regular problem of high-speed optimization in the absence of random disturbances, 

discussed in the first part of this paper. In this case K,, = K*, while 
‘ 0 
if l= y = VU, 

Pri i) {0 in all other cases. 

12) Plotting of Optimum Switching Curves 


In the plotting of switching curves we shall use the same sequence of operations as was used in the work [4] 
and [5] in the absence of disturbances. 


We shall determine the expression for the maximum of the ensemble average of the value for x and y over 
p ™ steps at yp = Xp= 0. Making use of the non-dependence between control response and disturbance, we may write 
[Ey],,, in the form of the sum 


[Eylm = [Ey + LE Ay (p) 
max = max [Ey + max[EAy(p)lm, i= 1,2,.. ..m. 


Uy 
he 
The first term attains the maximum at u;= 1; i= 1, 2,..., m; the second term does not depend or u;, and 
13) is equal to pm. Finally, we have 


max[Ey],, = gm. (14) 
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For the ensemble average of the x-coordinate, we find in the same manner as above that max [Ex], is 
attained at uj=1,i=1,...,m, and equals be | 


max [Ex),,= >) max [Ey]m= 


For the maximum of the ensemble average of the quantity xm, provided x9= yg= 0, [EyJzm= 0, we get the 


expression 
max [Ez] m = gm’. 
ui (15) 


The maximum is rcached at 
u, = sign [m — i] (im, signO=1), i=1,2,..., 2m. (16) 


It is now possible to determine the quantity ky, for the points on the x axis. Let us examine the points Xx, = 
= -gn’, Yo = 0, where m is an integer. It is readily conceived that the equality (13) can be satisfied for such points 
over 2m steps when using the control defined by expression (16). For such points K,,, = 2m. Since expression (15) 
determines the maximum value of [Ex]gm, then at xg=-q{m®+1 ),2>0, the duration of the transition process must 
be larger, i.e., in this case 


K,, > 2m. 
At the same time for x = -q(m + 1)”, according to (14), 
Kn >2(m +1). 
Hence, at -qm’ > Xp = -q(m + 1), we have 
2m < Km < 2 (m + 1). 


We select within these limits the points with an x, coordinate equal to 


— q (m’?+ 1),—q (m* + 2), — (m* + m), — q(m* + m + 1), (m* + 2m). 


We shall show that by adding one step it is possible to satisfy (13) for the points X» = -qim’ +1), O0<2 sm. 
The second equality of (13) is satisfied only in the case if at the additional position r, the control response u, = 0. 
The state of the additional position (the rth step) should be selected such as to satisfy the condition 


[EZ]om+1 — [Etlem = gl, 


point. From here we find the value of r. Making use of the transfer matrix, and taking into account that [Ey], = 
we write 


(qr) + Far +1)+ Fr—1)=al (r=). 


Consequently, at r= the condition 
= 0 


| 


Similarly, for the points = — q(m?+ l), m< l<2m ,weget = 2m+ 2,7, =m-+ 1, 
r, = 1, u,, = u,, = 0. Finally we have 


+ 2, — qm (m+ 1) >2>—q(m + 1)*. 


Fig. 7. Block diagram of optimum 
pulsed relay control system. 1) Con- 
trol member; 2) algorithm of u,° 
production from 3) transfer 
matrix; 4) disturbance. 


Fig. 8. Optimum switching curves in a 


second-order relay system in the presence 
The equations (14), (15), and (17) differ from of random disturbances. 


the corresponding equations in the absence of ran- 
dom disturbances in that they have the factor q in 


front of the values x and y. Therefore, after substituting the variables — = x/q, = y/q, the further analysis can 


be carried out in the same manner as shown in [4] and [5]. On the phase plane (€, 1) the equations of the switching 


curves (10) maintain their validity. In the (x, y) coordinates there occurs a displacement of the switching curves in 
conformity with the equations 


an 


i.e., in the fourth quadrant the switching curves are defined by the expression 


z= Fy (y)=— + — 9)- 


In the second quadrant of the phase plane the switching curves are displaced in an analogous manner. 
The duration of the process at optimum control is determined by the expression 


0 
Km = K°(2). 


The presence of random effects described by other types of matrices leads to similar results, i.e., to a dis- 
placement of the switching curves. 


From the fact that after the substitution of the variables — and 1, at a selected criterion for optimality, the 
problem under discussion reduces to the well-known problem of optimum control in the absence of disturbances, it 


follows that the block diagram of a second-order optimum pulsed relay system in the presence of disturbances does 
not change (Fig. 7), the change affecting merely the algorithm of control. 


It is of interest to point out another peculiarity of the optimum pulsed relay system in the presence of random 
disturbances. In an ordinary system without disturbances, the control law, in principle, becomes apparent as soon as 
the initial conditions are specified. On the other hand, in a system with disturbances, in determining the final 
minimum variance at given initial values, it is possible to obtain only the first value of the control response. The 
following commands can be selected only after the actual state of the system following on the preceding step has 
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become known, i.e., the magnitude of the subsequent commands depends upon the realization of the disturbance 
over the first steps. In this way, the optimum control-response sequence in such a system is of a random nature, 
dependent on the concrete realization of the disturbance. 


If a predetermined control-response set, for instance (16), is used as a control device, the equality (13) will 
be satisfied as before; the variance, PKn (i, j), however, will increase. 


The foregoing canbe illustrated by the following example. Let a system described by the equations 


Sy 
Tk = Yr +>» 
Ye+i — Yr = Ux, 
0, —1}, 
be subjected to the effect of a random disturbance described by the matrix shown in Fig. 5a, where q= 0.5. 


The requirement is optimum transfer from the initial conditions x» = 2, yg = 0 to the origin of the phase-space 
coordinates. 


For such a system, K,,(2) = K°(4) = 4. To evaluate the control results, use was made of the averaging of 50 
random -disturbance realizations (the random disturbance being given from a table of random numbers). The follow- 
ing results were obtained: 


For an optimum closed system, the mean values were equal to 


Zs = 0.32, Ya = — 0.36, 
while the squared values of x and y were 


a= 7.94, = 1.20. 


Thirteen trajectories terminated in the point (0, 0). 


For a control system with a predetermined sequence of control response, irrespective of whether the disturbance 
was concretely realized or not, in conformity with (16), we have 


ty = ux (=) = (4) 


=> —u,=1. 


In this case, for the same 50 realizations of random disturbance, we have the following results, 


7, =—0.46, = —0.12, 
= 23.18, = 2.96. 


One trajectory terminates in the point (0, 0). 


From this simple example it can be seen that while the values for x and y correlate well with the theoretical 
zero values in both cases, the final variance differs widely. Optimum control provides for smaller variance. 


Optimum Relay System 
"With a decrease in the repetition period (T -> 0), the pulsed relay system resolves to an ordinary relay-type 
system. Let the transfer matrix in such a system retain its previous form, and the distortions be independent as 

before. In this case, according to the law of large numbers, at given initial conditions, in the system there will be 
realized a single trajectory which differs from the phase trajectory in a system without disturbances. 


The expression for the optimum switching curve can be obtained by means of the limit transfer on the basis of 
the formulas (19) and (5). Inthe fourth quadrant of the phase plane we have for F; 


s(nT) 1 [v* (nT) (nT) 
= = — 3, [ners — 44 
As T + 0 we get 
s= 
At the same time 
lim [F, (y) — F; (y)] = 9,. 
i.e., the switching curves coincide. 
ce 
7” Using the equation of the switching curve in the second quadrant, we get finally 
50 ‘ 
follow- — sign v. 

As q-— 0, when the disturbance increases in strength, the switching curve approaches the s axis to an unlimited 
proximity (Fig. 8). The duration, T°, of the optimum transfer process for the initial conditions s9; vg = 0, is defined 
by the expression 

[50] 
= 27/ 

In the absence of disturbances (q = 1), we arrive at the well-known formulas of optimum control in a second- 

order relay system; at N = 1, we have 
urbance s= —'/pv*signv, T= 2V%. 
SUMMARY 

Optimum control in second-order pulsed relay (number) systems can be realized by plotting switching curves 
that divide the phase plane into zones in each of which the control response is constant. The number of discrete 
control -response levels is defined by the accuracy requirements imposed upon the optimum response. 

The optimizing algorithm of control in the presence of disturbances consists of the utilization of switching 
curves shifted with respect to their position in the absence of disturbances. The magnitude of the shift increases 
with increasing strength of the disturbance. 
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THE DESCRIPTION OF MULTIDIMENSIONAL LINEAR 
SYSTEMS IN MATRIX FORM 


I, Maty4$ and Y. 


(Pardubice, Czechoslovakia) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 7, 
pp. 876-884, July, 1961 

Original article submitted November 19, 1960 


The paper derives matrix relationships that describe the behavior of multidimensional linear 
systems. These relationships coincide formally with the generally known formulas for de- 
scribing unidimensional linear systems. The results are applied in the solution of the prob- 
lem involving the synthesis of shaping filters for the generation of m stationary random 
processes. 


Statement of the Problem 
Figure 1 shows the conditional representation of a unidimensional system. The output quantity u = u(t) of the 
system is the convolution of the input quantity v = v(t) and the pulse transient response W = W(t} 


u(t) = DW 


—oo 


which shall be expressed in the form 


where the symbol * denotes convolution. 
Then we shall use the definition 
Livy =V, L{uy = U, L{W)=Y, (2) 


where L is the Laplace transform operator. Evidently Y is the transfer function of the system. After applying the 
Laplace transform we obtain from (1) 


U =VY. (3) 


Fig. 1 


If the input quantity v = v(t) is a stationary random process, then the output quantity will also be random. We 
introduce the following statistical characteristics. 


Definition 1. The correlation and mutual correlation functions © for the stationary random process v and u 
are determined from the relationships 


) of the 


(1) 


(2) 


the 


(3) 


=@=9@(t) = M{ v (t) v(t + 
= 1 = H(t) = M{ v(t)a (t+ 
= = 4 (t) = M {u (t) v(t + 2), 
Ou =) = = M {u (t) u (t + 


where M is the mathematical expectation (averaging operator). It is not difficult to prove the fact that 


(4) 


n=9+W, (6) 
y= 1+ W =W (6) 


where W(r) = W(-r). 


Definition 2. The proper and mutual spectral densities § of the stationary random processes v and u will be 
defined using the relationships 


8,, = 5 = 5S (jo) = F 
8, = H =H (jo) = F (n), 


8, = H = H (jo) = F (n), 
8,,, = G = G (jo) = F {¥), 


where F is the Fourier transform operator. 


Note. Certain authors define the correlation functions by means of the relationship 


®,,, (1) = M (t+ u(t) 


Uru (j@) = kF (Dou (t))}, 


where the constant k acquires the values 1/29, 2/# etc. (cf., for example, [2]). Under these conditions it is 
evident that all the formulas derived in this paper remain valid. 


We shall apply the operator F to Eq. (5) and (6) and obtain 


H = SY, (8) 


G =Hy 


where Y(s) = Y(-s), s = jw. 


System with several inputs and outputs (Fig. 2) shall be called multidimensional. Assume now that a multi- 


dimensional linear system C is specified whose conditional representation is given in Fig. 2. The input quantities 
Vy Var.- ++» Vp form the row-vector 


v= v(t) = (t), % (0), 


(10) 


and the column -vector 
(2) 
(t) 


( 


= 
and the spectral densities 
= 
| 
and u 
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(the prime denotes transposition). Analogously, the output quantities uy, ug, . . . , Up, form the vectors 


u=u (t) = [u, Ug (t), Um 
u’ =u (t)’ 


Um 


Fig. 2 Fig. 3 


Figure 3 shows that portion of the system C which corresponds to the kth output Wj; = W;,,(t) is the pulse tran- 
sient response, and Yj, = Y\,(s) = L { Wj, } is the transfer function for the block between the ith input and the kth 
output of the system C. The characteristics Wj}, form the pulse transient function matrix 


Making use of the matrix and vector notation cited above, it is possible to represent the multidimensional 
system C (Fig. 2) in schematic form in just the same way as we do for a unidimensional system (Fig. 1). 


Our problem is to derive the matrix relationships that describe the behavior of multidimensional systems and 
are a generalization of formulas (1), (3), (5), (6), (8) and (9), which describe the behavior of unidimensional systems. 
The derived results demonstrate that matrix calculus is a very useful means for investigating multidimensional systems. 


Matrices When They are Subjected to Operators 
Definition 3. Assume that there exist the matrices A = [Ay], B = [Bj,] and the operator K which converts the 
elements A; of the matrix A into corresponding elements B;,, of the matrix B: 


(16) 
K {Au} = By. 


(12) 
(13) 
4 
Fig. 4 
Wim 
Was +s Wan 
and the characteristics Y;,, form the transient function matrix for the system 
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(13) 


(14) 


(15) 


The effect of the operator K on the matrix A is defined by the equation 


K {A} = B. 


Thus, applying the Laplace transform operator L to the matrix W, we obtain 


L{(W)=Y. (18) 


Analogously, we obtain the following results for the Laplace transforms of the vectors of the input and output 
quantities for the system C: 


V =L {vr}, 


(19) 
U' =L{u'}, V’ =L {v’}. 


Henceforth we shall apply various forms of operators to the matrices in accordance with Definition 3. 


Systems with Definite Input Quantities 


Definition 4. The symbol A B, where A and B are matrices whose elements Aj), By, ate functions, denotes 
the matrix 


C=AxB, (20) 


whose elements are determined by the relationships 


Cx => > Au * Bix. (21) 
j 


It is obvious that the matrix C is a “product” of the matrices A, B, where the “product” of the elements denotes 
their convolution. 


The output quantities . . 


- » Uppy Of the system C can be written in the known form 
n 


m= > Wi, k=1,2,...,m, (22) 
i=1 


In accordance with Definition 4 the system (22) can be written in matrix form as follows: 


(23) 


The Laplace transform of the convolution is the product of the transforms of the functions entering into the 
convolution. Thus, applying the Laplace transform to Eq. (23), we obtain the known formula 


U = Vy. (24) 


Formulas (23) and (24),which describe the behavior of multidimensional systems with definite input signals, 
evidently coincide formally with Eq. (8) and (10), which describe the behavior of unidimensional systems. 


The Serial Connection of Multidimensional Systems 

Assume that the system B can be represented in the form of a serial connection of two multidimensional ‘ 
systems (Fig. 4); the first of these (C,) has n inputs and q outputs, and the second (C2) has g inputs and m outputs. ; 
Expressing the vectors z and Z in terms of v and V, respectively, it is possible to write the following relationships : 
on the basis of the formulas cited above: 


= ve W,, u=2*W,, Z=VY,, U =ZY;, (25) 


(12) 
a7) 
| 
u, 
Up 
| 
kth 
= | 
u= ov, W. 
and 
ems 
| systems. 
(16) 
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where W;, Y; are the corresponding matrices for the system C,; Wg, and Y2 are the corresponding matrices for the 
system Cy. Performing a substitution, we obtain 


u= U* W,* W,, VY,Ys. (26) 
If in general the system C is a serial connection of k systems Cy, G, .. . , Cy, then we can write 
u=veW,eWye... *W,y, (27) 


Making use once more of the relationships (23) and (24), we obtain the following expression for the matrices 
W, Y of the system C: 


(28 
W = W,*W,e...* Ws, 


Y = Y,. (29) 


The Correlation Functions and the Spectral Densities 
The correlation functions of the input random processes form a quadratic matrix of order n, 


where in accordance with Definition 1 


= M (0) % (¢ + 1). (31) 


Aralogously, we denote the matrix for the mutual correlation functions of the input and output quantities by 


(32) 


and the matrix for the correlation functions of the output quantities by 


Yu 


[yp 


Ny (t) = M (0) uy (¢ + 1), 
Wie (t) = M uy (¢ + 1}. 


In accordance with Definition 3, these matrices may be written as 
(t) = M {v(0)’ v(t + 


(t) = M (v(t u(t + 1), 
(t) = M u (t + 


| 
Qu--- (30) 
Pan 
n= = ter |_| 
where 
(34) 
(38) 
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For the matrices of the spectral densities S, H, G we obtain 


(39) 
S = [Su] = IF = F = F (M (0)’ v(t + 
H = (Hu) = = F = F (M (v(t) + (40) 
G = [Gu] = LF = F Op) = F (M (t + 


5) 


(41) 
7) 
Systems with Random Input Signals 
In this section we derive matrix expressions which describe the relationship between the characteristics of 
" multidimensional systems (the matrices W, Y) and between the statistical characteristics of random quantities at 
- the inputs and outputs (the matrices 9, y, p; S, H, G). 
8) a) The matrices for the correlation functions. 
9 Substituting into (37), we obtain 
n(t) =M \ + — 2) W (2) dz} = | @ —2) W (2) de (2) 
—oo —oo 
0) from (23) and (36). 
Therefore, 
1) 
v= (x — x) W (2)dzx 
The matrix 7" can be written as 
= sq’ = W' 
3) Substituting in to (44), we finally obtain 
y= W’ (46) 
Relationships (43) and (46) describe the relationship between the matrices for the correlation functions and the 
34) matrix for the pulse transient function of multidimensional linear systems. For the case n = 1, m = 1 these relation- 
ships evidently coincide with the formulas for unidimensional systems. 
5) b) The matrices for the spectral densities. 
We shall now apply the operator F to the matrix relationships (43) and (46). For the matrix of the mutual 
spectral densities we obtain 
F (n) = F (p+ W) = F (q) F (W). (47) 


H = SY. 
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In similar fashion we can write 


F (p) = F (W' «@ = F(W’) F FW). (49) 


for the matrix G. 


Therefore the matrix G for the spectral densities of the output quantities is determined by the relationship 


G =Y'SY. (50) 


Formulas (48) and (50) describe the relationship between the matrix of the transfer functions Y of multi- 
dimensional systems and the matrices for the spectral densities of the input and output random quantities. For the 


case n = 1, m = 1 these matrix relationships evidently coincide with the generally known formulas that describe 
unidimensional systems. 


By comparing the matrices in the left and right sides of (48) and (50),we obtain known relationships for the 
spectral densities Hj}, Gy, (cf., for example, [2]). Thus these simple formulas are actually complex systems of 
equations. 


Shaping Filters ) 


In investigating dynamic systems whose input quantities are random functions we often use simulators. In that 
case the simulators must be used to generate the required family of random signals with the specified statistical 
properties. 


For stationary random signals this problem can be reduced to the following: It is required to determine the 
transfer function of a linear system with m inputs for which the desired random quantities appear at m outputs when 


m uncorrelated white noise signals are applied to its input. In that case it is obvious that the matrix S for the spec- 
tral densities of the input quantities is a unit matrix. 


Thus we obtain 


ry = G, (51) 


from (50), where G is the matrix for the spectral densities of the desired processes and Y is the matrix for the transfér 
functions of the shaping filters. 


Thus the problem of generating m stationary random processes reduces to solving the matrix equation (51) for 
Y. It is not difficult to prove that Eq. (51) has an infinite set of solutions rather than one solution. 


The problem of designing the shaping filters for m stationary random processes is treated in [3],where formulas 
are derived that represent a particular solution of (51). 


In the Appendix to this paper we describe certain properties of the general solution for Eq. (51). 


The basic property of the general solution of (51) consists of the following: If Y, is a particular solution and 
- ¥ is any solution of (51), then there exist quadratic matrices N, M of order m which satisfy the conditions 


NN =E, M'GM =G 


( denotes the unit matrix), which are such that 


(52) 


Y =NY,=Y,M. (53) 
Thus the matrices N and M are matrices for transfer functions. In accordance with (52), N is the matrix for 
the transfer functions which reconverts m uncorrelated white noise signals back into m uncorrelated white noise 
signals, and M is the matrix of transfer functions which converts the input random processes with the spectral den- 
sity matrix G into output processes with the same spectral density matrix G. 


| 


Making use of formula (29), we find from (53) that the system of shaping filters can always be represented in 
the form of a serial connection of the systems N and Y, or Y; and M in that order. This property of the general so- 
lution for (51) reflects the physical meaning of formula (53). It is not difficult to prove that the matrix NY,M is also 
a solution of (51); thus it is a matrix of the transfer functions for the shaping filters. 


For n= 2 we shall cite one form of the matrix N which does not coincide with the unit matrix. Assume hy, 
hg, hg are any functions that satisfy the conditions 


hih; = 1, i = 1,2,3. 


In other words, assume h; are the transfer functions for phase-shifting filters. Then the quadratic matrix 
(54) 


evidently satisfies the first condition of (52). 
Example. The spectral density matrix G for the desired processes is specified in the form 


1 
(i—s)(2+8) Ts 


2 i 


Using the method described in [3], we determine the matrix for the transfer functions of the shaping filters 


which is a particular case of (54), we find 


{ 
2 1 
5 0 


from (53). 


It is obvious that Y2 also satisfies Eq. (51), and therefore it is a matrix of the transfer functions for the shaping 
filters which generate the desired random processes. In certain cases we can use formula (53) to simplify the matrix 
of the transfer functions for the shaping filters. 


In practice the problem of generating random processes with a singular spectral density matrix G is also im- 
portant. The problem of solving Eq. (51) with a singular matrix G is analyzed in [4]. 


If we use r to denote the rank of the matrix G (r < m), then the solution of (51) can be sought in the form of 
a matrix Y with the dimensions s, m,:where s =r. Therefore this solution should always be sought in the form of a 
matrix Y with the dimensions r, m. 
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After ordering the quantities u; and correspondingly ordering the matrix G (cf. [4]), it is possible to seek the 
particular solution of (51) by the method described in [3]. The properties of the general solution in that case are 


completely analogous to the properties cited above for the case of a nonsingular matrix. 


SUMMARY 


The input and output quantities for multidimensional systems can be represented in the form of vectors v andy, 
The pulse transient functions W;,, and the transfer functions Y;,, form the matrices W and Y. The behavior of multi- 
dimensional linear systems in matrix form is determined from formulas (23), (24). Equations (28), (29) express the 
law for serial connection of multidimensional systems. 


If the input quantities v for the system are random, then the correlation functions for the input and output 
quantities form the matrices y, 7, », and the spectral densities form the matrices S, H, G. In this paper we derive 
the matrix expressions which describe the interrelationship between these matrices and a) the matrix for the pulse 
transient functions of the system; b) the matrix for the transfer functions of the system. Thus, for example, the 
matrix G for the spectral densities of the output quantities is given by the relationship 


G=Y'SY. 


The derived matrix relationships make it possible to describe the behavior of multidimensional linear systems 
using a very simple form which is the generalization of a description of unidimensional systems. Therefore matrix 
calculus is a very useful means for investigating multidimensional linear systems. 


The derived formulas are used for solving problems in designing shaping filters for m stationary random process 
This problem can be reduced to solving the matrix equation (51). The properties of the general solution of this equa- 
tion are described in the Appendix. It is demonstrated that the system of shaping filters can always be represented 
as a serial connection of several systems. 


APPENDIX 


The Solution of the Matrix Equation Y'Y=G 
The problem of determining the transfer functions of the shaping filters reduces to solving Eq. (51) relative to 
the matrix Y; here G is a Hermitian matrix (the matrix for the spectral densities of the desired processes). 


The solution of this equation can be sought in the form of matrices with the dimensions s, m, where s = m. 
Here we shall limit ourselves to the case of quadratic matrices Y (s = m). 


It is not difficult to prove that this problem does not have one but an infinite set of solutions among quadratic 
matrices of order m. One particular solution is cited in [3] where the process of computing it is demonstrated. 


The family of matrices Y [which are solutions of Eq. (51)] and the basic properties of these matrices are de- 
scribed by the following theorem. 


Theorem. The quadratic matrix Y is a solution of the equation 


=G, (55) 

where G is a nonsingular Hermitian matrix only for the case where it can be written as 
Y = NY\M, (56) 

where Y; is any particular solution of (55) and the matrices N, M satisfy the equations 
NN=E, MGM =G (57) 


(E denotes a unit matrix). 
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The proof of this theorem is simplified by Lemmas 1, 2. 


Lemma 1. If the matrix Y, is a solution of Eq, (55), then the matrix Y = NY,M, where N, M are arbitrary 
matrices that satisfy Eq. (57), is also a solution of (55). 


Proof. From the statement of the lemma it follows that 
=G. 


By transposing the matrix Y we obtain Y' = M'Y',N’, and therefore Y'Y = M'Y',;N'NY\M, The use of conditions 
(57) converts this equation to the form (55),which proves the lemma. 


Lemma 2. Assume Y;, is a particular solution and Y is an arbitrary solution of Eq. (55). Then there exist the 
matrices N, M which satisfy (57) and are such that 


Y = NY1, 58) 
Y = (59) 


Proof. From the statement of the lemma it follows that Yj¥,= Y'Y=G. Since G is a nonsingular matrix, 
it follows that Y, and Y are also nonsingular. Therefore ¥ = (¥')~*¥',Y;. The matrix N = (¥')"*Y', evidently satisfies 
the first condition in (57). Assuming M = Y }Y, we immediately obtain (59). It is not difficult to prove that M 
satisfies (57), and thus Lemma 2 has been proved. 


By combining the statements of both lemmas we obtain the proof of the theorem cited above. 


LITERATURE CITED 


1, F. R. Gantmakher, Theory of Matrices [in Russian] (State Technical Press, Moscow, 1954). 
2. -V. S. Pugachev, The Theory of Random Functions and its Applications to Problems in Automatic Control [in 
Russian] (State Technical Press, 1957). 


3. ‘1. Maty48 and Ya. Silh4nek, “A generator which produces random process, with a specified spectral density 
matrix," Avtomatika i Telemekhanika, 21, 1 (1960). 


4. Maty&s, J. Tvarovacf filtry pro generovEnf s skupiny stacion&rfch procesu s pfedepsanfmi statis 
tickymi charakteristikami, Aplicace matematiky, No, 4 (1961). 


e 
and u. 
ulti - 
he 
Jerive 
items 
trix 

m. 
de- 

| 

(56) 

|| 

177 


THE ANALYSIS OF PERIODIC STATES DUE TO LEVEL 
QUANTIZATION OF SIGNAL IN AUTOMATIC DIGITAL SYSTEMS 


Yu. M. Korshunov 


(Ryazan') 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 7, 
pp. 885-895, July, 1961 

Original article submitted December 19, 1960 


A method is presented of determining the amplitude and the frequency of periodic states due to 
level quantization of signals in ADS, the method being based on the describing function approach. 
Graphical evaluations are shown. The accuracy of the method is estimated by comparing it with 

numerical solutions of simultaneous difference equations. 


The level quantization of the signal in a closed-loop automatic system can lead to oscillations whose mag- 
nitude may considerably exceed the quantization step magnitude. In [1] an easy method was given of determining 
the highest possible error value caused by level quantization and conditions were found under which periodic states 
are feasible in the system. In this paper a method is presented leading to the determination of such periodic states 
by using the describing function approach [2]. 


1. Equivalent Block Diagram for Analyzing Periodic States in an ADS 

One can assume when investigating periodic states in an ADS that only the output quantitity is level quantized, 
as the output quantity is taken as zero. In order to obtain more general results it is useful to take the magnitude of 
one quantization step as the unity of measurement of the signal level, and also the time interval T» of quantization 
as the unit of time measurement, i.e., to introduce the nondimensional time rt =t/T». Under these assumptions one 
is able to study the periodic states due to level quantization via an equivalent relay -pulse system shown in Fig. 1. 
There the pulse element PE produces instantaneous signal value ¢(r) at discrete time moments rf = n, that is, it 
transforms the signal (rT) into a lattice-type function e[n]. The clamper Cl complements the holding of the 
measured level in the time interval between consecutive measurements. The release of discrete levels from the 
signals €.. is effected by a relay element RE with a multistep characteristic shown in Fig. 2a. 


It is convenient for our purpose to assume that the level quantization is implemented by single-step relay 
elements. A characteristic of such an element with dead zone Ay is shown in Fig. 2b. By putting A, = k - 0.5(k= 
=1,2,...), we can represent the output quantity of a multistep relay element in the form 


G=)G., 


which is equivalent to parallel single-step relays with dead zones Ay, dg, . . . , respectively. 


It is a particular feature of the relay-pulse systems that the switching (operation) can only take place in the 
relay at integral values of r. Consequently the period T of a periodic state must be an integer. The output of 
each relay is in the form of a rectangular train pulse of unit height, the pulses being separated by blanks. The 
duration of a pulse, as well as of a blank, will again be a whole number. 


It will be assumed in the sequel that the continuous part CP of the system contains at least one integrator. 
Moreover, in the case of a periodic state the signal G(r) cannot have a constant component. We shall limit our- 
selves to the symmetric state of the operation of a single-step relay in which the width of the positive and the nega 
tive pulse is the same for each relay and is subsequently denoted by 8. The set of values of § for all single-step 
relays, together with their period T, characterizes fully the form of the periodic state, 


2. Describing Function Approach in Relay-Pulse Systems 


~ The describing function approach can be used to investigate periodic states in nonlinear systems which occur 
in the absence of the input signal: 


a(t) = 0, e(t) = — v(t). a) 


oft) 


Fig. 1 


In[3] an account was given of the application of this method to an approximation of a general periodic state 
in systems with a nonlinear pulse element. 


Let a signal 


e(t) = sin(Qr D), (3) 


ing be sent to the input of a pulse element, where Q = 2n/T is the angular frequency. The complex amplitude of the 
tes signal is 
tes 
tized, At the output of the relay element a non-sinusoidal but periodic signal 
e of G(r) is obtained whose fundamental harmonic, denoted in the sequel by IT (rT), is 
tion 
one == ysin(Qt — a) = Asin Qr B cox Qr. (5) 
* The complex amplitude of the fundamental harmonic is 
ye * = A-+ jB. (6) 
ca In view of the relation (1), the quantity F can be represented in the form 
(7) 
(1) 
the Dy = = Ay + (8) 
" is the complex amplitude of the fundamental harmonic of the G(r) signal. 
The nonlinear element is characterized by an equivalent complex gain N; we define it as a ratio of the com- 
plex amplitudes of the fundamental harmonics of the input and the output signals respectively, the signals coming 
a from and to the nonlinear part only: 
e nega” 
ep 


(9) 
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TABLE 1. Limiting Values of |e] and the Values of I for Relays with Single-Step Charac- 
teristic 
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even cos ) cos ( tel ) 
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0<e< +492 ¥ 
odd ye cos (F iat) cos (= + ist ) 


T2 4 


—02< < vei? cos (> 0 con + 


Remark. v sin for even T, 


— 


In computations it is of greater practical value to use the reciprocal of this quantity: 


(10) 


N i 


The transfer function of the system's continuous part is denoted by W(q), where the appearance of the operator 
q = pT, is due to the introduced nondimensional time r. When the sinusoidal signal I(r) is the input to the con- 
tinuous part, the output is also sinusoidal, say v(r),and its complex amplitude is 


(11) 


V TW (Q). 


The condition governing the existence of a periodic state in the case of (2) can, in view of (10) and (11), be 
written in the complex form as 


W(jQ) 


(12) 
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The easiest way to obtain the points where the condition (12) is satisfied is to find the common points in the 
complex plane of the characteristics W(jQ) and of Z. 


We examine the form of the Z characteristic in the case of a single-step relay with dead zone A. The input 
signal (3) at the moments of the pulse element operation is 


e [n} = | e| sin (n® + Q), (13) 


where © = 2n/T is the angular frequency corresponding to r = 1. One can limit the investigation to the initial 
phase lying within theband from -®/2 to +0/2, as any other value of ¢ is reduced to these values by an appropriate 
choice of the initial ordinate. 


Fig. 4. Construction of the boundary line of 
the -E domain. 


The relay output signal is shown in Fig. 3. In order that the positive and negative pulses be oflength 5, the 
following conditions must be satisfied: 


elph>A, elgl< A, elrla—dA, els) >—A. (14) 


Generally speaking, the conditions (14) will be satisfied not only for one value of |¢| and ¢, but for some 
complete intervals of these values. To each oscillation mode G;, (17) correspond two limiting values, |¢| min and 


le | max’ Which depend on A and ¢; they determine in the complex plane the admissible domain of values of £ 
whose boundaries are 


Emin= 1 \min®? =| ® (15) 


_ It should be mentioned that to facilitate the computation of the Z characteristic one should form the domain 
of -E instead of the domain of the admissible values of E. The values of |€| min and of | | ,,,, for a range 


of values of T and 6 were computed in [4] in the case of a single-step relay characteristic, A summary of these 
results is given in Table 1. 


The output signal G(r) of a relay with a multistep characteristic is determined by the entire set of values 
54, each corresponding to a separate step of the characteristic; it is known how the domain of the admissible values 
of -E should be formed for each step. The common part in the complex plane of all the domains thus obtained 
gives the domain of admissible values of -E, which determines the mode of oscillations of G(r). The domain of 
admissible values of Z obtained with the aid of the formula (10) will correlate with the latter domain. 


The domains of admissible values of Z are needed for all possible combinations of § and at the same time 
for different values of T. It has proved advisable to combine groups of domains into one general domain in order to 
reduce the number of graphs. We shall construct the general domains of Z corresponding to the operation of one, 
two, three or more steps of a multistep characteristic. These will be called the first, second, third, etc. generalized 
domains. Thus when T = 6, there are six possible combinations of § corresponding to the operation of the first three 
steps of the characteristic (when reading from left to right the digits show the number of steps when 6 equals three, 
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TABLE 2. Modes of Periodic States 
Oscillation | 
T | mode of G(T) 
7 O10 | 0.970 0° 0.704 +24°04’ 
8 1120 4.185 6°10’ 3.875 +-20°50’ 


TABLE 3. The Values of e[n] Obtained by Graphical and Com - 
putational Means 


T = 8, mode 1120 


T = 7, mode 0010 


numerica 
method 


graphical 
"| method 


numeri 
method 


40.275 | 0 +4.384 1.353 
+0.681 | +0.637 | 1 

2 40.560 | 40.598 | 2 

3 40.018 | 4.0.05 | 3 

4 —0.537 | —0.572 | 4 

5 —0.690 | —0.704 | 5 

6 —0.323 | —0.265 8 


two or one respectively): 030, 021, 120, 111, 201, 102 (see Appendix 1). Six Z domains will correspond to these 
combinations. When T = 6 the third generalized Z domain is obtained as the external boundary of these domains. 


It should be noted that the generalized Z domains, corresponding to the operation of a large number of steps of 
the characteristic, are mainly to be found within the domains corresponding to the operation of a smaller number of 
steps. As the occurrence of periodic states with a large amplitude usually proves inadmissible, one is able in practice 
to limit the investigation to the construction of a small number of generalized Z domains, say from one to five. The 
fifth generalized domain will then correspond to the operation of five or more steps on the characteristic. 


Construction of Domains for -Eand Z. 
As can be seen from Table 1, the boundaries of the domains are formed by a straight lines whose equation can 
be written in the form 


(16) 


A 


The straight line mn is shown in Fig. 4. In order to construct such a boundary straight line it is enough to 
find the points where it intersects the ox and oy axes, or the ou and ov axes, the latter forming an angle with ox of 


(x —4) and (x—{) respectively. The axis of ox is a symmetry axis for the -E domain, because the 
expressions for |¢| min and |¢ | max remain the same when replacing ¢ by -9. 


One starts by constructing first the boundary lines of the -E domain for a single-step characteristic with dead 
zone 4. Then by putting A = A; = k- 0.5, and by giving k the values 1, 2, . . . , we obtain the boundary lines for the 
-E domain with a multistep characteristic. The set of § values for any of these domains is determined by finding 
all single-step characteristic domains to which the considered domain belongs. 


When T is even, one can compute for each set of § values the quantities y' and y", which are the sums of 
Y, Values corresponding to even or odd §, respectively. When T = 4Q and 0 < @ < 9/2,we obtain from Table 1 


= 54 ="4 = = *4 I—Lep’o = "4 ‘g = 1 (8 ‘006% = 
= = = ‘9 = 1 (@ = 5A ="4 ‘THOS = = ="A ‘sg =L 


(p = ‘ze0's—009'¢ ="4 = = = *A = 1 (2 = 84 ‘oze’s = 
= *A ‘gogo = = 1 (q = SA ‘ozg’e = = = = 1 7 jo sydein “Bry 


" | 
3 2 
actice 
2 SS 
~ 
can 
pe | | \ 
ad “TS | oJ 
nding \ 4 ~ 
f 3 q “Tit | 2 
el \ a 
783 


T=y'+y'e° = ye-ia, 


In this case the expression for Z assumes the form 


— (18) 


But when we have — E = + ye = 
(19) 


= —lel 
Z ei (20) 


We note when comparing the formulas (18) and (20) that the Z domain for -9/2<¢ < 0 is the reflection in 
the ov axis of the Z domain for 0 < y< 9/2. A similar result is also valid when T = 4Q + 2. The formulas (17)- 
(20) remain the same but for an exchange of y' and y”. 


% 4 


Fig. 6. The -E domain the case for single and multistep characteristics when 
T=6. 


To obtain the Z domain from the expressions (18) and (20) calls for a considerable computational work. It 
can be reduced appreciably by the introduction of a characteristic vector 


for each -E domain. 


The transformation of the -E domain into a Z domain transforms the vector P into the vector P» which has 
unit length and is directed along the ov axis; at the same time the relative position of the vector Ps with respect to 
the Z domain remains the same as that of the vector P to the -E domain. For even T, therefore, the Z domain is 
obtained simply by taking the -E domain together with its characteristic vector to the vector Pg. It is sufficient to 
perform it for 0 < y < 9/2. The refection of these domains in the ov axis gives the Z domains for -9/2< 9 <0. Af 
example of the construction of the Z domain for even T is given in Appendix 1. 


In the case of an odd T, besides the symmetry of -E domain with respect to the ox axis already considered, 
there exists also a symmetry with respect to ou axis. The Z domains are determined in this case from the following 


expressions : 


Z=—ltleie when <o< San <0, (22) 


(23) 


an 
and therefore 
b 
| 
@y 
184 


Here the quantity y represents the sum of y,, of separate steps of a multistep characteristic. We see from 
these expressions that when T is odd, the Z domains can be obtained from the -E domains which are situated in the 
sector between the ou and the ov axes, that is,when 9/4 < ¢ < 9/2, by reflecting them in the ov axis and by mag- 


nifying them afterwards y times. An example of a Z domain construction when T is odd will be found in Appen- 
dix 2. 


The generalized domains 1-5 have been constructed by applying the above procedures for T = 2, 3, 4, 5, 6, 7, 8. 
Their graphs are shown in Fig. 5; they were obtained by taking only symmetric states of the operation of each single - 
step relay into consideration. It should be borne in mind that for odd T the cases of nonsymmetric operation of single - 
step relays are also feasible, but only such in which the areas will have equal positive and negative signals at the 
output of each multistep relay. Such possibilities involve much more complicated computations being at the same 
time fairly uncommon, i.e., the Z domains associated with them being small in size. They will be found in most 
cases within the bounds of the generalized domains constructed in Fig. 5 as, similarly to domains associated with 
symmetric operation, they are situated close to the end of the characteristic vector. 


00= 1.0514 

0b=12364 
a oc =17018 
@ 42354 


Fig. 7. The -E domains for a single and multistep characteristic when T = 5. 


4. An Example of Computation and Verification of the Accuracy of the Procedure 


The case when the transfer function of the continuous part is 


1.6 
W(q) : q(t +2.59) 


is taken as an example. 


If there is no level quantization of signals, then the system is stable [5]. 


The amplitude-phase characteristic W(jQ) drawn on the graphs of the generalized Z domains is shown in 
Fig. 5. One notices that in the case of the generalized domain 1 the periodic states can occur in the system 
when T = 7, or in the case of generalized domains 1, 2, 3, and 4 when T = 8. The highest possible values of the 


oscillation amplitude are found independently from the maximum values of y for every domain according to the 
formula 


le| = | x. 


Thus, when T = 7, for the generalized domain 1,|¢ | = 0.726 - 1.210 = 0.878. When T = 8, the largest 
oscillation amplitude corresponds to the generalized domain 4, which has Je] = 0.926 - 4.437 = 4,110, 


In order to estimate the reliability of the method, modes of periodic states were found by using the Z-domains 


gtaphs; this was performed for every possible oscillation mode of G(r). Two of the nine possible periodic states are 
given in Table 2. 
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Numerical calculations were performed for the same states by solving the system difference equations with 
the initial conditions chosen in an appropriate manner. The values of ¢[n] obtained both ways are cited in Table 3, 
One notices that the graphical means used of the describing function approach give results sufficiently close to the 

exact numerical calculations, especially in cases of large oscillation amplitudes. 


APPENDIX I 


Construction of the Z Domains when T = 6 
This is the case of T = 4Q + 2 with Q= 1. Then 9 =#/3, 9/2 =/6. We find from Table 1 the equations of 
the bounding lines of the -E domain in the case ofa single-step characteristic for values of § = 0, 1, 2 and 3 when 


The -E domains constructed by using these expressions are shown in Fig. 6a. By putting A= A,=k-0.5 
and giving k successively the values 1, 2, . . . , we obtain the graph of the -E domains in the case of a multistep 
characteristic (Fig. 6b). We notice that the -E domains represent symmetric figures, either regular hexagons or 


triangles with the side a= 0.5/cos(@/2)="8/3. Any -E domain can be fully characterized by three numbers r, m 
and 1, and they show the number of steps of the characteristic working with § = 3, 2, or 1 respectively. 


The quantity I is determined for any rml domain by the expression 


r = (rein 38 + m sinOe } 


After a few rearrangements we obtain 


n 


The characteristic vector is obtained from the expression (21), 


n 3 ¥3 n 


It is easy to verify by calculation that the vector H is the radius vector of the center of gravity of the rml 
domain. In this case the non-symmetric figures, adjacent to ox or ov axes,should be considered as the halves of the 
corresponding symmetric figures. In this case it is preferable to use vector H instead of the characteristic vector P. 
The vector H is given by the expression 


i 2. 
3° 
ibe 


When the -E domains have been transformed into Z domains, the transformed vector of the center of gravity is 


where 
186 


Wf 
that is,the center of gravity of the Z domain is situated on the ov axis at a distance of 4/3 from the origin of the 
coordinate system. The construction of the Z domains is reduced therefore to referring the -E domain to a new 
center of gravity. This operation consists in amplifying the unit of measurement y times and preserving simultane - 


ously the positions relative to the tadius-vector of the center of gravity; the latter can be characterized by the angle 
a or B = (0/2)- a. 


The kth generalized Z domain is obtained as an external boundary of analogously transformed -£ domains 
between the straight lines ayby and ay4sb+3. Constructing in this manner one can neglect the triangular -E domains 
as the corresponding Z domains are within the hexagonal domains. 


APPENDIX II 
Construction of the Z Domain when T = 5 
In this case . The possible values of § are 0,1 0r2. For t 
the equations of the boundaries of -E domains are, in the case of a singlestep characteristic, 


4 


max cos +9) 


4 
2max~ cos (5 


+9) 


The -E domains constructed in accordance with these expressions are shown in Fig. 7a, By putting A= A, 
= k-0.5, and taking successively k = 1, 2,... , we obtain the graph of the -E domains in the case of a multistep 
characteristic. The domains and their reflections in the ov axis are shown in Fig. 7b. 


The Z domains are obtained from the -£ domains by contracting the dimensions y times. The value of y is 
easily found for all values of & for each domain. Thus, for the domain 21 we have 


4 


9 
cos 7-(2sind + sin = 3.044. 
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The paper formulates the basic problems involved in an investigation of step systems in non- 
steady -state modes; the methods of solving these problems are cited, and the formulas for the 
electromagnetic torque are obtained when the boundary conditions and the optimal parameters 
are taken into account for a step system with capacitive current commutation in the motor 


winding. 


Step systems are an effective means of automatizing the most varied engineering units and technological 
processes. 


The basic problems which arise in designing a step system are the following: 
1) achieving the required speed of response, which is determined by the maximum control pulse*frequency 
for which the rotor of the motor will enter into synchronism from a rest state without losing a step; 
2) determining the dependence of the electromagnetic torque of the motor on the frequency of the control 
pulses; 
| 3) finding the limiting motor torque. 
The latter problem reduces to finding the maximum static torque which must be added to the initial torque 


corresponding to a steady-state mode without causing the rotor to arrive at the nearest position of unstable equilibrium 
(this would lead to loss of a step). 


The solution of the problems posed above can be given only after determining the dependence of the elec- 
tromagnetic torque corresponding to nonsteady-state modes on the network parameters when the boundary conditions 
are taken into account and the equations for the dynamics of the system are analyzed. 


Determining the electromagnetic torque reduces in the final analysis to determining the current in the motor 
winding in transient modes. 


In [1] it was demonstrated that the power applied to a motor winding for arbitrary current and voltage wave 


shapes can be represented in the form of the following sum,by analogy with circuit theory, when harmonic signals 
are used: the resistive component 


= rit 4 de 


and the reactive component 


i dz di \ 


In these expressions i is the current in the winding; x = wl (¢) is the instantaneous value of the winding re- 
actance;r is the resistance of the winding circuit; y is the angle of rotation of the rotor; 9= wt; w is the angular 


| 
| 
| 
| frequency of the pulses. 


In studying the energy balance in the system it is possible to neglect the reactive power drawn by the winding, 
since its average value is equal to zero and it is usually compensated in the network. 


For the majority of motor designs which are used,the winding reactance can be expressed approximately by 
he following formula when there is no saturation: 


t= — cosa, 
where Xmin» Xmax tepresent the minimum (for a = 0) and maximum (a = 180°) values of reactance; a =21y; 7 is 
the number of pairs of poles; Ax = X 4, —Xmin: 


~ 


= 


Zi 


LLY YY 


Fig. 1. The electrical circuit for a stepped drive. Lg, tg represent the inductance 
and resistance in the circuit of the rectified current; r is the resistance in the motor 
winding circuit; L(g) is the inductance of the motor winding. 


For zero boundary conditions we can write the following relationship [1] for the resistive current drawn from 
the source and the electromagnetic torque for an arbitrary supply circuit: 


2@u 
sina + 2or ’ a) 


2Azle@u® sin a 
Me ™ [A2l@ sin a-} 2or]*’ 


where u is the voltage across the winding and 9 is the angular velocity. 


The angle of rotation and the angular velocity which appear in expressions (1), (2) are determined from the 
dynamics equation for the step motor when the remaining parameters are specified. However, the substitution of 
(2) into the dynamics equation leads to a nonlinear transcendental equation which in the general case is not solvable. 
The effect of the first term in the denominator of (1) and (2) at low frequencies and at the initial instant of motor 
start is small at high frequencies and can be neglected. In steady-state and almost steady-state modes it is possible 
to replace the velocity with the independent variable 9 » w/i. The error which is caused by this is not difficult 
to determine if we assume that at the initial instant the motor develops the maximum possible torque 
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then we substitute this expression into the dynamics equation for the step motor and solve the resulting nonlinear 
classical equation relative to the velocity by one or another known method. 


For a steady-state mode and a specific supply circuit the solution of this equation is cited in [1]. We shall 
study the nonsteady~-state processes in the circuit shown in Fig.1 for the case where thyratrons or silicon tran- 
sistors with thyratron characteristics are used. 


The circuit under study is one variant of networks with capacitive commutation which have found their greateg 
practical application in stepped drives. The method of investigation used in this paper can also be applied to the 
analysis of other similar networks when their specific nature is taken into account. 


The differential equations for the current in the motor winding for a conducting and nonconducting rectifier 
can be written as 


“4 (t) + cgi"(t) + cyt’ + = Ug, <A; (3) 
(t) + + cgi = 0, A<t<7, (4) 


where i’ is the time derivative of the current; T and 2 are the period and duration of the plate current i,. 
Introducing the initial conditions,we can write the following expressions in operator form: 


ap+bp*+cpP+U, 

Pi P— Pa) P— Ps) 

ap+ bp 
(pP— Pi) (p— Pp») 


(5) 


(6) 


where for A=t2=0, py, Pe. Ps are the roots of the characteristic equation 


p? CLL, + (rCL, + + p(L+L44+ rr€)+r = 9; 
a= (+ 0) + cgi’ (+ 0) + (+0); b= (+ 0) + cgi (+ 0); = i (+0) 


and c9,C;,0 are the coefficients of Eq. (3); for T =t =X, py, pp are the roots of the characteristic equation 
prLC + rpC + 1= 0; a= i'(A) + ci (A); b = i(A); and co, cy and o are the coefficients of Eq. (4). 


The equations written above are nonlinear, since L = L(g), and in the general case are not solvable. Moreover, 
in order to obtain the current corresponding to nonsteady-state modes,it is necessary to solve a third order equation 
for each of the instantaneous circuits while successively 
“accumulating” the results. This method leads to ex- 
tremely cumbersome expressions that are difficult to analyz 


even when the value of the winding inductance is averaged. 
Our objective is achieved more rapidly if we make use ofa 
abe ser discrete Laplace transform [2]. However, the direct ap- 
sf | ad plication of the theory of pulse systems in the problem 
under study is also associated with definite difficulties. 
‘ e} 4 For the purposes of further simplication we shall study 
a b certain laws governing the network. 
Fig. 2. Instantaneous equivalent circuits for a In [1] it has been demonstrated that there are no 
current source (a) and a voltage source (b). conditions and circuit parameters for which the duration of 


the rectifier plate current can exceed 120°. Moreover, ex" 
perimental verification was given of the fact that for Ly 


| 


= L(g) the rectifier plate current can be assumed rectangular within the limits of a step. The amplitude of this 


current is determined by the frequency of the pulses and the parameters of the circuit. Here we shall cite a theoretical 
verification of this proposition. 


In order to increase the efficiency it is expedient to reduce the resistors rq and r; the effect that this produces 
on the dynamics of the system requires a special analysis. 


We shall distinguish between two extreme cases where Lg = and Lg = 0. Such a separation achieves two 
things: 1) It simplifies the general solution of the problem, 2) It is necessary to clarify the expediency of 
one or another form of “ballast” (r, rg or Ly) which assures commutation that is optimal in the sense of achieving 
the best system dynamics. 


In the first case we arrive at a current source, and in the second case we arrive at a voltage source (Fig. 2). 


For a current source in the case of zero initial conditions at the instant t = 0,when we take into account the 
equations 


i, iy, + ip, Uc = ip (r + pL), ford0st=i, 


ic = pCuc 


= + = Cuc ()— for A St ST, 
uy = pli, — Liz (A) 


expressions (5), (6) acquire the following form after certain transformations: 


in 

= per $1" 

pLCiy + Cug (A) 
PLC + pCr+i 


Correspondingly, for a voltage source when the equation 


A<t<T. 


i= 


Ug = pra (ic + + + pL) 
is taken into account, in the interval 0 =t =A 


(9) 


It is evident that for the case under study Eq. (8) remains in force in the interval 0 =t =A. We shall assume 
that the steady-state current in the windings is the basic quantity for both sources; under these conditions 


i (00) = iy (00) = iy = lim pi, (p) = lim piy (p) = = 
po po d 


where ij and i, are the currents in the motor windings in the cases of current and voltage sources, respectively. 
Assuming i, = I, i(A) = 1, uc (A) = rl, we convert Eq. (7), (8) and (9) to the form 
{ 


ij pCcL+rcC . 
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+ +1)” 


i (oo) ico)’ 
In both cases the characteristic equation is written as 
-+ 2ptv + 1 = 0, 


where for Lg = © we have 


and for Lg = 0 we have 1 — k (1 —4?) 


Fig. 3. Family of curves for the relative values of current in the motor winding. 
Making use of the family of curves plotted in Figs. 3 and 4, the functions for t, and ty can be approximated 
in the form 
ty = (2.882 + 0.5E,) t, = (6.68? + 3.28,) +, (12) 


where t, is the time required for deionization of the rectifier after quenching; t, is the time during which the cur 
rents in the commutated windings overlap; 1 is the time constant for the winding. 


(10) 
(11) 
In practice the quantity € lies within the limits 0< € <1 for k = 1. 
é=0 
| : A\ 
\ 
| | 
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From (12) we find the approximate results 


The overlap time is found from the formula 


t, = (64 + 52,8 ts. (13) 


Analogously, we can obtain the functions for the case k # 1 andL,= 0. The variation of k and € , is bounded 
by the intervals 1 > k > 0, oc >€&, > &. . The nature of the variation of the current and voltage in relative 
units for various k is clear from Fig. 5. The rate of rise of the relative voltage upp drops for a reduction in k; 

therefore, for specified £4 and v; the lead time increases. 
However, under these conditions the rate of current rise in 

May "*ral the winding decreases; this, as demonstrated in [1], de- 

a” creases the electromagnetic torque on the motor shaft. 
Moreover, for small k the operating conditions of the rectifier 
deteriorate, since at the initial instant of conduction ( for 
t + 0) the current i, = Ug/rg may appreciably exceed the 
magnitude of its average value. The decrease in the rate 
of current rise cannot be compensated by decreasing the 
capacitance without causing a deterioration in the charac- 
teristic of the system, since a reduction of C is equivalent 
to reducing k [cf. formulas (10), (11)]. 


By increasing the capacitance it is possible to accelerate 
the current rise; however, when this is done, the rate of 
current decay decreases and this is turn leads to a deteriora - 
tion of the dynamic qualities of the system. Thus for im- 
proving the dynamics of the system,the case where Ly # 0 
is more favorable. 


Such a conclusion can be explained by the presence 
of a constant energy reserve in the inductance Ly, irrespective 
of the commutation of the current and the nature of the 


processes which occur in the system. The resulting relation- 
ships make it possible to draw another important conclusion. 


In existing motors the maximum deviation of the 
inductance from the average value comprises a magnitude 
of the order of + 20%, 


Expression (13) makes it possible approximately to 
determine the magnitude of the error caused by the assump- 
tion that the quantity L maintains its average value. Thus 
for ts > Oand > 2+ 10™ (ie,, for the in- 
terval which encompasses practically the entire range that 
is characteristic of the operation of a power drive), this error 
does not exceed 10% and increases when these intervals 
are altered (cf. Fig. 6). 


tan Note that the resulting relationships for k= 1 are valid 
rectifier after it is quenched. in practice for Lg * L, since even in this case the plate cur- 


rent can be assumed independent of the angle of rotation. 


nated 
| 


In determining the optimal magnitude of the capacitance from the average value of the inductance, the 
maximum deviation of the real lead time from the specified value is equal to 4 10%, A reduction of t, compared 


to the specified value can be compensated by a corresponding increase in the design capacitance by 2 
The derived relationships make it possible to study each motor stator with its supply circuits for a three-stator 


motor (or each phase for a single-stator motor) as a single pulse system whose pulse element forms rectangular cur- 
rent pulses with a variable amplitude that depends only on the pulse frequency and the parameters of the circuit, 


7 
f 
/ 
/ 
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Fig. 5. The effect of the quantity k on the rate of rise of the 
relative voltage and current. i; is the current in the connected 
winding; i is the current in the disconnected winding; the solid 
lines represent 4(r-+r4)//U, ; the broken line represents 
if + rq)/Uq, the dash-dot line represent + r4)/U gr 


When the above is taken into account, the expression for the amplitude of the rectifier plate current becomes 


2(U,— Au) 


where Ug is the rectified voltage and Au is the arc drop voltage of the rectifiers. 
The transfer function for the continuous section of the system is found on the basis of expression (7): 
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The poles of the transfer function are equal to 


hs = — B+ = Vot—B* (>); 
he = —B+4, where Vp (8 > @e). 


The transfer function for the pulse system under study can be written as [2] 


2 p (l—x) 
Po (9;,2) or 1) 
o(%.2) he — 


xSest, 


2 
(q, = e 


41,9 


T At 
where X= 


The output quantity of the pulse system (the current in the motor winding) is equal to 


1 


where Xq is the input signal amplitude (in the given case the amplitude of the plate current) and n= 1, 2,3,.... 
The expression for the electromagnetic torque in the nonsteady -state mode, when the boundary conditions 
are taken into account,for 6 > W,, can be obtained after performing certain transformations: 


(U , — du) sina 
M, (n, = [Az (K" (0, 


1 ch d+ 
x sh [6 (n + 1 + & — X) + + + 


(n+1+e—x) 


®, 1 — 2eF chy — + 
K* (0, = 14 
x e* sh +H — for OS <i, 
x sh [(8(e—X) for 


Introducing relative quantities, we obtain 
ared P 
cur- 
where 
2x 
q=—, pat, = 
(14) 
mes where 
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In a steady ~state mode 


(U , — Au)* sina 
My @ [Az + 2(r+r,)/ 


[ 1—2echdpe™ = 


x e* sh +H. — for 


2Azl (U,— Au) sina o, 1— 2 ch 


x sh [6 (© — X) + Pi — for 
In expressions (15) and (16), 


sh 8X 


5 


B 


= arc th 
¥s i—e*chd ’ 


—B (i—x) 84 — 
= arc th — 


The expressions for the case @, > 8 are obtained by sustituting @, for jé. 


By changing the parameter ¢ continuously, we obtain the value of the torque over any time interval for a 
specified control pulse frequency and all possible network parameters. 


In order to obtain the resultant torque of the shaft of the motor in a given case,it is necessary to sum the torques 
developed by each stator separately. Here we must keep in mind the fact that in the commutation interval for the 
currents in the motor windings,the resultant torque on the 
shaft of the motor is equal to the difference between the 
torques developed by each stator separately. This can easily 


be seen from Fig. 1. 
on When the expression found for the electromagnetic 
& torque is substituted into the equation of motion, it is pos- 
0} sible to determine the speed of rotation of the rotor for any 
mode of motor operation. 
5 : R The expression (14) can be written in the form 
5 we 20% sec 
Fig. 6. The error curves which were obtained Mein, e] =< j2 [n, e} sin a. 


for the assumptions L = L(¢). 


Note that the interval = = 1 corresponds to 
lap of the current in the windings. For an approximate deter- 
mination of the average torque on the motor shaft we can make use of the value i[n, €] only in the interval 0= 
ses 


Denoting the average value of the current in this interval by i; [n,e], we obtain the following formula for 
a shaft torque which is averaged with respect to the current: 


Mein, = jn, sina, 0<e<y. (17) 


Substituting expression (17) into the dynamics equation, we arrive at the classical nonlinear transcendental 
equation in descrete form: 


Ja” + ke’ + My, = In, ejsina, 0<e<x, as) 
where J is the moment of inertia of the moving masses coupled with the rotor; k is the coefficient of friction; Mg, 
6) is the static resistance torque; a’, a" are the first and second time derivatives of the angle. 
In particular, for a large current pulse duration and Mgr= 0, the movement of the rotor in the step motor is 
described by the equation for a mathematical pendulum as the stepping process is performed: 

M, sina = Ja" + ka’; (19) 
the solution of this equation, as we know, yields attenuating fluctuations of the angle relative to the equilibrium 
position. 

In expression (19) we can assume that 
2 
vee, 
20 
where I,, is the steady-state value of the current in the winding. 
Solving Eq. (18) by known methods of numerical integration or by some other method, it is possible to deter- 
mine the speed of response (i.e., the maximum frequency which the motor can handle without losing a step) as 
well as the limiting torque which the motor can develop without 
dropping out of synchronism. 
one In particular, in order to determine the conditions governing 
I ques the appearance of resonant modes which cause the step motor to “stall” 
le when the resistance forces which limit the maximum rate of operation 
of the step motor are taken into account, it is possible to use [3]. 
sly For a constant supply source and constant network parameters the 
motor has a dropping characteristic M, = f(w). In practice it is possible 
I=. to achieve any slope of the characteristic M, = f(w) for various fre- 
quencies by changing the supply voltage correspondingly. Especially 
:- LAW —. great possibilities in this regard are offered by the networks described 
ry {omm i C in [1, 5],where a choke instead of a resistance is used as the “commuta- 
———s 8 tion ballast.” This conclusion derives from the fact that the choke 
12 AAA! 54 produces virtually no voltage losses and it is therefore possible to 
BR produce a sufficient voltage margin at low frequencies; this margin 
Pi c is realized when the frequency is increased, The derived expressions 
wm for torque permit easy determination of the required law governing 
cel the variation of this voltage for variable frequencies in order to achieve 
exe the specified characteristic. 
In that case the de supply for the amplifier can be expediently 
= — realized using the networks in [4]. These networks assure a high power 
coefficient, a lighter mode of operation for the rectifiers, and con- 
for tinuous current for large regulation. 
A variant of one of these networks is shown in figure 7. The experimental verification of the correctness of our 
theoretical assumptions is given in [1]. 
The author thanks A. A. Bulgakov and Ya.Z. Tsypkin for their useful advice. 
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PRECISION TRANSISTOR INTEGRATOR 


Yu. L. Kurkin and N. S. Kurkina 


Translated from Avtomatika i Telemekhanika, Vol. 22, No. 7, 
pp. 907-913, July, 1961 
Original article submitted January 5, 1961 


A circuit for an integrator composed of transistors is proposed. This circuit is shown to be in 
better agreement with transistor qualifications for amplifying devices than the customary 
Miller circuit, and, consequently, provides better utilization of the transistors. 

A computation of the integrator’s operation is given, an expression for the transfer func- 
tion is derived, and limitations imposed by the finite time constant of the amplifier are pointed 
out. 

Errors introduced by the diodes used for protection against circuit overloading are dis- 
cussed. 


In modern electronic prototype and simulator devices, integration of potentials over time intervals is widely 
used. The basic means usually employed for integration consists of an RC circuit loop (Fig. 1a). An RC circuit, 
however, does not provide sufficiently accurate integration, because there is a retroactive action of the output volt- 
age across the capacitor C upon the input current of the circuit. This is clear from the equation 


Therefore a circuit for improvement of the integrating properties of the original RC circuit by means of a 
voltage amplifier with a high gain (K,) has been proposed. This scheme is called the Miller circuit (Fig. 1b) and 
is described by the equation 


t 
i 
Vou + \ at = — 


It is evident from this equation that 


t 
i 
0 


In Fig. 1c an integrator circuit is shown, using an amplifier with a voltage gain of Ky = 1. 


The principle of this circuit has been known for quite a long time, but,employing vacuum tubes, it has not 
obtained widespread application because without providing any special advantages it requires an additional separate 
power supply. When transistors are used in the circuit, this shortcoming is less important [1, 2], and the advantages 
are appreciable since with a gain of Ky = 1 it is easy to obtain a high input resistance rj » f,. 


Transistors by their nature are current amplifiers and have low input resistances. A high input resistance r; 
is obtained by introduction of a series-parallel feedback at the expense of only a low voltage amplification gain 
Ky. Therefore, when using transistors it is disadvantageous to build a Miller integrator according to the customary 
vacuum tube circuit, requiring use of the highest possible values for the gain K, of the order of 10° to 10°, A 
transistor integrator designed according to the block diagram of Fig. 1c is more acceptable, since with Ky = 1 it 
is easy to obtain a high input resistance (rj » f¢). 
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Fig. 1. a) Quasiintegrating circuit; b) Miller integrator (Ky - ©); c) integrator with 
voltage amplification gain Ky = 1. 


Fig. 2. Amplifier circuits: a) With gain Ky=1; b) with K, =1. 


It is essential that the obtained high input resistance be determined only by the resistance R, of the feedback (Fig. 2) 
and by the current amplification in the circuit, and does not depend upon the small input resistances of the tran- 
sistors which vary with temperature and other conditions. It is well known that current amplification is one of the 
most stable parameters of a transistor. The circuit of an amplifier having a K, = 1 and a high input resistance of 
100 meg is presented in Fig. 2a. The integration error in a device according to Fig. 1b is proportional to 1/Ky 
and in a device according to Fig. 1c the error depends upon how close the condition K, = 1 is fulfilled. Since 
the transistor is basically a current amplifying device, the use of some small minimum of current in the input to 4 
‘transistor integrator is unavoidable; this current is the larger, the smaller the amplifier gain (all other conditions 
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being equal). The presence of this current is equivalent to the connection of a leakage resistance r; in parallel with 
the integrating capacitor (Fig. 1c) and brings about an error similar to that of the Miller integrator. In the Miller 
integrator this source of error can not be eliminated. In the scheme proposed here (Fig. 1c), however, the leakage 


through r; can be compensated for by a small increase of the voltage gain Ky in the amplifier. 


The amplifier whose circuit diagram is given in Fig. 2b provides stable operation with any predetermined 
voltage gain K,,; among others, values of Ky slightly higher than one are also feasible, and these are the ones that 
are necessary to compensate for the discharge of the condenser, as is evident from the following equations: 


aR, agR, 


+ 


For the circuit of Fig. 1c with leakage it can be shown that 


that is,the time constant of the integrator is 


The needed value of K,, = 1+ — can be obtained by proper choice of the ratio Rg/Re. 
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Fig. 4. Input resistance of the amplifier 
amplifier. versus input signal voltage. 


Let us point out that the circuit remains stable even when a little overcompensated. In that case the integra - 


tion error changes its sign with respect to the error of the Miller integrator. 
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The value of K,, must not change throughout the entire range of operational voltages, that is, the amplifier 
must possess a high degree of linearity. The amplitude characteristics of an amplifier built according to Fig. 2a 


is presented in Fig, 3. 


It is not difficult to see from Fig. 1c that the impedance of the device looking in at the points A and B 
equals zero within the entire working range of voltages: the circuit at these points possesses infinite capacitance. 
Certain characteristics of practical interest for this circuit are shown in Fig. 4 and Fig. 5. 


Thus this integrator operates with positive feedback and is distinguished from known integrators of this type 
by the fact that as an amplifier component a three transistor current amplifier with a K, = 1 is used [3]. A circuit 
diagram of a complete integrator of this type is shown in Fig. 6. 


Transistors T;, Tz, and T; (Fig. 2b), all of type 1114, form a current amplifier having a nonlinear temperature 
characteristic; this is compensated for by transistor Ty. The standard operational characteristics of the consecutive 
amplifier stages are obtained by use of the proper resistances R,, Re, and Rg. By proper choice of Rg, an adjustment 
of the integrator to zero is also obtained. It is essential that R, and Ry be connected not in the power supply 
line, but in the collector loops of the respective transistors: it is not difficult to see that insertion of R, in the power 
supply line would be equivalent to connecting it in parallel with the integrating capacitance C. The potential in- 
crements at the collectors of T,; and T; during the integration process are practically equal to zero because the col- 
lectors of T, and T; are connected to zero potential (ground) through the low ohmage input circuits of the subsequent 
transistors. Therefore, R, and Rg do not limit the time constant of the transistor integrator. 


Obviously, it is possible first to consider the circuit (Fig. 6b) for R,; > », and then account for any Ry; by 
sustituting ReRti/ (Re + Ry) for Re. We thus obtain for the circuit of Fig. 6b the following: 


(1 + pCR,) 
R 1+pCR, 


C= 


From (1) we derive the following: 


= 
P) = i 


P+ (K+ICR, 


t 
(K+1) om, | (3) 


There is an approximation formula for the transfer function H(t) when t is within the usual operating range, 
that is,when t « (K + 1)CR,. This formula is as follows: 


K t 1 eK 


Let us point out that to a first approximation this does not depend upon Re, thatis, upon the load. 
The relative error during time t is 


2 CRY(K+1)* 


From this we derive 
| 
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1) 


3) 


4) 


5) 


In this way, when the acceptable error is given, the maximum of permissible integration time can be cal- 
culated. It is not difficult to see (Fig. 6b) that the possible increase of Re is limited by the output resistance of the 


current amplifier,which according to the circuit diagram 
is connected in parallel to R,. 


R,meg at U,. =3v The drift of the integrator can be estimated when 
100 the drift of the current amplifier is known (ig; in Fig. 6b). 
Solving equation (1) taking into consideration iy,,we obtain 
Fig. 5. Input resistance of the amplifier versus out RTH \R + 
temperature. 


Thus iy, directly enters into the expression under 
the integration sign along with U;p/R, thatis, ig, limits 
the maximum feasible value of R. 


Characteristics of such an integrator obtained in practice are given in the following table. 


Range of volt- Zero drift in Decay of charge in | Characteristic time | Integration error 


60 sec with maxi- | constant of integrator] for 60 sec in tem- 
Foci 8 Mag eae mum of #5 v in with C= 8 uf, sec | perature range of 
oe output, % of initial 15 to 45°C, in % 


+5 to -5 +15, slow drift | 2-5; with 24.5 v in 1200 0.29 
20 output,decay is only 
0.3% 


Let us point out that when output voltages exceed #4.5 v, the linearity of the amplifier is broken, and with 
an output of +5 v the internal discharge (decay) of the integrator already reaches 2-5% for a 60 sec period. 


Considerably better results can be obtained when a four transistor amplifier is used [4] which has a very high 
input resistance (of the order of 1000 megohms) and a voltage amplification gain exactly equal to unity. 


It should be remembered that the figures of the above table pertain to the use of the integrator without any 
special means for stabilization of the zero; that is, these figures indicate only minimum values basically feasible 
with this scheme. 


With the circuit of Fig. 6, practically any customary method for zero stabilization may be employed [5], 
and this will result in improved precision. 


The upper frequency limit of the integrator is determined by the frequency characteristics of the transistor 
amplifier. 


Since high precision of integration is required, it is reasonable to examine the problem in a simplified manner 
for frequencies at which the falling off of the transistors’ effectiveness has only started. Then the transistor can be 
regarded simply as an inertial link with the following frequency characteristics: 


(@ < @). 


Here w; = w,,(1 - @) is the cutoff frequency of the transistor with grounded emitter, and w, its cutoff frequency 
with grounded base. 
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It is known from automatic regulation theory that in a system composed of three stages possessing inertia, the 
critical gain is the higher, the more widely the time constants of the stages differ from each other. In practice it 


is sufficient to take two equal transistors and choose the third one with a much higher frequency characteristic, that 
is, to fulfill the condition ws >» w,;. For a frequency band near w, the transfer function of the amplifier can be 
written down as follows: 


(6) 


Fig. 6. a) Circuit of complete integrator; b) block diagram of in- 


tegrator; c) simplified block diagram with K; representing current 
amplification of three transistor amplifier. 


Substituting (6) into (2) we obtain the transfer function for the complete integrator in this case: 


Ke} 
KC 


| | |p 


Applying to the input a unit step, we obtain the transfer function for this case: 


Y (p) - 


t 
H (t) = 


20s V Kt) sin (o, V Kt) 
@,CR, (K +1) (K +1) VK |. 
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The first member in the curly brackets is the usual transfer function of the integrator when the amplifier has 
an infinite pass band. The second member represents the damped oscillations which are superimposed upon the in- 
tegration curve when the pass band of the amplifier is limited by the frequency w. 


For protection of the transistor integrator against overloads, silicon diodes are used which are connected in 
parallel with the integrating capacitance (lower left of Fig. 6a). 


When the signal at the integrator output exceeds a certain predetermined value, one of the diodes becomes 
transmitting. The integrator then is changed into an amplifier having only small gain, and the overload is prevented. 


The capacitance of the silicon diodes introduces a perceptible nonlinearity, since it changes with applied volt- 


age according to the formula 
qEN 
SV 8x (Ty ' 


where S is the junction area, N the density of the donors in the base, and , the contact potential difference. 


When a required precision of integration is indicated a priori, this effect sets a lower limit to the feasible 
integrator capacitance C, since we must have C > Cy/m, where m is the maximum admissible relative error. 


Leakage and polarization of integrator condensers have been neglected, since special condensers of type PGS 
were used in which these effects have been reduced to a minimum. 


In conclusion we shall point out that the described connection of the amplifier in the scheme makes it possible 
to carry out also other linear operations with given transfer functions. For these purposes, instead of the integrator 
capacitance, suitable combinations of resistors and capacitances are connected in the circuit. 
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SEMICONDUCTOR DEVICES 
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Translated from Avtomatika i Telemekhanika, Vol. 22, No. 7, 
pp. 914-918, July, 1961 

Original article submitted October 3, 1960 


A method is presented for finding the "reversal" conditions of semiconductor relay devices 
from a matrix of the equivalent circuit. The application of the method is shown by an 
analysis of specific circuits of semiconductor relay devices. Expressions for the "reversal" 
conditions are derived in terms of a system of parameters from the fourpole amplifier ad- 
mittances, parameters of the load, of the signal source, and of the positive feedback net- 
work between the collector of the output and the base of the input transistors. 


1. General Remarks 

The problem of finding the "reversal" conditions of semiconductor relay devices leads, in many cases, to the 
determination of the characteristic equations of the system and to the subsequent analysis of stability with the help, 
for example, of the Routh-Hurwitz criterion. The characteristic equation can be found either by the method presented 
in [1] or by the method presented below, for which theoretical bases are included in [2]. 


The essence of the method is as follows. An impedance or admittance matrix is constructed in a canonical co- 
ordinate system for the equivalent circuit corresponding to the conditions of equilibrium. A characteristic equation 
is obtained by equating to zero the determinant of the matrix (see the demonstration in the Appendix). Further, 
criteria for stability are used. In a number of cases it is desirable to get an expression for the reversal condition 
from the fourpole parameters of the amplifier, the parameters of the load, of the signal source, and of the feedback. 
For the solution of this problem, it is expedient to use the method of separating coordinates with the use of the 
matrix device [2], thus permitting the fourpole parameters to be found. This method is convenient because the 
fourpole parameters of an amplifier can often be found experimentally. The method is examined below in exam- 
ples of the analysis of specific circuits of semiconductor relay devices. 


2. Multivibrator with a Collector-Base RC Connection 

The general circuit of a multivibrator is shown in Fig. 1. Here Ry and Rp are, respectively, the resistances 
of the control signal source and of the load; Crp is the feedback capacitor (in the present case Rpp = ©); T, and Ty 
are the initial and final transistors of a nonlinear amplifier. The equivalent circuit of the relay device is shown in 
Fig. 2. Here Cpy and Cy are, respectively, parasitic capacities at the input and output of the amplifier; T;,2 is a 
transistor in which the parameters of transistors T; and Tz are embodied. 


In order to construct the admittance matrix of the circuit in Fig. 2, we will make the following simplifying 
assumptions: 

a) We will neglect the delay of transistor T, with respect to the delay of transistor T;,since the power of 
transistor T, is usually less than the power of T2. 

b) Wewill assume that the parameters Y2) ana ¥2) of the first and second transistors are equal to zero. This 
is permissible since in general « Yo; and « =1/Rey, Yo = 1/Re), where are, respectively, 
the output admittances of the first and second transistors when the input is short circuited . 


c) Wewill assume that the amplifier has the properties of a detector, i.e., v2), y?) 0. 


d) Wewill neglect the capacity of the collector junction, assuming that the delay of transistor T, depends 
only on the diffusion process and is characterized by the time constant Tg. 


ack. 


vely 


Taking into consideration the above mentioned assumptions, the circuit matrix in a canonical coordinate 
system will have the form 


{ 2 


2) Yew + (Cop + 
y? 
where Yu = = Ya = = % rye’ = = and fate 


sistor T, and the index 2 to Ty. 


-€, 

y Ne 

Fig. 4 


During the determination of the parameters of transistor T, in the active region at the saturation limit (cutoff) 
of the first stage, the parameters of transistor T, are established in the active region at the cutoff limit (saturation) 
of the second stage. As studies have shown, in the majority of cases it is sufficient to fulfill the “reversal” condi- 
tions obtained below for the limiting parts of the active region. 


For the performance of a circuit with alloy junction transistors,we can assume Cpy, Cpp w 0. The character - 
istic equation obtained by equating to zero the determinant of tice matrix (1) will in this case have the form 


A (p) = + Ys + Yu) p* (2) 
+ + Yu) +Cpp (¥1 + + Yu — =0. 


Applying the Routh-Hurwitz criterion to equation (2),we get the condition for instability corresponding to self- 
excitation (reversal) of the multivibrator: 


Yu + t+ Yu @) 


For the perfornrance of a circuit with high frequency transistors,we can assume tg = 0. The characteristic 
equation will in this case have the form 
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(p) = (Crp + Cpe) + pel + + Cy + Yu) 
+ Crp + Ya + Yu-— Ya) p+ + Yu) = 0. 


The condition of self-excitation of the multivibrator is found from the inequality 


(4) 


C, 
Yau > + + Yat Yu 65) 


For a discontinuous treatment * tT, - 0, Cys, Cpe 0, Crp > and the condition of self-excitation will 
be of the form 


Ya Yu. (6) 


Expression (6) can be obtained from the inequality 
A’ < 0, (7) 


where A’ is the determinant of the circuit matrix for the equivalent circuit of Fig. 2 after short circuiting the 
feedback capacitor. 


In order to find the determinant in Eq. (7), it is necessary in the end result to realize the transition in the 
limit Ypg- «. In finding the conditions for “soft” self-excitation in a multivibrator with two feedback condensen, 
there is no need for the transition in the limit. 


3. Bistable Semiconductor Relay 
A bistable semiconductor relay operating from a continuously varying input signal is obtainable from the cir- 
cuit of Fig. 1 if we make Crp = 0. With the assumptions mentioned above the circuit matrix is obtained from 


Eq, (1) after equating Crp to zero. 


For the performance of a circuit with alloy junction transistors (Cp;, Cpz ¥ 0),the characteristic equation will 
have the form 


+ + + Yes Ip + (Yo+ Yep + (8) 
— YoY rp = 0. 


The reversal condition is found from the inequality 


YoYpp > Yo¥ pp + + Yur) (Vo + Yop). (9) 


Similarly it can be shown that in the performance of a circuit with high frequency transistors, the reversal 
condition will be given by expression (9). Expression (9) can be obtained from the inequality A"< 0, where A" 
is the determinant of the circuit matrix obtained from the equivalent circuit (Fig. 2) with Crp = 0,when the delay 
of the transistors and the parasitic capacitors are neglected. 


We note that from Eq. (9) it follows that the reversal condition is indepndent of the transistor delay. This 
result can be considered as the theoretical basis of the calculation method reported by the author in [3], where the 
delay of the transistors is not considered. From the other direction, experimental studies described in [3] appear 
to confirm the conclusion indicated. For Crp #0 and tg = 0, as follows from an analysis of the coefficients of 
the complete characteristic equation, this conclusion is wrong if Crp exceeds the value given by the expression 


* By discontinuous treatment an idealized condition is implied in which a circuit “reversal” occurs instantaneously, 
and this is possible when small parasitic circuit elements are neglected. 


4) 


5) 


(8) 


max_ Yep + (12) 
SUMMARY 


The method presented for investigating the reversal conditions of semiconductor relay devices makes it very 
easy to find relations among the circuit parameters from the matrix of the equivalent circuit. The method can be 
extended also to relay devices built of other amplifying elements. The relations derived can be used to analyze 
relay devices having more complicated amplifiers. For this the values of the fourpole admittance parameters of the 
amplifier are required. The latter can be found experimentally. 


In conclusion the author takes this opportunity to express thanks to B. S. Sotskov for valuable advice while doing 
this work. 


APPENDIX 


Development of the Characteristic Circuit Equations 


We will show that for the determination of the “reversal* conditions of semiconductor relay devices, it is suf- 
ficient to investigate with the Routh-Hurwitz criterion the equation obtained by equating to zero the determinant 
of the admittance matrix of the circuit in a canonical coordinate system. In the generalized scheme of a relay 
device shown in Fig. 3 it is always possible to refer both the feedback impedance Zp, and the load impedance 
Z, to an internal fourpole parameter, and to obtain the scheme shown in Fig. 4. 


Let the expression for the transfer function of the system have the form 


(D 
A 
K 
Upon equating the denominator to zero,we get the characteristic equation of the system 
B(p) = 0. qi) 


If only one of the roots of Eq.(II) has a positive real part, the system will be unstable. In relay systems under 
the specified conditions there develops,in this case,an avalanche process of “reversal."* The presence.of a root with 
a positive real part can be detected, for example, by the use of the Routh-Hurwitz criterion. Equation (II) can be 
obtained immediately by equating to zero the admittance matrix of the circuit. 


Actually, the admittance matrix of the circuit in Fig. 4 with the terminal connected to the base has the 
form 


Yio 
2 Yu + Yoo 


The voltage transfer function at no load is written in the form** [2] 


Ko 


* In the general case the instability condition can correspond not only to an avalanche reversal but to an oscillatory 
type of transfer from one state to the other, a composite (avalanche with oscillations in certain regions)process and 
80 forth. 

** When determining the transfer function, the point of connection of the signal source is the input, and the collec- 
tor of transistor T, is the output. 
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where Ajs(p) and Ags(p) are determinants obtained from the circuit of Fig. 4 after crossing out the appropriate 
rows and columns. Thus, in order to find the reversal condition, it is necessary to examine the equation Ag3(p) = 0, 


The determinant Ags(p) is obtained from the basic matrix of the circuit in Fig. 4: 
i 3 


—Y;, 


by crossing out the third row and the third column. Yet connecting terminal 3 to the base 0 also leads to the 
diminished (abridged) matrix obtained by crossing out the third row and the third column. And connecting terminal 


3 to the base leads to the basic circuit of the relay device; consequently, the matrix of the latter is identical with 
the Ags(p). 


Similar demonstrations can be carried out for multipole schemes (more than two poles) and also for an im- 
pedance matrix in a canonical system of coordinates. 


LITERATURE CITED a 
£. V. Zelyakh, On the Stability Analysis of Vacuum Tube Circuits and of Transistor Circuits, Electrosvyaz , 
No. 7 (1960). 
V. P. Sigorskii, Methods of Analysis of Electric Circuits with Multipole Elements. Izd-vo AN USSR (1958), 
S. V. Kulikov, Design of a Semiconductor Relay. Radiotekhnika, No. 4 (1960). 


| Yu + | 
| 
1, 
2. 
3. 
810 


with 
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The paper determines the relationships between core parameters when saturable -reactor magnetic 
logic elements have a maximum gain under conditions of maximum power output and limitations 
on the amount of heating; the conditions required for a stable characteristic are taken into account. 


Three basic saturable -reactor logic repeaters are quantitatively compared for operation at 
50 and 400 cps. 


Among the most common two-cycle logic elements, it is expedient to isolate two groups of networks. The 
first of these includes networks in which a reversal of the magnetic state of the cores in the controlling element 
causes a simultaneous reversal of the magnetic state of the controlled element. Such elements perform the function 
of repeaters. The other group of networks includes those where the reversal of the magnetic state of the core in the 
controlling element does not permit reversal of the magnetic state of the core in the controlled element; the non- 
reversal of its magnetic state, which is due to saturation, causes a reversal of the magnetic state of the core in the 
controlled element. These networks can be called networks with self-exclusive reversal of the magnetic state of 

the cores. They perform the function of inverters. 


One of the basic parameters which characterize magnetic logic elements is their gain (i.e., the maximum 
number of elements which can be controlled by one logic element). 


The selection of all the parameters in a logic element network must be based on obtaining the specified 
number of controlled elements for stable operation with any number of elements within limits up to the maximum 
number, while simultaneously taking into account the requirements of economy and small size. 


This paper is devoted to determining the relationships between the core parameters and gain of elements and 
to selecting the optimal parameters of <lement networks for the purpose of obtaining the maximum gain when the 
conditions required for a stable static ~:aracteristic of the elements are satisfied. 


Two-Cycle Elements Which Have Simultaneous Reversal of the Magnetic States of 
the Cores and Are Supplied from a Voltage Source ~ 

The most frequently used two-cycle logic elements with simultaneous reversal of the magnetic states of the 
cores and supply from a voltage source include networks which are based on a repeater cell which is decoupled 
by controlled diodes (Fig. 1) or by resistors (Figs. 2 and 3). 


We shall study the controllability conditions for networks with simultaneous magnetic state reversal for the 
cores in the controlling and controlled elements. 


We shall call that half-cycle of supply voltage which causes conduction of the rectifiers in the controlling or output 


circuit the controlling half-cycle; the half-cycle which cuts off the rectifiers in the controlling or output circuit is 
called the non-working half-cycle. 


* This paper was presented at the All-Moscow Seminar in Contactless Magnetic Elements on October 12, 1960. 
** We shall study logic elements with one core having two windings. 


| 


For a minimal signal at the output,the core of element | in the networks shown in Figs. 1, 2, 3 will be de- 
magnetized, and then during the next controlling half-cycle it will be magnetized. The maximum signal at the 
input of the element I, conversely, will not permit demagnetization of the core; therefore during the next controlling 
half-cycle the core of the element | will be saturated during the entire half-cycle. 


When the core of a controlling element is magnetized, its output voltage is minimal, and the core of the con- 
trolled element II must be demagnetized simultaneously with magnetization of the core of element I due to the 
supply voltage in the controlling circuit. 


If, however, the core of element | is saturated during the controlling half-cycle, then its output voltage must 
be sufficient to prevent a reversal in the magnetic state of the core in the controlled element. 


For an increase in the number of controlled elements,the conditions for simultaneous reversal of the magnetic 
state of the cores do not deteriorate as the number of elements increases, since an increase in the number of elements 
causes a decrease in the minimal output voltage from the 
controlling element. But the conditions for simultaneous 
“non-reversal of the magnetic state" of the cores as the 
number of controlled elements increases deteriorates, since 
a reduction in the maximum output voltage of the con- 
trolling element leads to a situation where (beginning at 
a certain limiting number of elements) the output volt- 
age is no longer sufficient to avoid at least partial de- 
magnetization of the cores of the controlled elements, 
That number of controlled elements for which demag- 
netization of the cores does not yet occur when the core 
of the controlling element is saturated shall be assumed 
to be the maximum possible number. 


In [1] it was demonstrated that in order to obtain 
a stable characteristic for logic elements of the types 
under study, it is necessary that the core of the controlled 
element begin to magnetize later than,and reach complete 
demagnetization sooner than,the magnetization of the 
controlling -element core begins and ends. In this paper 
we do not determine the exact quantitative relationships 
between these quantities, since in practice they may vaty 
over rather wide limits and depend on the quality of the 
rectifiers and the core material. We shall demonstrate 
only that the choice of the optimal relationships in the 
networks shown in Figs. 1, 2, and 3 is directly related 
to obtaining a stable characteristic for the elements; thus 
this choice must be performed while taking stability into 
account. 


We shall make a more detailed study of the network 
cited above under conditions where the cores have a rec- 


tangular hysteresis loop. 
Fig. 3. Network of repeaters with resistor decoupling A Repeater with Decoupling by Means of 
and a common supply voltage. Controlled Rectifiers (Network l) 

For a minimal signal Us in the control circuit of 

the controlling element (Fig. 4a),the core of the con- 

trolling element and the cores of the controlled elements must experience a reversal of their magnetic states during 
a controlling half-cycle. Neglecting the voltage drop across the resistance of the windings and the forward resistance 
of the rectifiers caused by the current that reverses the magnetic state of the core, it is possible to write the following 


expressions: 


U 
a) 


[Je 
" u’ 
Fig. 2. Network of repeaters with resistor decoupling. 
R 
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U’ = kee, (2) 


where U,,;, and U, are the average values of the voltages required for complete reversal of the magnetic state of 
the cores during a controlling half-cycle (Fig. 4b). In order to obtain a stable characteristic (cf. [1]) the voltage U, 
is produced by a source that generates a sinusoidal voltage U' that is cut off at a» at the beginning of the controlling 
half-cycle. We assume that the voltage U' is chosen in such a way that the core is completely demagnetized for 
sat = ® - Gg. Thus, the coefficient kc, which is associated with an increase in the core demagnetization rate in 


ust comparison to its magnetization rate, characterizes the stability of this network. 
etic 
-ments 
n the 
1eous 
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» Since 
on- 
g at 
ts, 
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ned 
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olled Fig. 4. a) Repeater network with decoupling by a diode loaded with 
mplete several repeaters; b) relationship between the voltages in the output 
he circuit of the controlling element and the input circuit of the con- 
aper trolled element during a working half-cycle. 
ships 
y vary 
f the For identical cores, 
ate U. w 
the Te he @) 
ated 
3; thus 
y into Toyo (3") 
_— where Iyjy is the current required for reversing the magnetic state of the working circuit of the element; 1,,, is the 
ape: current required to reverse the magnetic state of the controlling circuit of the element. 
In order for the coupling rectifier V to remain open during the entire controlling half-cycle, the instantaneous 
i values of the currents from the source U" through the resistances Ry, and R, during this half-cycle must be greater 
— than the currents which reverse the magnetic states of the corresponding elements: 
t of 
(4) 
during 
sistance (4") 
ollowing 
The voltage u" is usually produced by an ac source and varies sinusoidally. Therefore expressions (4) and (4') 
(1) cannot, rigorously speaking, be satisfied during the entire half-cycle (Fig. 5), since the magnetizing currents are 


constant quantities for cores with a rectangular hysteresis loop. In order for the equality conditions to be satisfied 


ing 

7 
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in expressions (4) and (4") during the greater portion of the half-cycle, the average value of the current from the 
sinusoidal source must be greater than the average value of the magnetizing current (Fig. 5). 


We shall denote the ratio between the average values of the currents from the source U" and the magnetizing 
currents by the coefficient m; then 


(5) 


In practice it is sufficient for the condition m > 2 to be satisfied. 


It is evident that if relationships (1), (2), (5), and (5") are satisfied in the controlling element and in each 
controlled element, then when the magnetic state of the controlling element is reversed the number of controlled 
elements does not affect their operation. Thus for a minimal signal 
the cores will have their magnetic state reversed irrespective of the 
number of "attached" elements. 


f ut, For a maximum signal at the input of the controlling element, 
its core will not reverse its magnetic state, and the voltage 


(6) 


w 
"4 \ Ve will appear across the rectifier V, where i,,,, is the current passing 
0 x wt through the working winding of the saturated controlling-element 
Fig. 5. Shapes of the curves for the mag- core. 
————eee Under these conditions the cores of the controlled elements 


must not demagnetize in a properly operating network. This con- 

dition will be satisfied if the maximum value of the current in the 

circuit of each controlled element is less than its magnetization 
current: Igy = Ieysts Where 1-11 is the current for reversing the magnetic state on the static hysteresis loop. 


In terms of average values this condition is written as 
(7) 
Ic 


where kg = . 
kq 
Based on this condition, it can be assumed that 


when the resistance of the control circuit is small. 


When condition (7) is satisfied, the current i,,,, is equal to the difference between the currents from the source 
U® and the currents in the control circuits of all the attached elements, i.e., 


(9) 


We shall introduce the substitution U"/Uw; = kggz and shall express U' in terms of U" by using (1) and (2): 
Kkwkgat 


Rey 
=u’ (8) 
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Making use of (5), we obtain 


Fig. 6. a) Repeater network with resistor decoupling and a load consisting 
of several repeaters; b) relationships between the voltages in the output cir- 
cuit of the controlling element and the input circuit. 


If Kiar > 1 (i.e., if the diode V is supplied from a current source), then expression (10) is simplified and be- 
comes 


(10°) 


the substitutions 


(11) 
and 


(2) 


| 
U" ke ke 
5') 
From expressions (8), (6), (9°), (3), (3°), (2) and (1) we find the allowed number of attached elements: | 
d 
nent, 
6) 
wk R ud 
ng | 
u u 
b u 
ts Wk 
NQ 
n- NS u 
the KX y Rt 
Pi 
(8) ky 
ource ic) | akg) 
(9) These expressions can also be represented in the form of power ratios computed in average values if we use 
Finax’ 
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where Py is a quantitity proportional to the power expended in reversing the magnetic state of the core (with an 
accuracy of up to the shape coefficient of the hysteresis loop), and P,,,, is the maximum output power of the con- 
trolling element computed in average values. 


Making use of (10°), (11), and (12), we obtain 


4 
max 
naa 


(mag) He 


We shall demonstrate that the quantity Uwk/tlwky = 4(Pmax /Py) depends on the dimensions, shape and 
material of the core and is independent of the network relationships. We know that 


(14) 
(15) 


(16) 


(a7) 


It is obvious that the quantity defined by formula (17) does not depend on the relationships between the net- 


work parameters; henceforth, it will be assumed constant in the analysis of the networks. 


Returning to the expressions for the number of attached elements (10") or (13), we find the optimal value of 
kw/Ke for which n has its maximum value. For this purpose we formulate the expression 


On 0, 


is) 


The expressions cited above were obtained without considering the heating of the cores. However, beginning 
with certain values of output power from the elements, the heating of the working winding of the cores does not 
permit the use of the elements in the maximum -output mode. 


(13) 
ke 
Owe 
I 
whl “wy 
r= Pp 
; From this we obtain 
(i) 
C/opt a(™— 
which for large kgaz yields (kw/ke)opt = 2. Under these conditions, 
4rl 2 
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(13) 


(18) 


The conditions governing maximum output for the network under study are obtained for Kw/Ke = (kyw/Ke)opt- 
In this case the losses in the core of the controlling element are equal to the output power. If this amount of loss 
is impermissibly great based on the heating conditions, then in order to obtain the maximum possible n,the element 
must have a smaller value of ky/k,. This new value of (kw/Ke)g1; may be determined as follows. 


From (6), (8), (1), (2), and (3) we find the relationship between kw/Ke and Ly»: 
k 
If we denote the current density by A, then 


wk 
ik 
where kf = 9/2, and 1,,),/2 is the average value of the working current over a cycle. 
Making use of (14), (16), and (20) we find that 


(22) 


If we substitute the value kw/Ke = (Kw/Ke) opt into (21), then we can find the value of App, for which the 
element has the maximum power output. 


Thus, for kj,, > 1, kw/Ke = 2,we have 


(23) 
Opt 2pl, 


If this value Aopt = 4,11» Where 44); is the allowable current density, then we can choose kw/Ke = (Kw/Keopt 
and the number of “attached” elements will be equal to Dopt in accordance with (18). 


However, if Appt >Aaii , then the quantity kw/kc must be determined from the allowable current density 
given by (21) and the number of attached elements can be determined from formulas (10), (10°), or (13). 


In the case where ky/k, is not chosen in accordance with expression (21),we will have a situation where for 
a specified current density the network will not operate in a limiting mode,and the number of attached elements 
will be less than the number found for formulas (10), (10°), or (13). 


A Repeater with Resistor Decoupling (Network II) 


~ The network for such a repeater with a load that also consists of repeaters is shown in Fig. 6a. 


Figure 6b shows the distribution of the voltages among the network elements when the magnetic state of the 
cores is reversed. 


For a sinusoidal voltage this network (due to the presence of the resistor R,) can operate without any special 
cutoff in the supply voltage U' for the control circuit, since before the cores of the controlled elements begin to 
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demagnetize (i.e., for wt <a), the flux in the core of the controlling element increases by the amount A, 
where AD = Wye ; after complete demagnetization of the controlled-element cores (i.e, 
for Wt > Osa, ) the flux in the core of the controlling element will continue to increase up to +s, and the relation- 
ship Ao = (u—u’)dot. will again be valid. The ratio / (@;— in essence determines the 


stability of the element characteristic, since the greater its ratio,the greater the slope of the characteristic Uour/ tip 
for the element. 


In the network under study this ratio is determined by the area (Fig. 6b) bounded by the curves u and u’ from 
0 to a» and from a ,,, to #, and by the area bounded by the curves u and i,R which is proportional to Uy . It 
depends both on the value kw = Uy /U, and on the choice of the coefficient 


(24) 


which determines the angles a and @sq;. The stability increases with an increase in k,, and k,. In the network 
under study it is not expedient to choose kw less than 2; however, a growth in k, also causes a growth of the losses 
in the resistor R and increases the power drawn by the logic elements. 


Just as in the preceding network, for a minimum signal U, in the control circuit of the controlling element,the 
cores of all the elements in the network shown in Fig. 6a must simultaneously be subjected to a reversal of the 
magnetic state. Neglecting the voltage drop caused a) by the currents which flow through the winding resistance 
and rectifiers and cause the reversal of magnetic state and b) by the current I w_y, it is possible to write the following 
relationships for the controlling half-cycle: 


5 
u’ = ug + ig RK. (26) 


As in the case cited above [formulas (3) and (3")], 


If the parameters of all the elements are chosen in such a way that Eq. (25) and (26) are satisfied, then, when 
the core of the controlling elements is subjected to a reversal of magnetic state,all the remaining elements will 
be subjected to a reversal of magnetic state irrespective of their number; this follows since it is a fact that, just as 
in the preceding network, the elements do not affect one another when their magnetic states are reversed. 


For a maximum control signal the core of the controlling element remains saturated during the entire con- 
trolling half-cycle. 


For proper operation of the network,the cores of all attached (controlled) elements must also remain in a 
saturated state. The condition for this is that the maximum current flowing through the control winding of each 
controlled element must remain less than the demagnetization current Icyst on the static hysteresis loop, i.e., (cf. 
formula (7)] 


c,sat d 


2 


Neglecting the voltage drop across the control windings of the cores and across the rectifiers, we obtain the 
following equations for the saturated state of the cores (the subscript “sat " denotes the saturated state of the cores): 


wk,sat 


U 
“wk_ }. 
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tion - where the current i, sar is the current in the working winding of the saturated core of the controlling element. 

Making the transition to average values and substituting expression (7) into (27) and (28), we obtain the limit- 
ip ing value for the number of controlled elements: 
from (29) 

(24) From (25) we find the following result for the average values of the currents and voltages: 
tk 
a where k, = ‘on > 1 (cf. Fig. 6b); from (25), (26), and (3) we obtain 
U—U'=Ufi—). 
1ce 
From (29) we obtain 
(25) 
(26) 
while taking (30), (31), (24), and (3) into account, or if we use (11) and (12) we obtain 
ky — A 

ynen An analysis of these expressions demonstrates that n will increase as k, increases fork, > 1. Fork, <1, 
ll expression (33) does not have any physical meaning, since in such a case the core cannot be fully demagnetized. 
a It is evident that the maximum value of n is achieved for 
kw opt = 2. (34) 
: For a sufficiently large kp and ky = kw, opt »expressions (32) and (33) become simplified and are written as 
(ef. 
the 
cores): 


If we again introduce limitations for the heating, then k,, must be chosen in accordance with the following 
reasoning. 


) 
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From expressions (27), (28), and (31) we obtain the relationship between ky and L,,,, which is analogous to the 
relationship for the preceding network: 


= 
Performing the same transformation, we again obtain 


A x 
opt pl; 10° 


If Aopr< Man where 4,1) is the allowable current density, then it is possible to choose k,, = 2; the number 
of attached elements will be equal to nope in accordance with (32). If App > 4,1) . then the quantity k,, must be 
determined from the allowable current density in accordance with (21),and the number of attached elements can be 
determined from (32). 


However, if k,, from (21) is found to be less than two, then (as was shown at the beginning of this section) the 
stability of the network deteriorates. Therefore it is expedient to choose k,, > 2. This leads to a situation where 
the number of attached elements must be less than the limiting number which we derive from the condition of 
maximum output. Thus, since the number of attached elements will be less than the limiting number, the voltage 
across the coupling resistor will exceed u' and the rectifier V, in the control circuit will be cut off. This means 
that coefficient kg in expression (32) is equal to ; therefore, the limiting value of n decreases to the value deter- 
mined from (32) for kg = ©,and the allowable value of n will be less than the limiting value by a factor of A/Agy, 
where 4 is the value of the current density obtained from (21) for the chosen value of k,,. Thus 


They ke 


Substituting A into this expression from (21), we obtain 


Making use of (14), (16), and (22), we find 


kok, 


where Ama, tepresents the maximum allowable ampere turns on the working winding of the core on the basis of 
heating considerations, and A,, is the number of coercive ampere turns of this winding. 


In expression (36) the value of k, must be chosen to be the same as ky while considering stability require- 
ments. 


where 
C pi, 
and for ky = Ky opt + 
A 
U. (le 
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the 


of 


The choice of k- also determines the choice of the relationship between the supply voltage and the voltage 


Uwk across the working winding which causes a reversal of magnetic state. We shall use the substitution k,,, = 
= U/Uwk- On the basis of expression (31), 


U Ut Owe 


+ 


A Repeater with Resistor Decoupling and One Supply Source (Network III) 
The network for such a repeater with a load that also consists of repeaters is shown in Fig. 7a. 


Figure 7b shows the distribution of voltages among the network elements when the magnetic state of the cores 
is reversed. For a sinusoidal supply voltage this network may be stable due to the presence of the resistances R and 
Re (this is analogous to the situation prevailing for the preceding network). 


Fig. 7. a) Repeater network with a common supply voltage and resistor 
decoupling which is loaded with several repeaters; b) relationship be- 
tween the voltages in the output circuit of the controlling element and 
the input circuit of the controlled elements during the working half- 

cycle. 


The conditions governing a reversal in the magnetic state of the cores of the controlled elements when a core 
of the controlling element reverses its magnetic state are written as 


(37) 


u— U Witt be R, 
(38) 


The satisfaction of these conditions, as in the preceding networks, depends not on the number of controlled 
e’ements, but only on the proper relationships between U, R and R¢ for specified cores and a definite relationship be - 
tween the number of turns on their windings. 


nber 
t be a 
= a, sar wt 
on 
“wk. 
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The condition governing the possibility of controlling n logic elements when the core of the controlling elemey 


reverses its magnetic state is also based on the relationship 


Jem Seu an 
or, in terms of average values, 
4 
and on the relationships between the parameters of the working and controlling circuits: 
rl = Re 


[y= eset 
from (39) and 
RoR 
~ 
from (40) and (41). 


From (37) we find the following result in average values in accordance with Fig. 7b: 


ka 


where = Ipy/Ic- 
From (38) we obtain (in average values) 


In our subsequent analysis we also introduce the coefficient 
c _ Sad 


which together with k,, determines the stability of the characteristics for the elements. 


(39) 
(40) 


(41) 


(42) 


(43) 


(44) 


(45) 


The choice of k,, determines the angles 0 and og,, and the termination of the process involving the reversal 
of magnetic state for the cores of the controlled elements. As k, and k, increase, the stability also increases; how- 


ever, the losses in the resistances R and R,, increase (i.e., the power drawn from the supply sources increases). 


Substituting the relationships (7), (3), (43), (44), and (45) into (42), we find 


U ky ky 
wig 

2h, ke 


(46) 


. Making the transition to average values, we obtain 
cH 
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The value of n increases as k, increases and has a maximum as a function of k,, for 


Ky,opt 


For large ke we have kw opt = 2. 
We shall now proceed to a determination of the heating limitations. From (41) and (44) we obtain 


Making use of (20), (14), (16), and (48), we find 


2k 


where C is determined from (22). 


From this we can find Aop, if we make the substitution kw = kw, opt - In the case where > Say, the choice 
of ky may be made in accordance with (49) relative to 4,)), and the number n can be determ in accordance 
with (46). However, based on stability considerations, the value of kw cannot be less than two. But since under these 
conditions the current density for the limiting (based on the conditions governing non-reversal of the magnetic state) 
number of elements becomes greater than the allowable value, itis necessary to retreat from the limiting number 
of elements and connect a number such that the current density does not exceed the allowable value. This number 
can be determined from (46) for a new value of kq which is now no longer equal to Icy / Icyst but is determined 
from the following expression on the basis of (49): 


(50) 


We shall use numerical examples to compare the networks studied above for characteristic operating modes. 


In order to compare the networks for operation at a frequency of 50 cps we perform a computation for cores 
having the following data: material 50 NP, Hg = 0.2 amp/cm, Bry = 15000 gauss, / , = 8.3 cm, 7, = 5 cm, 8 = 
= 0,18 cm’, s = 1.8 cm’, Kgp = 0.3. 


In accordance with (17), 


C = 0.0021 cm/amp. 
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In accordance with (22), : 


We shall compute the number of elements which can be controlled by one element (i.e., the gain of an element) 
for identical values of k,,, and values of k,, close to the optimum values for each network. 


a. Fornetwork I, (kw/ke)ope=2- For stable operation it is sufficient that k, = 1.5, and therefore we choose k, = 
= 1.5; kw = 3, m = 2, kq = 2, kg, = 2. Then according to expression (10), 


nw = 15. 


b. For network II we choose ky = 3, kq = 2; then k, = 4. For these values k, 4/3 (this relationship is 
determined by graphical plotting in a manner analogous to that shown in Fig. 6b). In accordance with expression 
(33), 


n 65. 


We shall verify the conditions governing the heating in such networks. For both networks the current density 
for k,, = 3 is determined from expression (21): 


A = 315 amp/cm*, 


For cores with the selected dimensions, A,; = 400 amp/cm*. Thus, networks I and III can operate for an 
optimal value of ky, i.e., for maximum output. 


c. For the network Ill we have k, = 4 and ky op, = 1.8 from (47) for kg = 2. However, based on the stability 
conditions we should not choose ky to be less than 2; we choose k,, = 2, and k, = 4/3. Then from (46) we have 


For this network the current density is 


100 amp/cm?, 


opt 


for the chosen value k,, = 2, i.e., Aopt<4ay; and the network can operate for the chosen value of kw. 


The examples cited above clearly show the advantages of network II over networks I and III, which were com- 
puted for optimal parameters with respect to the maximum gain. The gains for network II are more than four times 
the gains for networks I and III. 


Here the power loss in the coupling resistors in network II is less than the power loss in the resistors R, in net- 
work I and in the resistors R and R, in network III. 


In order to compare the networks for operation at a frequency of 400 cps,we shall perform the computation 
for cores with the following data: H, = 0.4 amp/cm, Bm = 14500 gauss, 7, = 6 cm, = 4m, s,= 0.1 cm’, % = 
= 1.2 cm’, kgp = 0.3. 


In accordance with (17), 


4P 
max _ 500. 
p 


C == 3.85-10¢. 


First we shall compute the number of elements which one element can control for identical values of Kgat 
and values of kw which are optimal from the viewpoint of maximum power output. 


| In accordance with (22), 
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a. For network I(kw/Ke)opt = 2, where from the stability conditions ke = 1.5; therefore, kw = 3, k, = 3, 
Kat, = 2: 
Making use of (10), we find 


" opt< 36. 
Under these conditions we obtain 
Agpt = 1350 amp/cm* 


from (21). 
Since Appt > Aall. we find that the allowable value (kw/Ke)aii for a limiting mode of operation is 


“all 


ky att = 1.76. 


for k,, = 1.5. 
Under these conditions we have 


n = 16. 


Thus, when the heating effect is taken into account,one element in this network can control only 16 elements. 


b. For network II we have kw opr = 2, kq = 3; thus, according to the stability conditions,we have ke = 4, 
and kgat = 2.5. When heating is not taken into account, 


n < 142. 


Under these conditions App = 1350 amp/cm”. We shall find the allowable ky for which A,y) = 400 amp/cm’. 
Based on (21), 


i 
0.846, k,, = 1,2. 


However, for the value k,, = 1,2, it is not possible to obtain a stable characteristic for the elements. From 
the stability conditions for this network, it is expedient to choose ky no less than two, In accordance with (36), we 
obtain 
n= 22, 
c. For network III we choose ky = 3, as before. From the stability conditions we choose k, = 4, and then 
kw, opt = 1.6. Therefore, we have 
n= 22, 


Under these conditions 


oh = 380 amp/cm®, 


Le., Aopt < 4,), and the network can still operate in a limiting mode. 
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SUMMARY 
From a comparison of the derived formulas, whose use has been illustrated by means of several examples, we cap 
draw the following conclusions. 


1. The optimal choice of parameters is determined in the same way as it is in linear networks or by the con- 
dition governing the matching of the input and output resistances of the elements; this condition leads to a mode that 
is close to maximum power output from the controlling element or to a limitation of the output power by the heating 
effect in the controlling element. 


In the first case (for maximum power output) the ratio k,, between the number of turns on the working and 
controlling windings is confined within the limits 2 to 3; in networks II and III this makes it possible to achieve a 
stable characteristic for a ratio k, (the ratio between the control circuit supply voltage and the voltage which produces 
a reversal in the magnetic state of the controlling winding) that is confined within the limits 4 to 8. In network 
I stability is achieved by choosing k,, in the range 1.2 to 1.5. 


In the second case, for a limitation based on the heating effect,the value of ky in network I is chosen in ac- 
cordance with requirements governing the operation of a network with the limiting number of elements, and the 
stability is assured by choosing k, within the range 1.2 to 1.5; in networks II and III the value of k,, must be chosen 
from the condition governing the achievement of a stable characteristic, and thus it can be no less than two. Under 
these conditions, the limiting operating mode for networks II and III cannot be achieved, and the number of con- 
trolled elements must be determined from the allowable current density on the basis of expressions (36) and (51), 
which take into account the deviation from the limiting operating mode. 


2, Of the networks studied above,the most useful are networks II and Ill,which yield the maximum values of 
n. Network II may be preferable to networks I and II for operation with heating limitations. 


3. The derived formulas make it possible not only to determine the gain n for specified core parameters, but 
to determine the core dimensions for a specified n, Here from expressions (10), or (32) and (34), or (36) and (51), it 
is possible to determine the value of Pryay/Py or Amax/ Ay; then from formulas (17) or (36) we choose the corre- 
sponding dimensions and other parameters for the element. 


Here in the cases of networks II and III it is necessary to find Pryay/Py, or Amax/ A, from the formulas which 
determine n on the basis of the conditions of maximum power output and from the formulas which take heating into 
account; then we choose the largest value. 


4. Since the chief purpose of this paper was a comparison of various networks on the basis of the maximum gail 
for stablé operation, we did not take into account certain factors which would have complicated the analysis substan- 
tially. These factors include fluctuations in the supply voltage, a spread in the parameters of the cores and resistors, 
the effect of the actual rectifier characteristics, and variation of the demagnetization current with the demagnetiza- 
tion rate. Therefore,when the formulas cited above are used for computation purposes, we should treat the values of 
n derived above as the upper limit. 
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The paper describes a transmitting unit for a pulse-frequency telemetering system; the unit is 
based on semiconductor elements. The unit employs a new type of dc to ac converter — 4 
magnetic modulator with transverse excitation. The range of variation for the output frequency 
- from the unit is 5 to 15 cps. Due to the high sensitivity of the unit (the input signal power is 
10-* w) the unit can be used for transmitting readings from dc transducers which have a low 
output power. 


The pulse -frequency telemetering unit "ChTI-1" designed in the Institute of Automation and Remote Control 
of the USSR Academy of Sciences is designed for measuring low dc voltages. The “ChTI-1" transmitter produces 
rectangular pulses with a duty ratio of two and a repetition frequency proportional to the measured voltage. The 
pulses from the output of the transmitter can be transmitted directly over communication line wires or over fre - 
quency division channels. The increased noise immunity requirements of modern telemetering transmission using 
frequency channels gives rise to the problem of narrowing the frequency band in pulse-frequency telemetering. For 
this reason the range of variation for the pulse repetition frequency is chosen between 5 and 15 cps. A further lower- 
ing of the upper frequency by modern engineering means, including the use of transistors, would lead to an im- 
permissible reduction of the operating stability of the unit. 


For telemetering transmission over frequency division channels,the pulses from the transmitter modulate the 
generator producing the sub-carrier frequency. The chosen frequency range (upper frequency 15 cps) makes it possible 
to achieve extremely narrow-band transmission (this therefore results in a high noise immunity). The “ChTI-1" 
unit is matched to a simple narrow-band frequency division unit [1] which was developed in the Institute of Automa- 
tion and Remote Control of the USSRAcademy of Sciences; this makes it possible to transmit telemetering data over 
a band 40 cps wide in any section of the total range. 


We can use any one of numerous frequency metering networks as the receiver in the pulse ~frequency tele- 
metering unit “"ChTI-1." For example, a number of modern designs of transistorized capacitor frequency meters 
are available which have a stabilized supply voltage [2]. Because of this the network for the receiver is not cited 
in this paper. 
The Block Diagram for the Transmitter 

The "ChTI-1" transmitter operates according to the compensation method and is essentially a rectangular 
pulse generator with a pulse duty ratio equal to two; it is controlled by the measured voltage and is encompassed with 
a large negative feedback. The transfer coefficient for the closed system is approximately equal to 1/6 for large 
values of k6 .° Therefore the stability of the transfer coefficient will be determined basically by the feedback cir- 
cuit which is much easier to make highly stable than is the direct channel. The block diagram for the unit (Fig. 1) 
consists of a low-pass filter 1, a dc amplifier 2 containing a modulator -converter that converts a dc signal to an ac 
signal, an ac amplifier and a rectifier, a pulse generator 3, and a block 4 which performs the reverse conversion 
of a frequency to a voltage. 


* Here k is the transfer coefficient for the direct channel; 6 is the transfer coefficient for the feedback circuit. 
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The position of the filter in the direct channel permits a substantial simplification of the amplifier circuit 
(the phase -detecting stage* can be eliminated) due to the fact that there is no inertial section in the feedback 
circuit. This is possible, of course, in the case where voltages which do 
not change sign are measured. When there is no phase -detecting stage it 
is necessary to choose the filter parameters correctly in order to achieve 
stable operation of the network. We shall compare the expressions for the 
transfer functions of the device in two cases: 


1) The filter is located in the direct channel: 


4 


2) the filter is located in the feedback circuit: 


Here ; (p) = . Sa ear is the transfer function for a two-section RC filter; ke is the amplifier gain; 
kg is the voltage to frequency coefficient. The transient functions which determine the nature of the transient re- 
sponses are equal to the following expressions in the first and second cases, respectively: 


kk 
(1— 


where p; and ‘p, are the roots of the equation ap’ + bp + c + kpkg B = 0. 


4n 


Fig. 2, Response of the unit with the filter Fig. 3. Response of the unit with the filter 


in the direct channel to a rectangular pulse. in the feedback circuit to a rectangular 
a) Aperiodic mode; b) oscillatory mode. pulse. 


Making use of the relationships derived above, we plot the curves for the variation of the frequency at the 
output of the unit when it is subjected to a rectangular input signal (Figs. 2 and 3). The difference between the time 
characteristic for the output frequency and the time characteristic for the voltage at the modulator input involves 
only a constant coefficient. Therefore, if we choose the system parameters in such a way that the response to a 
rectangular input signal has the shape of curve a in Fig. 2, the voltage at the modulator input will always be of the 
same sign. Thus when the filter is located in the direct channel and the transient responses are aperiodic,the system 


can operate without a phase -detecting stage. 

* The phase-detecting stage is usually employed in networks which amplify a direct current by converting the direct 
current to alternating current, since the phase of the ac current at the output of the converter depends on the polarity 
of the direct current at its input. 


| 


In the case of a two-section low-pass RC filter, it is easy to derive relationships from which it follows that 
stable operation of the network requires the time constant of the first section to be much greater than the time 
constant of the second section. In view of the fact that the filter must achieve good filtering of very low frequencies 
(from 5 cps up), it has appreciable inertia. However the inertia of the filter has little effect on the speed of re- 
sponse of the system, since the “hard” negative feedback reduces this inertia by a factor of (1 + kB). 


One of the advantages of the block diagram we have chosen is the fact that due to compensation of the input 


signal by the feedback voltage in the measured -voltage circuit® it is possible to increase appreciably the input 
resistance of the unit. 


Circuit Diagram 


The circuit diagram for the "ChTI-1" transmitter is shown in Fig. 4, The device which converts the dc 
imbalance voltage at the filter output into a dc voltage is a new type of magnetic modulator with transverse excita- 


tion [3]. The magnetic fields of the excitation current and the measured current are spatially separated in the core 
of the modulator M (the windings W, and W2, respectively); therefore, it is possible to achieve a sharp increase in 
null stability. The transfer coefficient of the modulator is equal to approximately 0.8-0.9 and does not depend on 


wide variations in the voltage and frequency of the supply source and in the temperature of the ambient medium. The 
modulator is practically inertialess. 


The ac amplifier consists of two paired stages of transistors "P13B" and "P14" (these are the transistors T,, T, 


the and Ts, T, in the circuit shown in Fig. 4). The first transistor in such a stage uses a common-collector connection. 
she time The input impedance of the stage is sufficiently high, and therefore the stages can be connected together by means 
nivel of RC coupling. Such interstage coupling permits us to achieve temperature stabilization easily. 

oa The rectified and filtered output voltage from the amplifier controls the frequency of a blocking oscillator 

re of the (the transistor T;). The use of a blocking oscillator as a voltage-to-frequency converter is based first on its high 

system 

e direct * In contrast, for example, to magnetic flux compensation in the cores of null-elements in pulse-frequency tele- 

‘polarity Metering systems of the type “ChIS-D-1" and “ChIS-D-2" developed by the Central Laboratory and Experimental 


Workshop of Electrical and Measuring Instruments [2, 4]. 


sensitivity with respect to the controlled voltage (in order to change the frequency from minimum to maximum ye 
require a voltage of 2 v across a 50 kohm resistor), and second on the possibility of achieving temperature stabiliza- 
tion in simple fashion (temperature stabilization is achieved with a silicon diode and a shunting resistor in the base 
circuit; stabilization with respect to the supply voltage is achieved by using a silicon stabilizing diode). 


A trigger consisting of transistors of the types “P2B? “P25; or "P26" is used to shape the output pulses; the 
primary windings of the differentiating transformers are connected in the collector circuits of these transistors. The 
voltage from the secondary windings of the transformers is rectified by diodes and becomes a pulsating voltage with 
a dc component that is proportional to the repetition frequency of the pulses; this voltage is then applied to the inpy 
of the unit in order to compensate the input voltage. The voltage which acts at the input of the unit as a result of 
the conversion of the lowest frequency (5cps) is compensated by a special tributary voltage that is stabilized by a 
silicon stabilizing diode. 


The *ChTI-1" unit was subjected to laboratory tests whose basic results are reflected in the engineering charac- 
teristics of the unit. 


The Engineering Characteristics of the "ChTI-1" 
The entire scale of the telemetering unit (i.e., the range of variation of the transducer voltage which corre- 
sponds to the full range of variation for the frequency) is 25 millivolts. The range of variation for the pulse repeti- 
tion frequency is 5 to 15 cps. The input resistance of the transmitter is of the order of 50 kohm. The output signal 
consists of dc rectangular pulses that produce a voltage of 20 v across a resistance of 800 ohm. The gain of the 
open-loop system is 25. The basic error caused by the nonlinearity of the feedback circuit does not exceed # 1%, 
The additional errors which appear when the temperature of the ambient medium varies from 0 to +40°C and the 
supply voltage varies by + 15% do not exceed +1%. The time required for the frequency to reach a steady-state 
value is 0.5 sec. The unit has the dimensions 255 x 215 x 140 mm. 


The Distinguishing Features and the Field of Application of the ""ChTI-1" 
The pulse-frequency telemetering unit “ChTI-1" has the following basic distinguishing features in comparison 
to analogous units such as the "ChIS-D-1” or the "ChIS-D-2” [2, 4). 


1. The use of a new type of magnetic modulator instead of modulators designed according to the magnetic am- 
plifier principle has made it possible to achieve a considerable simplification of the circuit. 


2. Due to the rationally chosen block diagram,the phase-detecting stage is eliminated. 


3. Compensation in the measured -voltage circuit makes it possible to provide a high-ohm input for the tans 
mitter; because of this the sensitivity is increased (the power required to control the transmitter is 10~* w), and the 
range of application of the device is expanded. 


The telemetering unit "ChTI-1" can be used for transmitting readings produced by a wide class of dc transducet 
that have a low output power. These include various transducers in which thermal elements are used: temperature 
transducers, electric power transducers, etc. Moreover, a small output voltage also characterizes individual tension 
transducers and other bridge networks such as gas analyzers. An example where a gas analyzer is used as a transduce 
in a telemetering system is a system for telemetering the methane concentration in mines. 


LITERATURE CITED 

1. Nz. V. Pozin, On the Transmission of Telemetered Values over Telephone Channels. Transactions of the Sessio 
of the USSR Academy of Sciences on the Scientific Problems of Production Automation, October 5-20, 1956 
{in Russian] (USSR Academy of Sciences Press, 1957). 

2. A.M. Pshenichnikov, A Static Pulse-Frequency Telemetering System Based on Magnetic and Semiconductor 
Elements. Scientific-Engineering Conference on the Problems of Design and Production of Remote Control 
Equipment and Communications Channels for Remote Control [in Russian] (Leningrad, 1958). 

3. _‘F. 1. Kerbnikov and M. A. Rozenblat, "Magnetic modulators with transverse excitation,” Avtomatika i 
Telemekhanika, 19, 9 (1958). 

4. A.M. Pshenichnikov, “A static transmitting unit for a pulse-frequency telemetering system,” Avtomatika i 

Telemekhanika, 18, 5 (1957). 


$; the 


Ze with 


sult of 


ic am- 
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The paper cites tables and formulas which facilitate the computation of the dispersion (as a func- 
tion of time) and the correlation function of the output variable for a linear stationary dynamic © 
system when its input is subjected to a random process with a known correlation function. Cer- 

tain possible methods are cited for determining the limiting accuracy characteristics for the oper - 
ation of the system: the steady-state correlation function and the steady-state dispersion. 


1. Determination of the Correlation Function of the Output Variable 


The determination of the correlation function of the output variable of a linear stationary dynamic system is 
usually performed using the formula [1] 


K = \ett, 5) 8 (ti, 81) Kx (s, 81) dsds,, (1) 


where Ky (t, ty), Ky (t, ty) are the correlation functions for the input and output variables,respectively, and g(t, 8) 
is a weighting function. 


Formula (1) is applicable when the transfer function for the system has an arbitrary form; it does not assume 
that the random input signal of the system is stationary and makes it possible to take into account the time for which 
the system operates. The shortcomings of (1) are the necessity for the preliminary determination of a weighting func- 
tion and the complexity which is involved in performing the double integration. 


This paper presents a different method for determining the correlation functions of the output variables of a 
linear stationary dynamic system. This method does not require the determination of a weighting function and 
eliminates the necessity of integrating. The computations are reduced to performing a series of simple operations. 


The basis of the method is the relationship [2, 3] 


(u, us) = W (u) W ®,, (u, uy), (2) 


in which #,,(u, Uy), ®jy(u, uy) respectively denote the results of two-dimensional Laplace transformations of the 
correlation functions for the signals at the input and output of the system; W(p) is the transfer function for the 
system. 


First we shall assume that a stationary random signal with the correlation function K,(r) = K,( | t-t|) is 
applied to the system input. Under these conditions 


(u) + F, 
u + uy 


(3) 


®,, (4, %) = 


where F,(p) is a unidimensional transform of the function K,(r), 
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p+a 


+e 


—Bsi 


co 
(p) = \ at. 
(4) 


Table 1 shows the transform for the functions K,(7T) that are frequently used in solving practical problems. 
It is evident that the function F,(p) can in many cases be represented by the proper fraction 


(5) 


where M(p), N(p) are polynomials of degree m and n, respectively, relative to the variable p(m < n). We shall 
assume that the function W(p) can also be represented by the proper fraction 


W (p) = (6) 


The quantities q and r-shall be used to denote the order of the polynomials Q(p) and R(p). 
The position of the poles of the function W(p) [the zeros of the function R(p)] will be characterized by a syste 
of quantities A; for 1 <i <r; the position of the poles for the function F(p) [the zeros of the function N(p)] will be 
characterized by a system of quantities A; forr+1<i=r+n. During the first stage of our investigation we 
shall assume that all poles of the function W(p) and F,(p) are simple and that the position of the poles for the fune- 
tion Fx(p) does not coincide with the position of any of the poles for the function W(p). 


We shall represent the fanction W(p) and the product of the functions W(p)F;(p) in the form of a sum of simp 
fractions: 


r-+n 
W (p) F = W = 


M (hy) 


TABLE 1 
Transform F;(p) Original Ky{p) 
+ 
i 
where 
832 


(4) 


®,, (u, wy) = 


F (dy hj) = 


Making use of this representation for transforming the formula 


W (u) W [F,, (u) + 


which is a corollary of formulas (2) and (3), we obtain 


u+u, 


r+n r 


i=1 j=1 
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W(p)= >) 


f (Ais 45) = 


WF. 
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p+d,. 


cu+d 


cu, +d 


w +au+b’ uy + au +b 


Formula (10) in the region of the originals corresponds to the formula 


T = min (t, t). 


r+n 


i=th+1 


in various allowed combinations. 


(4 eit 5) 


If all the quantities \; are real then the result is final, However, if among the quantities \; there are complex 
quantities then additional transformations will be needed to eliminate the imaginary values. 


We can avoid performing additional transformations if we formulate the corresponding table of originals and 
transforms for the case of a two-dimensional transform. Starting from the formulas 


4 
—A 


which are equivalent to formulas (7),and for b,>—-(1 <-s < &) representing the expansion of functions W(p) and 
W(p)F,(p) in simple rational fractions, we can conclude that in the transform column of this table it is necessary 
to write not only the functions F(Aj, 4j) but also the functions which have a structure analogous to the functions 


F(A, Aj) and contain terms of the form 


The originals corresponding to these transforms can be found using formulas (12), (13) if we first represent each 
newly introduced tabular transform by the sum of terms of the form F(A, 4). 


Table 2 is a correspondence table for the originals and transforms which assure the computation of the proba- 


bility characteristics for the variable at the output of a stationary dynamic system when there are no multiple roots 
for the function W(p) and the product W(p)F,(p). 
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The method of tables for determining the correlation function of the output variable assumes that the uni- 
dimensional transform of the correlation function for the input signal is found; then formula (9) is used to find the 
two-dimensional transform of the correlation function for the output variable, the transform for the correlation func- 
tion of the output variable is converted to a sum of tabular transforms, and finally the transition is made into the 
region of originals. 


In accordance with Eq. (12) each pole of the transfer function in the expression for the correlation function 
of the output variable corresponds to n components. Comparing the values of the individual components, it is possible 
to establish which of them are the most important (usually it is those corresponding to the smallest values of the 
parameters ;), and thus it is possible to clarify the possibility of replacing the equation that describes the behavior 
of the specified dynamic system by another simpler equation. 


It is evident that this same method for finding the correlation function of the output variable can be extended 
to apply to the case where a nonstationary random input is applied to the system under conditions where the trans- 
form for the correlation function of the signal is represented in the form of a linear combination of tabular transforms, 
In particular, the problem of investigating the process involved in the passage of a signal through serially connected 
sections of a dynamic system whose input is subjected to a stationary random process reduces to this case. 


If the function W(p) or the product W(p)Fx(p) have multiple poles, then it is impossible to determine the 
correlation function for the output variable solely by means of Table 2. 


2. Determining the Dependence of the Dispersion of the Output Variable on Time 
In order to find the dispersion of the output variable of the system, it is necessary to set t = t, in Eq. (12)-(14). 
Under these conditions r = 0, T = t, and therefore 


r+n yr. 


Dy (t) = K,(t,)=2>5 5 


it j--1 


Formula (15) determines the dependence of the dispersion for the output variable of the system on time for a 
known position of the poles of the functions W(p) and Fx(p). However, very frequently the position of the poles of 
the functions W(p) and F,(p) is unknown, and therefore it is expedient to transform formula (15) to the expression 


t 
D, (t) = x, (u) x, (u) du, 


r+n r 
>) Aye’, = (17) 


Comparing Eq. (17) with relationships (7), we arrive at the conclusion that the transforms W(p)F,(p) and W(p) 
in the region of originals corresponds to the functions x,(t) and x(t). This proves the possibility of using the following 
method for computing the dependence of the output-variable dispersion on time: 


1) Using the tables for correspondence of the formulas when the Laplace transform is used, we find a unidimen- 
sional transform for the output-variable correlation functions F,(p). 


2) Using the same tables, we determine the originals x;(t), x(t) corresponding to the transforms W(p)F (Pp) 
and W(p). 


3) Using formula (16), we compute the dispersion Dy(t). 
A corollary of formula (16) is the formula 


co 
D, = 2( 21 (w) 24 (w) du, 
0 


which expresses the steady -state value of the output-variable dispersion. 


We find the unidimensional transform of the function Dy(t). In accordance with the basic properties of a 
Laplace transform [4], this transform is determined by the formula 


joo 
Nip) = (2) W (p— #)ds, as) 
1P e—Joo 


where c is the abscissa of the absolute convergence for the products x4(t) x(t). 
When the input of a definitely stable dynamic system is subjected to a stationary random signal, c < 0; there- 


fore formula (19) establishes the value of the function /\(p) at the point p= 0. Applying the final-value theorem 
[4], we obtain 


nded 
D, = li sale 
= lim pit (Pp) 2) dz. on 
>cted 
The validity of formulas(19), (20) is proved for the condition that the poles of the transfer function W(p) and 
the products W(p)Fx(p) have a multiplicity equal to unity. The cases where the functions W(p) and W(p)F,(p) have 
ov multiple poles can be treated as limiting cases with respect to the case treated previously. For the transition in the 
Tid), limit A; Aj, the function W(p) and W(p)Fx(p) do not lose their continuity relative to the variable p. Therefore 
it can be assumed that formulas (19), (20), as well as their corollaries, are valid for an arbitrary multiplicity of the 
poles of the functions W(p) and W(p)Fx(p). This fact proves the possibility of using the method cited above for 
determining the function Dy(t) and the value of Dy when the form of the transfer function W(p) is arbitrary and the 
(15) transform of the correlation function F,(p) is practically arbitrary. 
If the position of the poles of the functions W(p) and F,(p) is known and the multiplicity of each pole is also 
known, then it is possible to use the theorem of residues [5] in computing the steady -state dispersion Dy. 
“a In the general case we shall assume that for 1 = i = k the quantities ); characterize the position of the poles 
sol of the function W(p), and the quantity mj determines the multiplicity of the poles Aj. | Performing summation over 
the poles of the function W(-p) that are located symmetrically with respect to the poles of the function W(p) relative to 
(16) to the point p = 0, we can use the theorem of residues to obtain 
Dy = 2 [W (2) W (—2) F, (— 
=1 z 
(17) 
r 
Q(—A,) M (— Q(A,) 
D,=2 (22) 
| W(p) RO) 
ollowing 
is a particular case of formula (21) and determines the steady -state dispersion of the output variables of the system 
anit when all poles of the transfer function W(p) have a multiplicity equal to unity (for k = r). 
3. Determining the Steady State Correlation Function of the Output Variable 
(p) From table 2 it follows that when a stationary random signal is applied to the input of a linear stationary 
| dynamic system the output signal theoretically cannot be assumed stationary. If the transfer function of the system 
has a pole at the point p = 0, or if it has poles at the points p = 4 A (A is an arbitrary number), then the quantity T 
is included in one of the components of the function Ky(t, t,) as a multiplier; therefore the output signal does not 
become stationary from a practical point of view even for the case where the preliminary operating time of the 
system is large (a large value of T). The signal cannot be considered stationary for a large value of T at the out- 
as) put of an unstable dynamic system, since in this case an unlimited increase in certain of the components of the 


function K (t,t) occurs due to the positive nature of the real part of at least one of the parameters A; for T+ ©, 
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If the dynamic system is definitely stable, then for a sufficieatly large preliminary operating time of the system 
(large T) the nonstationary parts of the components of the function Ky(t.ty) become small in absolute magnitude; 
therefore we“ assume that the output signal is practically stationary. 


The c. -lation function and dispersion of the output variable of a definitely stable dynamic system approach 
their limits for an unbounded increase in the preliminary operating time for the system (when T + «): they approach 
the values of the steady-state correlation function and the steady -state dispersion [1]. In many cases involving the 
solution of practical problems it is expedient to find precisely these characteristics of the operational accuracy of 
the system. 

At the present time we usually make use of the following formulas [1] in determining the steady -state charac- 
teristics for the operational accuracy of a linear stationary dynamic system: 


S,@) = 1W (jo) @), (23) 


+00 +20 
K,(t) = Sy (o)cos@tde, Dy = S,(o)do, (24) 
—co 


where Sx(w), Sy(w) are the spectral densities of the signals at the input and output of the system, respectively. 
Finding the steady~state accuracy characteristic by means of Eq. (23), (24) is usually associated with computing 


complex integrals. 
We shall demonstrate that from relationships (12), (13) we obtain relatively simple formulas that express the 
results of the integrations indicated by (24). 


If the system is definitely stable (all the quantities Xj, \j have negative real parts) and if the preliminary 
operating time for the system is sufficiently great (the value of T is large), then on the basis of Eq. (13) the function 
f(A 4, Xj) can be expressed with a sufficient degree of accuracy by the formula 


1 t 


whose components are independent of the nonstationary quantity T. 


We shall substitute (25) into formula (12) and take into account the fact that for large values of T the correla- 
tion function of the output variable depends solely on the value of r. After simple transformationswe obtain 


r+nr 


i=1 j=1 j=1 t=] 


But, in accordance with (7), 


r+n r B 
i = 
= W (—A,) F.(—A,), j W (—A,), 
> j j > ( 


j=1 


and therefore we can write the following formula if we take (5), (6), and (8) into account: 


r 
K = 


| 
| 


A corollary of formula (27) is the formula 


D, = K,()= >} | 


oach (—%) Wa) 
proach MA) 
+ 
harec- which expresses the steady-state dispersion of the output variable and is equivalent to formula (22), Usually formula 
(28) is more complex than formula (22). 
In the particular case where white noise is applied to the system input, the function Fx(p) has no poles. Under 
(23) these conditions we obtain the following equation from formulas (24), (27): 
r 
this formula is valid when the multiplicity of all the poles of the function W(p) is equal to unity. 
. The use of Eq. (27)-(29) is possible if the position of the poles of the functions W(p) and F,(p) is known. If 
this is known, then the steady-state characteristics for the operational accuracy of the dynamic system can be found 
$ the by another method. 
Starting from formula (27), we find the unidimensional transform for the steady-state correlation function Ky(r) 
uary of the output variable. Denoting this transform by Fy(p), we shall prove the validity of the equation 
nction 
F, (p) + Fy (—p) = W (p) W (—p) IF, (p) + F, (—p)) (30) 
(25) by means of simple transformations for an arbitrary location of the poles of the functions W(p) and F,(p). Relation- 
ship (30) is a generalization of formula (23), since in accordance with the definitions of the functions F,(p), Fy(p), 
Sx(w), Sy(w), 
orrela- F, + (— jo) = (@), Fy (jo) + Py (— jo) = 2nS, (@). (31) 
Originating from the proof of formula (30), the following method for determining the steady-state correlation Ky(r) 
(26) and the steady-state dispersion Dy can be proposed: 
1) We find the transform F,(p) of the correlation function for the input signal. 
2) Werepresent the function 
—p) M (p) 
W W (— = 
in the form of a sum of simple fractions: 
Ss AY 
W (p) W (— p)F x = 
then, using the method of indeterminate coefficients, for example, we find the functions S(p) and S;(p). 
en 3) In accordance with the formula 


S(p) (p) 
Fy (?)= Rip) Np) + Rip) 


we determine the function Fy(p). 
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4) Using the known transform Fy(p) we determine the original Ky(r). 


5) Using the initial-value theorem [4], we compute the steady -state dispersion of the output variable without 
resorting to tables or to integration on the basis of the formula 


D, = K, F,(p)- 


In determining the originals and transforms we can use the tables of operational calculus formulas with the 
appropriate correction; these tables differ from the tables of Laplace transform formulas in that the transform is 
multiplied by the argument p. 

Examples 

1. White noise is applied to a system that is described by the equation 


(p? + 10p + 100) = (0). 


It is required to determine the correlation function for the signal at its output. 
In this case 


i 


(tu + wy)(u® + 10u + 100) (u? + 10u + 100) 


o,, (u, = 


In accordance with the data in Table 2, the transform ©,y(u, uy) corresponds to the correlation function 


K,, (t, th) = Ry (7) e**cos 5 + Re (7) 5 


Ry (1) = 1.6-10-* {3 — [2 — con (20 
R, (T) = 4.6-407* [vs sin (20 V3T+ 


after transformations. 
2. A nonstationary random signal with the correlation function 


K, (t,t) = D0) 


is applied to the input of an aperiodic section. 
It is required to determine the correlation function for the process at the output of the section. 
In this case 


k 
WO= 


and in accordance with the data in Table 2, 


i 
®,, (4, “)= + 


and therefore 
where we have 
840 


ke 


ke 
®,, (uu) = (—a, —a) + F (—A, — F (—a, 


Performing the transition to the region of originals by means of the formulas in Table 2, we find the function 
Ky(t, ty)- 
3. A stationary random process with the correlation function K,(r) = D,e~*7 is applied to the input of a sec- 


tion with the transfer function W(p) = 1/(p +a). It is necessary to determine the time dependence of the output - 
variable dispersion. 


In this case 


D, 


In making use of the tables for operational calculus formulas [6] with the appropriate correction, we find the 
originals x,(t), X(t): 


(t) = (t) = 


On the basis of formula (16) we obtain 
t 


0 


Assuming t + ©, we find the steady-state value for the dispersion of the output variable 
D, 
Dy = Ta’ 


4. White noise is applied to the input of a section with the transfer function W(p) = 1/(p + a)’. It is required 
to determine the steady-state dispersion of the output variable. 
In this case the function W(p) has one pole A; = -a with the multiplicity m, = 2, and the function F,(p) has 
no poles. Therefore we find 
d i i 


by making use of formula (21). 


5. We shall study the problem which was cited and partially solved in [1]: a random signal with the correla- 
tion function 


K,,(t) = D, 


is applied to the input of an oscillatory section. 


7 - 


It is necessary to determine the steady -state correlation function and the steady -state dispersion of the output 
signal. 


In this case D 
x 


and therfore Q(p)= 1, R(p)= p’ + ap + b, M(p) = Dy, N(p)= p+c. If b > a*/4, then the poles of the function W(p) 
are located at the points), —a + 42 —a — /B,» where a= — The pole of the function 
Fx(p) is located at the point As = -c. Making use of formula (27), wé obtain 

D. 


© 
Ky (t) = (@ Fac + — ac +b) te + [ + 
x (ct — a? 36) sin pr |} 


after a series of transformations. 
Assuming r = 0, we find 


(a + ¢)D, 


Dy = pac 


If it is necessary to determine only the steady-state dispersion, then it is expedient to make use of the equation 
W (p) W (—p) F,, 


D p+a a+ec 


= 


Therefore 


D 
x p+a a+e 
p+attac—p 
+ ac + b)(p* + ap 


Making use of the initial-value theorem, we find 


(2+ ¢) D, 
Dy = lim PF y(P)= ac +b) 
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An industrial remote control system with a two frequency combination code is described which 
uses a tone range of frequencies. The apparatus is designed for the remote control of widely 
distributed groups of objects and provides transmission and reception for a large number of 
signals. There is described a contactless converter of discrete values of rotational angle into 


a frequency combination code which is based on scanning the indications by a static induction 
method. 


The SRP-3 remote control system for lumped objects was developed by IAT AN SSSR jointly with VNIIKAneftegaz 
and served as the basis for a number of the following developments. 


The transmission of signals in the SRP-3 system is accomplished by means of a two frequency combination 
code using a tone frequency range which ensures its universality with regard to communication channels,i.e., the 
feasibility of operating on physically separate wire lines, frequency multiplexed channels, and radio relay lines. At 
every control point (CP) there is provision for the remote control of objects, telesignalling their condition, telemetry, 
and telephone communication. From a dispatching point (DP) to every CP up to 24 signals can be transmitted, and 
in the reverse direction 36 signals for information and emergency alarms, and 5 of telemetry. There is provision for 


connecting to oe DP up to 10 CPs, so that the total capacity of the system in numbers of signals is 240 issuing from 
a DP and 410 signals coming in to a DP. 


The number of frequencies n employed for the transmission of N signals by means of a two frequency code with 
sequential transmission of frequencies is determined from the expression 


n =0.5(1 V4N + 1). 


SRP-3 apparatus was developed for the remote control of interconnected pumping stations in a system for 
maintaining the stratum pressure in the petroleum industry [1], but it can be used for the remote control of any 
groups of objects. In the course of industrial operation of an experimental assembly, some changes, which are de- 
scribed in the present article, were introduced. 


A block diagram of one-terminal sets for a dispatching and a control point is shown in Fig. 1. There are 
similar sections in the DP and CP half assemblies: the amplifiers of the receive and transmit channels, the coders, 
the decoders, and the control sections. 


The normal state of all one-terminal sets is a stand-by condition for reception. When it is necessary to transmit 
a signal, the control section switches amplifiers from receive to transmit. At a DP the control section is switched 
on from the selector section, and at a CP when one of the information generators is operated. The signal is encoded 
by a coder which causes voltages of appropriate frequencies to be generated. The sequence and duration of the trans- 
mission is determined by the control section. 


A signal received at a CP is fed, after amplification, to a decoder, on the output of which control devices are 
connected. 


The transmission of information signals and emergency alarms proceeds in a similar fashion during transmissions 
from information generators. 


=| 
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During a telemetry call the control section makes a measurement, switches the amplifier to transmit, and 
connects a measuring generator (for the frequency method of measurement). 


A CP coder is shown in Fig. 2a. The networks of the coder control are connected by the information generato, 
contacts. While a signal is passing, voltage is removed from the control networks, and that, in conjunction with the 
capacitor connected to the relay and their self-blocking, delays a new signal until the completion of the preceding 
one. Contacts of the A and B relays (a single relay operates simultaneously in each group) connect networks to the 
RC oscillators (Fig. 2b) thus creating voltages of appropriate frequencies. 


Fig. 1. Block diagram of one-terminal sets at dispatching and control points 
of the SRP-3 system. SS) Selector section; C) coder; CS) control section; TA) 
transmit channel amplifier; RA) receive channel amplifier; D) decoder; SCS) 
signal control section; MS) measuring section; CD) controlling device; OC) 
object controlled; MCS) measuring control section; TSG) telesignalling infor - 
mation generator; TG) telemetry information generator; MG) measuring in- 
formation generator. 


L 


L 
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Fig. 2. Transmission sections for signals from CP to DP, a) Coder; b) oscil - 
lators; c) control section. 


While relays of groups A and B are operated, the control section is connected and performs the following func- 
tions. 


Relay TR, cuts off the voltage from the coder control networks and operates relay TR,,which, in turn, connects 
the amplifier of the transmit channel. Relay TR, carries out the sequential transmission of the oscillator frequencies. 
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It is possible to use a single oscillator having the networks connected to it by turns as determined by the contacts of 
relays A and B. The TRr relay limits the total time of a signal transmission and restores the circuit to its original 
state. 


The decoder consists of a series of frequency selectors, equal in number to the number of frequencies in the 
encoded signals which are accepted by the particular decoder. Series resonant LC networks are used as frequency 
selectors, and the possibility of using K-type filters (Fig. 3) is also provided for. The passband of the selectors is 
150 to 170 cps; the attenuation for frequencies in the cutoff region is 4 to 5 nepers. 


Tr to Try 


Fig. 4. Receiving sections for signals at a CP. a) Contact decoder; b) contactless 
decoder; c) interference protection circuit. 


The decoder networks were made in two types; contact and contactless. The contact type (Fig. 4a) is dis- 
tinguished by the small number of contacts. The output of each frequency selector consists of a relay with two 
make contacts. When two selectors are operated, a circuit is completed to the appropriate control device. 


The contactless type (Fig. 4b) differs from the contact type in that transistors, biased to cutoff, replace the 
telays. Voltages from frequency-selector outputs switch on the transistors so that common paths are produced. 


In order to increase the noise immunity of the system there is added to the signal combination a lead fre- 
quency which precedes the transmission of code frequencies. The lead frequency pulse passes into a selector ahead 
of the decoder which has relay F,, for an output. Contacts on F, turn on a relaxation oscillator (Fig. 4c). The in- 
coming pulse is integrated by network RyC;,whose time constant ensures that the interference protection relay IPR 
operates only on pulses of a specified duration which is appreciably longer than the duration of the expected inter- 
fering pulses. During either one or a series of short interfering pulses,the circuit reverts after each pulse to the initial 
condition by discharging the capacitor through the small resistor Rg. Hence interfering pulses do not accumulate and 
cannot cause the apparatus to operate. 


Upon the arrival of a working pulse of appropriate duration relay IPR operates, locks itself up, and connects, a 
battery to the output circuits of the decoder. 
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At the expiration of the time required to receive the signals,relay IPR reverts to the initial condition. 


The dependability of the interference protection arrangement was tested with various kinds of real and ar- 
tificial interference, in particular during operation over radio channels. The tests showed extremely high reliability 


of operation for the arrangement. 
For telemetry transmission no additional control or coding sections are necessary. Provision is made in the 
SRP-3 system for frequency, time pulse, and code telemetry (frequency combination code). 


For the numerical registration of the parameters measured in VNI[KAneftegaz,a numerical converter arrange- 
ment was developed which has an output to an electric typewriter with solenoid controls, thus providing printed 
results in numerical form for the three telemetry methods mentioned. 


For the measurement of discrete steps of angular rotation in IAT AN SSSR and VNIIKAneftegaz an information 
generator was developed in which the value of a rotation angle is transformed into a frequency combination code by 
the static induction method. When using this information generator, the telemetry signals are no different from the 
transmission of control or information signals. 


The information generator consists of several stages, the number of which determines the accuracy of the 
measurement, In each stage there is a counting wheel of nonmagnetic material with magnetically permeable strips 
distributed around the periphery of the wheel. Over each stage there is located a contactless scanner comprising 
five magnetic circuits with windings which are completed by the magnetically permeable strips of the counting 
wheel. The positions of these strips causes the completion of two magnetic circuits out of five for every position 
of the counting wheel, which thus produce a code for a specific number. 


A secondary winding on each magnetic circuit is connected to a transistor controlling a fixed frequency oscil- 
lator. The scanners of all the stages are connected to the five oscillators by turns. 


In the information generator the counting wheel of the first stage is brought up to an integral number of units. 


The SRP-3 remote control apparatus was introduced into the interconnected pumping stations which maintain 
the stratum pressure in the petroleum industry of the Tatar Council of National Economy. The remotely controlled 


' interconnected pumping stations have been converted to “locked up” operation without the constant presence of duty 
personnel. 


Based on the SRP-3 apparatus, remote control systems of other objects have been developed for future use. 
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The paper studies the principles of designing dispatcher monitoring systems which make use of digital 
techniques for teletransmission and reproduction of information; the case where information is processed at the dis- 
patcher point is also treated. 


At the present time the transmission of data to a dispatcher point and data processing at the dispatcher point 
are performed almost exclusively on an analog basis. However, the modern state of the techniques available for 
digital conversion and reproduction provides wide avenues for perfecting methods for transmission, reception, re- 
production and processing of information in dispatcher monitoring systems. 


The most rational engineering solution is achieved if the design of a dispatcher monitoring system is sub- 
ordinated to the main problem, which is to reduce the volume of information on the technological process to the 
required minimum. Such a proposition requires 4 certain clarification. 


At the present time the generally accepted tendency toward expanding the sphere of monitoring and control 
which is encompassed by a dispatcher point has led to an increase in the complexity of the functions which are in- 
volved in dispatcher monitoring and control. Therfore the necessity of a special selection of the essential data, and 
their automatic classification and processing, is understandable. In the final analysis such a development logically 
leads to the application of control machines which exclude the intermediate participation of a human operator be- 
tween the monitoring system and the data processing on the one hand, and between the elements which control the 
production processes on the other hand. At the present stage of development of monitoring systems, it is still easy 
to speak of eliminating the dispatcher; therefore as the volume of input information increases, the requirements 
imposed on the monitoring system involve an ever greater degree of processing for this information [1, 2]. One of 
the results of processing input data by the system is the operative information which is used by the operator to affect 
the process directly. It is possible to indicate the following stages of development in formulating the functions 
which govern the preparation of the operative information. 


1. The dispatcher performs simple processing of the received information, estimates the the situation, and 
adopts operative decisions directly on the basis of the readings of the telemetering instruments and telephone com - 
munications. 


2. To the degree that the volume of the monitored production process increases, the continuous readings of 
a large number of instruments are replaced by signals with information on the output of the monitored parameters 
within the allowed limits. Such a system not only economizes monitoring equipment but also solves the important 
problem of filtering the information. 


3. When signalling a small number of limiting values at each monitored point proves to be inadequate, 
automatic computing devices are used for processing the data. Using computers, generallized indices of the operating 
mode of the monitored units are determined and reproduced on the output indicating meters. Thus the dispatcher is 
freed from performing the operative computing operations, and at the same time the quantity of data monitored by 
the dispatcher is reduced. 


4. A further automatization of the information processing permits the formulation of the control demands to 
be automatized also. Here the direct participation of a human operator in the operative control of the process is 
completely or partially eliminated, 
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Besides the operative information which allows a dispatcher to adopt immediate decisions on direct manip- 
ulation of the process, it is usually necessary to acquire statistical information on the process as well. In contrast 
to operative information, statistical information is subjected to averaging or comparison over a definite time interyaj 


(a day, a month, etc.) and must be retained for a prolonged period; it always requires the use of a memory device in 
which it is recorded. 


7 


Digital recording can be used both in the case where the recorded information is designed for subsequent 
readout by a human operator and in the case where it is designed for machine readout. In both cases the digital 
form of recording has appreciable advantages. Recording in tabular form is somewhat less convenient than analog 
diagrams obtained by means of automatic recorders, but it permits us to achieve an appreciably greater accuracy in 
terms of both the actual measurement and the synchronization of measurements at a large number of points. 


Moreover, when compact recording of many parameters is required, digital printing on one blank is definitely 
the most practical method. For further processing of statistical data in digital computers, only the digital form of 
recording is a practical method (on punch cards, punched tape, or magnetic tape). 

Digital Reproduction of Readings 

For purposes of obtaining information in the form of instrument readings both analog and digital indicators are 
equally practical. However, analog instruments are in themselves considerably simpler. The advantage is with 
digital indicators only in those cases when an appreciable noise level in communication channels makes it neces- 
sary to transmit information in coded form [8, 4]. Pulse-telemetering systems, as we know, are among those with 
the highest noise stability. If we also take into account the fact that the noise stability can be appreciably in- 
creased by using codes with a high noise immunity, then the advantages of pulse -code systems over telemetering 
systems with an analog signal is obvious in those cases when noise immunity is the most important property of the 


system. The comparatively high cost of pulse-code telemetering equipment in one metering channel is lowered 
when we design a multichannel system with a common coding unit. 


When a digital code arrives from a communication channel, the conversion to analog form at the receiving 
end is, as a rule, completely unjustified. Such a conversion leads to the loss of a number of advantages which a pulse 
_ code telemetering system has; these include the absence of receiver errors, objectivity of readings, and the presence 

of memory. 


We should especially note the memory of digital instruments. In all those cases where information arrives at 
the instrument at certain time intervals, it is necessary to memorize the readings over a comparatively short period. 
When the duration of the interval is to any extent appreciable, the solution of this problem by analog techniques is 
extremely difficult. Only one practical method for memorizing readings in analog form is known—the use of a 
system with a motor or with a mechanical means for fixing the indicator. However, data in digital form can be 
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memorized using the simplest means: practically any elements which have two distinguishable states (relays, triggers, 
etc.) can be used. 


The necessity for memorizing the readings in a multichannel telemetering system with time division of the 
channels may in itself require conversion to the digital form of reproducing readings if the monitored object is inertial 
and the cycle of measurements is greater than the time constant of the reproducing instrument [1]. 


The advantages of digital indicators also include the fact that they yield a reading that is easily reproducible 
by the operator for very small dimensions of the indicators proper; this is due to the rational use of the indicating 
phase area of such instruments. Thus it is possible to mount such instruments directly in thememo-circuit of the 
monitored object. Moreover, digital reproduction eliminates subjective errors in taking readings. At present, 
satisfactory designs for digital indicators are already available (projection indicators, indicators with multielectrode 
gas tubes, etc.), and work on perfecting them is continuing [5, 6]. 


The choice of the repetition periodicity of the digital measurements for each parameter is determined by the 
upper limit of its frequency spectrum. In accordance with the Kotel'nikov theorem, the behavior of the parameter 
with respect to time can be established in single-valued fashion from the values of individual measurments that are 
repeated at a frequency which exceeds the upper boundary of the frequency spectrum for the parameter by more than 
a factor of two. In the case of digital data processing in which the pulse train is not passed through the filter, this 
requires multiplication of each ordinate by a function of the form sin at/at and addition of the products. Here it 
is possible to determine the behavior of the measured quantity exactly in the interval between the kth and (k + 1)th 
measurements only after processing the results of a number of subsequent measurements: from the (k + 1)th to the 
(k + n)th. The number n is the greater, the higher the required accuracy of computation. Therefore in order to clarify 
the behavior of the measured parameter with respect to time from the values of individual measurements taken in 
accordance with the Kotel'nikov theorem, we require both cumbersome computations and a time delay. 


If a time delay is impermissible and computations are undesirable, it is necessary to choose a higher quantiza - 
tion frequency than the theorem requires; then horizontal steps are used to approximate the behavior of the function 
in the intervals between measurements. Under these conditions the relationship between the quantization frequency 
fq and the upper boundary fh, for the frequency spectrum of the function depends on the allowable approximation 
error. Figure 1 shows the analytically derived dependence of the mean-squaze error § of the step approximation 
of the frequency ratio fg/f, for the case where the function is random, stationary, and has a uniform energy spec- 
trum from 0 to fh. 


It can be demonstrated that a digital converter can be designed for a number of quantization levels such that 
the mean-square level quantization error will be of the same order of magnitude as the error of approximating 
the dependence of the parameter on time. This makes the choice of the number of quantization levels dependent 
on the time quantization frequency. 


The transition to digital reproduction of the parameters is possible in dispatcher systems for inertial industrial 
objects whose cutoff frequencies are sufficiently small. It is expedient to establish an upper limit for the average 
frequency with which the readings of the digital indicator change on the basis of the possibilities of human perception. 
Thus far there is insufficient experimental data, but on the basis of preliminary results it is possible to assume that 
this frequency is equal to 0.2 cps. Having specified this value, it is possible to use the function shown in Fig. 1 to 
determine the digital reproduction error for any value of the upper boundary for the spectrum of the monitored 
process. The corresponding relationship is shown in Fig. 2. 


In those cases where it is sufficient to provide the operator with averaged readings of the process flow during 
comparatively long time intervals, it is possible to use digital reproduction devices which incorporate intermediate 
conversion to a frequency followed by calculation of the number of oscillations during a fixed time interval. 


Returning to the assertion that it is expedient to perform digital conversion of signals received at the dis- 
patcher point from multichannel analog telemetering systems with time division of the channels, we should in- 
dicate the comparatively simple engineering realization of this problem. The signal from a frequency or pulse- 
time telemetering system is in itself already prepared for digital conversion. In fact, it is sufficient to compute the 
number of oscillations of a frequency telemetering signal over a time interval, or the number of oscillations of a 
fixed-frequency oscillator over the duration of the signal from a pulse-time telemetering system in order to ob- 
tain the numerical value of the parameter. Of course, here it is necessary to make a proper choice of the fre- 
quencies and the time required for a reading. The structure of the scaling circuit is determined by the nature of 
the chosen code. 
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Fig. 3. A circuit for digital conversion with the introduction of scale coefficients by means of volt- 
age dividers. T) Transducer; VD) voltage divider; C,,C,) commutators; DC) digital converter; Dec.) 


decoder; DI) digital indicator. 
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Fig. 4. Circuit for digital conversion with the introduction of scale coefficients by shifting the 
places of the binary code (the notation is the same as in Fig. 3). 
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We shall now study certain problems involved in digital conversion for the purposes of both indication and 
recording. 


Since the dispatcher must obtain the digital value of the parameters in a decimal system, it follows that the 
scaling circuit must be designed either in accordance with the decimal system or in accordance with a system which 
allows subsequent single conversion of the resulting code to a decimal code. In our opinion this problem can be 
solved most economically by using a binary-decimal scaling circuit with a simple matrix decoder which converts 
four binary stages in each decade of the scaling circuit into one of the digits in a decimal number. The digital re- 
production is achieved most simply in percentage form, since for various scales corresponding to the variation of the 
measured parameters,all of these scales are reduced to one nominal at the input of the digital converter. In that 
case the resulting digit expresses the relative value of the parameter. In individual cases the transition to a percentage 
scale is permissible. However for the most part we encounter the requirement of reproducing the values of parameters 
in an absolute digital expression. Under conditions where the dispatcher point does not contain computer devices 
which permit us to multiply the number by an individual coefficient for each scale, we are left with two possibilities. 


It is possible to introduce the scaling coefficient even before introducing the analog parameter into the digital 
converter. This can be done conveniently using a simple voltage divider if the parameter is represented in the form 
of a voltage. In that case we obtain a definite relative quantization error for digital conversion. Therefore the 
digital converter must be designed in such a way that for a minimum nominal for the input analog quantity the 
quantization error does not exceed the specified value. Under these conditions, increased requirements with respect 
to the number of digital places are imposed on the digital converter as compared to the case of reproduction in 
percentage form. 


We shall clarify the above on the basis of the example and illustration cited in Fig. 3. Assume that the 
nominals for the measured parameters have various magnitudes: from 25 to 100 units. If the quantization error for 
digital reproduction must not exceed +0.5%, then in order to produce the appropriate numbers in percentages we 
require a coding device incorporating two decimal stages; for producing the appropriate numbers and absolute values 
we require a digital indicator with three decimal stages for a scale of 25 units. Therefore the same number of stages 
must be contained in an over-all digital converter. 


Of course, if we allow a lower accuracy for the reproduction of the absolute numbers when lower nominals are 
used, then the digital converter does not become more complicated. 


The other possibility for using absolute digital readings without making use of arithmetic units at the output 
can be simply realized when the digital converter produces a binary code. Shifting a binary code n places to the 
left is equivalent to multiplication by 2". This makes it possible to introduce coefficients which are multiples of 
the numbers 2, 4, 8, 16,etc. after digital transformation. Under these conditions the relative magnitude of the 
quantization error remains constant, Since the nominal for the absolute scales do not always differ by a factor of 
2", it follows that an additional correction can be performed using a divider at the input analog end of the system. 


We shall study the example illustrated by Fig. 4. Assume that we have a binary coding device with seven 
stages. The number of quantization levels in this device is equal to 128. If the range of variation for the analog 
parameter at the converter input is chosen within the range of from 0 to 100 quantization levels, it follows that the 
coding will be performed with a normalized quantization error equal to + 0.5%, Assume that these 100 conditional 
units of the analog parameter correspond to an absolute nominal of a hundred units for measured parameter No. 1, 
an absolute nominal of 400 units for measured parameter No. 2, and an absolute nominal of 750 units for parameter 
No, 3. In all cases the coding device produces the binary code corresponding to the number 100 for the nominal 
value of the parameter; this code is equal to 1100100. For the first parameter this code must be retained, and for 
the second and third parameters it must be altered. In order to multiply this number by 4 (for parameter 
No. 2) it is sufficient to shift the entire code two places to the left for transfer of the code to the cells of the digital 
memory for parameter No. 2. We obtain the binary number 110010000 which is equal to 400. For this purpose we 
need only two extra memory cells. Increasing the number of code places to 9 does not mean 4 variation in the 
relative quantization error, since the number will be obtained with an accuracy of up to the third place counting 
from the lowest place. 


As far as parameter No. 3 is concerned, a voltage divider with a division coefficient 750/800 = 15/16 must 
be placed after the corresponding transducer. Then the nominal value of the measured parameter will correspond 
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to (15/16) - 100 conditional units of the analog parameter at the input to the digital converter; with an accuracy of 
up to integers this will be equal to 94 units, or 1011110 in the binary code. Making a shift of three places to the 
left for transfer of the code to the memory cells (i.e., multiplying by eight), we obtain 1011110000, This corre- 
sponds to the number 752. The accuracy with which parameter No. 3 is reproduced will be slightly lower than the 
accuracy in the first two cases. 


For scales with a suppressed zero the introduction of scale coefficients is extremely complicated, since in this 
case it is necessary to subtract the particular number corresponding to the parameter in question from the coding 
result. 


Summing up what was said above on the relationship between analog and digital techniques in performing the 
simplest function of a dispatcher monitoring system (i.e., for the simple reproduction of data), we should note that 
digital techniques are superior for transmitting information over noisy channels, reproducing signals in multichannel 
low-frequency telemetering system with time division, and in all other cases where a memory of long duration is 
required. 


Digital reproduction is also convenient to use at a dispatcher point for individual parameters which are trans- 
mitted in analog form. In the case under study the convenience of such a solution as far as the dispatcher is con- 
cerned is obvious, since it is easier for him to make use of a single mnemo-circuit with digital indicators and a single 
recording blank. Furthermore, the converse is also probable: For individual parameters whose spectrum contains 

comparatively high frequencies, it is necessary to use point indicators irrespective of the method by which the signal 
is transmitted to the dispatcher point. Moreover, if it is important to make it possible to achieve a convenient 

estimate of the flow of the process over a certain time interval, then analog automatic recorders which trace out a 


graph are incorporated. 


The Equipment for Centralized Monitoring and Processing of Data 

In the analysis of a dispatcher monitoring system which we performed above,we did not take into account 
the problem of filtering the information. Only when this system is augmented by equipment which has become 
known a: centralized monitoring equipment does it become fully usable from the point of view of the minimum- 
information criterion. 


The equipment for centralized monitoring can be designed on the principle of comparing parameters with 
settings both in digital and analog form. The latter method is realized considerably more easily, but it is applicable 
only for transmission of the parameters to the dispatcher point in analog form. 


The design of centralized monitoring systems using the digital form evidently makes it possible to obtain a 
universal dispatcher monitoring system which is simpler in design and more successful from the point of view of 


equipment usage. 


When the received information requires more complex processing associated with the performance of mathe- 
matical operations, the advantages of complete conversion of dispatcher monitoring systems to the digital basis are 
even more clearly manifested. Only digital computers have the necessary flexibility and accuracy, and therefore 
it is precisely these devices which must be used in universal data processing systems. It is logical to design the 
entire dispatcher monitoring system in such a way that the output data of the monitoring system is introduced directly 
into a computer without any intermediate conversions. 


Thus, digital conversion at the earliest possible stages of information transmission is rational here. This 
principle is employed most appropriately for transmission of signals over a communication channe! in coded form 
(i.e., in a pulse-code telemetering system). In those cases when the use of comparatively complex and extensive 
coding equipment at the transmitting end is not justified by cost-engineering considerations (for single-channel 
systems), the analog telemetering signals which are received at the dispatcher point should be subjected to digital 
conversion. Such a converter may also service local measurements. 


The Transmission of Commands 

We shall study methods for transmitting commands from the dispatcher point. Our use of quantized signals for 
this purpose is based on reliability consideration. Thus it is true that the means for transmitting the control signals 
will also be basically identical to those for digital transmission of information from the objects to the dispatche. 
This provides broad possibilities for the unification of the engineering means for measurement and control, and for 


the realization of the block principle in designing equipment for dispatcher systems. 


We should make a separate analysis for the case where smooth remote control signals are transmitted (the 
remote control process). Here we can introduce the following additional arguments in favor of the coded form in 
choosing between transmission in coded form or analog form. 


1. Even for manual transmission of the control pulses (for which there is no direct coupling between the digital 


system for processing the information and the actuating elements) the use of coding is justified by a considerable 
increase in the noise immunity of the commands. 


2. For automatization of the command function on the basis of digital processing of the information the reverse 
transition to the analog form of the signal in the final stages of the system is not rational. Therefore it is more ex- 
pedient to send quantized control signals over the channel. If under these conditions we make use of actuating 
mechanisms that are controlled by discrete commands (for example, step motors), then digital methods will become 
most effective. 


The argumentation presented above in favor of introducing digital techniques into all the phases of reception, 
processing, and transmission of information from the dispatcher point makes it possible to conclude that this is the 
most rational path for achieving the transition to complete automatization based on the use of control machines. 
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